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1. Introduction 



In generic analytic nonlinear differential systems in the complex plane, 
we study the position and the type of singularities formed by solu- 
tions when an irregular singular point of the system is approached 
, along an antistokes direction 1 . Placing the singularity of the system 

{Sj ' at infinity we look at equations of the form y' = f (x , y) with f ana- 

lytic in a neighborhood of (0, 0), with genericity assumptions; x = oo 
is then a rank one singular point. We analyze the singularities of 
those solutions y(x) which tend to zero for x — * oo in some sectorial 
4^ ' region, on the edges of the maximal region (also described) with this 

property. 

After standard normalization of the differential system, it is shown 
that singularities occuring in antistokes directions are grouped in 
nearly periodical arrays of similar singularities as x — » oo, the loca- 
tion of the array depending on the solution while the (near-) period 
and type of singularity are determined by the form of the differential 
system. 

This regularity in type and position of movable singularities has 
been observed previously in various examples of nonlinear systems: 
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1 In the sense stemming from Stokes original papers and the one favored in ex- 
ponential asymptotics literature, Stokes lines are those where a small exponential 
is purely real; on an antistokes line the exponential becomes purely oscillatory. In 
some references these definitions are interchanged. 
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Painleve equations ((BUI, ESj- ^El) third order nonlinear equations 
(HH, |22j) nonintegrable Abel equations ([22], EH) 

among others. 

We show these features are rather universal and find a formalism to 
calculate them (asymptotically). 

When f is meromorphic and satisfies some estimates the singular- 
ities in the arrays are generically square root branch points. 

The mechanism of singularity formation is elucidated by expo- 
nential asymptotic analysis, which also provides a general and effec- 
tive calculation tool for determining the type and position of singu- 
larities. The present method generalizes that of ^Bl- The analysis 
yields two-scale asymptotic expansions of solutions, valid in a re- 
gion which includes the directions along which y — ► and extend- 
ing, on appropriate Riemann surfaces, into regions where the solu- 
tions typically develop singularities. The expansions have the form 
y ~ F(s;£(aO) = E^o^^m^W), where £(x) = Ce~ Xx x a - A, a 
and the functions F m are uniquely determined, modulo trivial trans- 
formations, by f; the constant C depends on the solution y. F m 
satisfy a recursive system of equations, typically simpler than the 
original system. In particular, for all first order equations and for 
Painleve's PI and P2 equations, the solution of the recursive system 
is expressible by quadratures. 

The method can be interpreted as a trans asymptotic matching 
technique in that the expansion F of y matches (and is fully deter- 
mined by) its Ecalle transseries 2 in a sector where y is asymptotic 
to a power series. The constant C in the definition of £ is one of the 
constants beyond all orders in the transseries of y. In some instances, 
the technique provides a connection method even in nonintegrable 
systems (in which case, the connection data are path-dependent). 
The constant C becomes thus accessible by classical asymptotics and 
determines the position of singularities of y. 

The expansion y ~ F(x;S t (x)) satisfies Gevrey-type inequalities, 
and thus produces exponentially accurate estimates of y (see 

Some examples are outlined. In the first one, a nonintegrable Abel 
equation, the method provides a description of the exact type of all 
but finitely many singularities of solutions in a sector, and of the as- 
sociated the Riemann surfaces. The connection between these com- 
plicated Riemann surfaces and the numerically observed chaoticity of 
solutions [2S1 is briefly discussed. 

As other examples we consider the Painleve equations Pi and Pn 
for which we use the technique to express the asymptotic distribution 
of poles near the antistokes lines of the so called truncated solutions. 



2 In our context these are algebraically determinable formal combinations of se- 
ries in x~ x and small exponentials, and generalize classical asymptotic expansions 
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The position of the poles only depends on an exponential asymptotics 
quantity, the constant beyond all orders. 

2. Setting 

We adopt, with few exceptions that we mention, the same conditions, 
notations and terminology as ^Hj; the results on formal solutions and 
their generalized Borel summability are also taken from |13j . 
The differential system considered has the form 

y' = f(x-\y) yeC\ xgC (1) 

where 

(i) f is analytic in a neighborhood V x x V y of (0, 0), under the gener- 
icity conditions that: 

(ii) the eigenvalues Xj of the matrix A = — < f^-(0, 0) > are 

I % J i,j=l,2,...n 

linearly independent over Z (in particular Xj ^ 0) and such that 
(hi) argAj are all different. 

(In fact somewhat less restrictive conditions are used, namely 
those of HI §1.1.2.) 

By elementary changes of variables, the system can be brought 
to the normalized form ^5] , [IS] 

y' = -iy + iiy + g(x-\y) (2) 

where A = diag{Aj}, A = diag{<x,-} are constant matrices, g is ana- 
lytic at (0, 0) and g(x _1 , y) = 0(x~ 2 )+0(\y\ 2 ) as x — > oo and y — > 0. 
Performing a further transformation of the type y i— > y — a k x ~ k 
(which takes out M terms of the formal asymptotic series solutions 
of the equation) , makes 

g(|x|-\y) = 0(x- M - 1 ;|y| 2 ;|x- 2 y|) (x -+ oo; y - 0) 

where 

M > max3f?(aj) 

3 

and 0(a; b; c) means (at most) of the order of the largest among a, b, c. 

Our analysis applies to solutions y(x) such that y(x) — > as 
x — ► oo along some arbitrary direction d = {x £ C : arg(x) = eft}. A 
movable singularity of y(x) is a point x E C with G where 
y(x) is not analytic. The point at infinity is an irregular singular 
point of rank 1; it is a fixed singular point of the system since, after 
the substitution x = z _1 the r.h.s of the transformed system, ^| = 
— z~ 2 i{z,y) has, under the given assumptions, a pole at z = 0. 
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2.1. Classical versus exponential asymptotics 

In order to understand the properties of solutions of (j2j) for large 
x, one way is to find formal asymptotic solutions, then use asymp- 
toticity relations to deduce information about the true solutions from 
the formal ones. It is easy to see that equation © admits a unique 
asymptotic formal power series solution |51| 

oo 
r=2 

The coefficients yo, r of yo can be computed recursively by substitu- 
tion in (j2J) and identification of the coefficients of x~ r ; the series yo 
is a formal solution, and is usually divergent. Its Borel summability 
was shown, in a more general setting, by Braaksma [§]. 

Given an open sector of the complex x-plane, of angle less than ir, 
there exists a true solution of ^ which is asymptotic to © in that 
sector (as |x| — ► oo) [HJ- This solution is not unique in general. 

To illustrate the way different solutions with the same asymptotic 
series (in a sector) can be distinguished consider the simple linear 
equation f'(x) = —f(x) + x _1 with the general solution f(x; C) = 
e~ x Ei(x) + Ce~ x where Ei(x) = P f*^ t~ x e l dt. Any solution f(x; C) 
has the same (divergent) power series asymptotic expansion in the 

right-half plane: f(x;C) ~ fo(x) = Yl r >o r] - x ~ r l ^ or x ~~ * °°> > 
0. The parameter C which distinguishes different solutions multiplies 
the term e~ x which is much smaller than all the terms of the asymp- 
totic series /o : C is a constant beyond all orders. 

The theory of linear equations with an irregular singular point is 
well developed and there are comprehensive results; we mention the 
works of Babbit and Varadarajan 1 , Balser, Braaksma, Jurkat, Lutz, 
0, [H], |5) EH) Balser, Braaksma, Ramis and Sibuya [Hj, Deligne [T7| . 
Jurkat [SJ , Katz [HS] , Levelt [SHI , Levelt and Van den Essen j37j , Lutz 
and Schafke gU, Manin [HH], Olver 00], Malgrange EH, Ramis gl], 
Ramis and Martinet 13 , Ramis and Sibuya jH] , Sibuya |46j , Turritin 
|49j . and others — see also (HOI and the references therein. 

For linear equations there exist fundamental systems of solutions 
in terms of which one can speak of exponentially small terms. A 
formal analogue for nonlinear equations is represented by formal ex- 
ponential series. 

An n-parameter formal solution of @ (under the assumptions 
mentioned) as a combination of powers and exponentials is found in 
the form 

y(x) = Yl CV Akl 2; Qk s k (i) (4) 

ke(NU{0})" 
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where §k are (usually factorially divergent) formal power series: so = 
yo (see (J3J) and in general 

oo ~ 
r=0 

that can be determined by formal substitution of (J3J) in (J2J); C € C n 
is a vector of parameters 3 (we use the notations C k = YYj=i ; 
A = (Ai, A n ), a = (a 1} a n ), |k| = fcj. + ... + k n ). 

Note the structure of Q : an infinite sum of (generically) divergent 
series multiplying exponentials. They are called formal exponential 
power series [5T] . 

Formal solutions (Jlj) of differential equations ((21) were introduced 
by Fabry [2^ and studied extensively by Cope [TT] . 

From the point of view of correspondence of these formal solutions 
to actual solutions it was recognized that not all expansions (J1J should 
be considered meaningful; also they are defined relative to a sector 
(or a direction). 

Given a direction d in the complex x-plane the transseries (on 
d), introduced by Ecalle ^HJ, are, in our context, those exponen- 
tial series (Jlj) which are formally asymptotic on d, i.e. the terms 
c k e -A-kx x a,k x -r ( with k G ( N u { }) n , r € N U {0}) form a well 
ordered set with respect to 3> on d (see also |H3)- (For example, this 
is the case when the terms of the formal expansion become (much) 
smaller when k becomes larger.) 

For linear systems any exponential power series solution is also 
a transseries: it consists of n power series multiplying exponentials 
since §k = for |k| > 2. 

In the nonlinear case if a formal exponential power series (Jlj) sat- 
isfies the condition Cj = if e~ x i x -f* as x — > oo, x € d then Q is 
a transseries on d. In fact, it is clear that (jlj) is a transseries on (any 
direction of) the open sector St ra ns defined by 

Straus = {x € C ; if Cj^ then X(\ jX ) > , j = 1, n } (6) 

This sector may be empty; it may be the whole C if all Cj = 0; 
otherwise it lies between two antistokes lines, and has opening at 
most TT. 

If yo is divergent (which is generic) then the terms containing 
exponentials in (j3J (i.e. terms with |k| > 1) are much smaller than 

3 In the general case when some assumptions made here do not hold, the general 
formal solution may also involve compositions of exponentials, logs and powers 
|19|. The present paper only discusses equations in the setting explained at the 
beginning of S|3 

4 We note here a slight difference between our transseries and those of Ecalle, 
in that we are allowing complex constants. 
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all powers of x in yo and cannot be denned by classical asymptotic 
inequalities in the Poincare sense. Hence their designation: terms be- 
yond all orders. 

Transseries and their correspondence with functions are the sub- 
ject of exponential asymptotics, which developed substantially in 
the eighties with the work of M. Berry (hyperasymptotics), J. Ecalle 
(the theory of analyzable functions), and M. Kruskal (theory of tower 
representations and nice functions) (see references). 

From a historical point of view we must stress the importance 
of the fundamental work of Iwano. Generalizing earlier results of 
Malmquist he proved in wide generality that locally meromorphic 
systems of differential equations have expansions, which for the class 
(J2J) discussed in the present paper have the form 

y(s) = 0o(*) + E C k e- A - k ^ a ' k k (x) (7) 

keN" 

convergent, and with 0k analytic, in appropriate sectors [2S], |27| - 

Later, Ecalle introduced a very large space of expansions, relevant 
to differential, difference, integral and other equations. The funda- 
mental space where formal solutions are sought is purely algebraic — 
the transseries (for example, see the expansion © for solutions 
of ©)• Then a general procedure (based on Borel summation), inde- 
pendent of the equation where the expansions originate, is outlined 
to associate functions to formal expansions. As a consequence the 
summation procedure can be used in a broad class of problems and 
yields a complete isomorphism between formal expansions and a class 
of functions (analyzable functions). For differential equations this pro- 
cedure shows that all the terms 0k in (jZj) can be in fact obtained from 
00 by a form of analytic continuation (Ecalle's resurgence relations). 
Also, Stokes phenomena can be described in detail. The paper [P3*] 
proves this procedure in the context of differential systems with a 
rank 1 singularity. 

Braaksma has recently extended the theory to nonlinear difference 
equations [TO] . 

The present paper studies the solutions in a region where the ex- 
pansions ® and those of .13; diverge. It is shown that in this region 
the solutions of generic systems actually do have singularities (see 
Theorems El and |3J) grouped in regular arrays. Thus, a posteriori we 
know that expansions of the form (jJJ) cannot converge there. Never- 
theless, asymptotic representations, in terms of functions themselves 
singular derived in the present paper hold in this region (see l)26jl ) and 
they enable finding the singularities of y. The region where the se- 
ries (ffi))) is asymptotic to solutions and the region where (J7J) , or (jTU)) , 
converge do intersect so the two expansions of the same solution can 
be matched in this region (see Theorem 
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2.2. Further notations and results referred to 

We recall that the antistokes lines of ([2| are the 2n directions of the 
x-plane iAjR + , —i\j~B& + , j = l,...n, i.e. the directions along which 
some exponential e~^ x of the general formal solution (jlj) is purely 
oscillatory. 

In the context of differential systems with an irregular singular 
point, asymptoticity should be (generically) discussed relative to a 
direction towards the singular point; in fact, under the present as- 
sumptions (of nondegeneracy) asymptoticity can be defined on sec- 
tors. 

A first question is to determine which are the solutions asymptotic 
to the power series solution @, and to find their regularity. 

Let d be a direction in the x-plane which is not an antistokes line. 
The solutions y(x) of ((2J) which satisfy 



are analytic for large x in a sector containing d, between two neigh- 
boring antistokes lines and have the same asymptotic series 



(see Add endix 16 . 1 1 for more precise statements and details). 

A sweeping correspondence between general transseries and the 
class of analyzable functions has been introduced in the monumental 
work of Ecalle HBJ-EOJ- 

In the context of © , a generalized Borel summation CB of trans- 
series (jlj is defined in • The rest of this section states some results 
of ^H! needed in the present paper; more details are included in Sec- 
tion ^16.21 of the Appendix. 

The formal solutions Q are determined by the equation ((2j) that 
they satisfy, except for the parameters C. Then a correspondence 
between actual and formal solutions of the equation is an association 
between solutions and constants C. This is done using a generalized 
Borel summation CB. 

The operator CB constructed in can be applied to any trans- 
series solution (@J of (J2J) (valid on its open sector Straus, assumed 
non-empty) on any direction d C Straus and yields an actual solution 
y = CBy of (J2J), analytic in a domain S an (see (|137j0 . Conversely, any 
solution y{x) satisfying Q on a direction d is represented as CBy(x), 
on d, for some unique y(x): 



y(x) ^0 (x G d; \x\ — > oo) 



(8) 



y(x) ~ yo (x 6 d; \x\ -> oo) 



(9) 
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y(x) = C k e- A ' k ^ M - k y k (^) 

k>0 

= C k e~ A - k:c x M - k £fiy k (x) = CBy(x) (10) 

k>0 

for some constants C G C n , where Ma = L^QyJ + 1 (L'J * s the integer 
part), and 

oo 

y k (*) = £^kfe (a' = a-M) (11) 

r=0 

(for technical reasons the Borel summation procedure is applied to 
the series 

y k (x) = x W s k (x) (12) 

rather than to s k (x) cf. ©,©). 

In any direction d, CB is a one-to-one map between the transseries 
solutions on d and actual solutions satisfying Q, see |I3J, Theorem 
3. 

The map y ^ CB(y) depends on the direction d, and (typically) 
is discontinuous at the finitely many Stokes lines, see [T3J, Theorem 
4. _ 

For linear equations only the directions AjlR+, j = l,...,n are 
Stokes lines, but for nonlinear equations there are also other Stokes 
lines, recognized first by Ecalle (the complex conjugate directions to 
Pj ;k M + cf. (|35|): see ^Hl)- CB is only discontinuous because of the jump 
discontinuity of the vector of "constants" C across Stokes directions 
(Stokes' phenomenon); between Stokes lines CB does not vary with 
d. 

The function series in (jlUj) is uniformly convergent and the func- 
tions y k are analytic on domains S an defined in (|137|) (for some 6 > 0, 
R = R(y(x),5) > — see Theorem d of flOl . 

2. 3. Heuristic discussion of transasymptotic matching 

There is a sharp distinction between linear and nonlinear systems 
with respect to the behavior beyond Straus- 

In the linear case there are only finitely many (at most n) nonzero 
y k in ([TO]) , and (flT)|) holds in a full (possibly ramified) neighborhood 
of infinity, except for jumps in the components of C, one at each 
Stokes line (see jlHlj also jlj, ^3], [SO] and the references therein). 
The map CB is continuous at the antistokes lines, and thus the trans- 
series y of y is the same on both sides of an antistokes line. What 
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changes at such a direction is the classical asymptotic expansion of 
y, because classical asymptotics only retains the dominant series in 
y, and exponentials in the transseries y exchange dominance at antis- 
tokes lines. From the point of view of exponential asymptotics, where 
transseries are considered rather than just the dominant series, be- 
havior of solutions of linear equations at antistokes lines is relatively 
simple. 

In the nonlinear case however, generically all components yk are 
nonzero and, beyond Strans (for example, in the notations of § El 
below, for arg(x) > vr/2), (|TU)) will typically blow up because of a 
growing exponential (e~ x in this example). 

The divergence of ljlt)|) turns out to mark an actual change in the 
behavior of y(x), which usually develops singularities in this region. 
The information about the singularities is contained in (|1U|) . 

The key to understanding the behavior of y(x) for x beyond S an 
is to look carefully at the borderline region where converges 
but barely so. Because of nonresonance, for arg(x) = ir/2 we have 
?R.(\jx) > 0,j = 2,...,ni. 5 By (JHZl) all terms in (@J) with k not a 
multiple of ei = (1, 0, ...0) are subdominant (small). Thus, for x near 
iR + we only need to look at 

yW(x) = Y J C k l z- kx x kM ^y kei {x) (13) 

The region of convergence of (|13|) (thus of (|10|l ) is then determined by 
the effective variable £ = Cie~ x x ai (since yfc ei ~ yke 1 -fl/x k ^ ai ~ Ml ^). 
Convergence is marginal along curves such that £ is small enough 
but, as |x| — ► oo, is nevertheless larger than all negative powers of x. 
In this case, any term of the form C k e~ kx x kMl y kei .Q is much larger 
than the terms C[e~ lx x lMl yi ei]r x~ r if k, I > and r > 0. Hence the 
leading behavior of yM is expected to be 

yW(x) ~ ^(de-V^s^o = F (0 (14) 

fc>0 

(cf. (fTTj) ): moreover, taking into account all terms in Sfc ei we get 

OO OO OO p /^\ 

r=0 fc=0 j=0 

Expansion (|15|) has a two-scale structure, with the scales £ and x. 

It may come as a surprise that each Fj is a convergent series in 
^ (though the whole expansion (|15j) is still divergent). 

5 In the notations explained below in 321 =0 for j > n 1 . 
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It turns out that the reshuffling 1)15)1 is meaningful and yields 
the correct asymptotic representation of yW, and therefore of y, be- 
yond the upper edge of S an . In fact, (fT5|) extends © right into the 
regions in C where y is singular, as near as (under mild assump- 
tions) 0( e - co ^-N) of these singularities. Once these two scales are 
known and once the validity of ()15)1 is proved for our class of systems 
(Theorems 1 and 3 below), it is easier to calculate the Fj by direct 
substitution of (|T5|) in ((21) and identification of the powers of x (see 
Remark and ffi).9JI . The exact form of the second scale £ is decisive 
for the domain of validity of the expansion, see H5.2I 

To leading order we have y ~ Fo (see also (fTl)) ) where Fo satisfies 
the autonomous (after a substitution £ = ) equation 

£F(,(£) = iF (£)-g(0,Fo) 

which can be solved in closed form for first order equations (n = 1) 
(the equation for Fq is separable, and for k > 1 the equations are 
linear), as well as in other interesting cases (see e.g. $5.21 ^5.3(1 . 

Assume that Fo(£) has an isolated singularity at £ = £ s . Then 
y(x) must also be singular near x s , if £(x s ) = £ s - Indeed, it is not 
difficult to see (see H4.6j) that there must exist some #(£) analytic at 
£ s so that that / g(t)F (t)dt = 1 on a small circle around £ s . Taking 
x s large we must have by {HJ <f(l + o(l))g(£(t))y(t)dt = 1 + o(l) 
on a small circle around x s . In many instances one can refine these 
arguments to see that the singularities of y(x) and Fo(£(x)) must be 
exactly of the same type. It is clear, on the other hand, that x s form 
a nearly periodic array of points as |a? s | — > oo (see Theorem |2J). 

In the following we will make rigorous these intuitive arguments 
and then proceed to explore further properties and consequences. 

3. Main results 

Let d be a direction in the x-plane (not an antistokes line). Consider 
a solution y(x) of (j2J satisfying Q hence Q. Let (fTfl|) be its rep- 
resentation as summation of a transseries y(x) (see @) on d. Let 
Straus be the sector of validity of y(x) see (jSJ). 

For simplicity we assume, what is generically the case, that no 
pj^k (see 1)35)1) lies on the antistokes lines bounding S trans- 

We assume that not all parameters Cj are zero, say C\ ^ 0. Then 
Strans is bounded by two antistokes lines and its opening is at most 

7T. 

Notations. It can be assumed without loss of generality (possibly 
after a linear transformation in x and renumbering the coordinates 
of y) that 

(a) Ai = 1, and 
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(b) Cj = for j > m (where n\ G {2, ...,n}) 

(c) arg(Ai) < arg(A 2 ) < ... < arg(A ni ) 

(d) Strans is bounded by iM+ and the direction arg(A ni x) = — ir/2 
(which are antistokes lines associated to Ai and A n J, and Strans is 
contained in the right half-plane. 

The solution y(x) is then analytic in a region S an (see ()137)) ). 

The singularities of y(x) that we find are related to the two an- 
tistokes directions bounding Strans- We will formulate the results for 
the direction iM + (and similar results hold for the other direction 
arg(A ni x) = -it/2). 

The locations of singularities of y(x) depend on the constant C\ 
(constant which may change when d crosses Stokes lines). We need 
its value in the sector between iM + and the neighboring Stokes line 
in Strans- Let d! C Strans be a direction in the first quadrant and 
consider the representation (|1U|) of y(x) on d'. 6 From here on we will 
rename d! as d. 

Fix some small, positive 8 and c. Denote 



Z = £(x) = Cie- 



(16) 



and 



£ = < x ; arg(x) G 



Also let 



7T . 7T 

-2 +5 <2 + * 



and 



»(Ajs/|x|) > c for all j with 2 < j < m} (17) 



S Sl ={xe£ ■ <*i} 



(18) 



The sector £ contains Strans > except for a thin sector at the lower 
edge of Strans (excluded by the conditions 5R(Aj-a;/|x|) > c for 2 < 
j < tii, or, if ni = 1, by the condition arg(x) > — f + 5), and may 
extend beyond iM+ since there is no condition on K(Aix) — hence 
3?(Aix) = 3?(x) may change sign in £ and 5^. 

Figure 1 is drawn for n\ = 1; £ contains the gray regions and 
extends beyond the curved boundary. 



3.1. Asymptotic behavior of y(x) in Ss 1 

Theorem 1 (i) There exists 5\ > so that for |£| < 5\ the power 
series 

oo 

F m (0=J2^ykev,m, rn = 0,1,2,... (19) 

fc=0 

6 Ci does not change at the possible secondary Stokes lines dj,k, |k| > 1 lying 
between R+ and iR+. 
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I x I =x 




arg (x) =-jt/2+8 

Fig. 1. Region T> x where l)2ri|l holds, when n\ = 1. The dark gray subregion is Ss 1 . 
Curves like the spiraling gray curves surround points in X (close to singularities 
of y) generate the region T> x . The picture is drawn with ni — 1,A = jQ,a = 
— h,5i = 3 • 10 6 ,a;o = 40. In this case Strans is a sector where | arg(a;)| < 5 — 0. 



converge (for notations see ^jj, |2J) ; \lb\) and for an estimate of 5\ 
see Proposition^. 
Furthermore 



y(x) ~ x~ m F m (t(x)) {xeS Sl , x^ oo) 



(20) 



m=0 



uniformly in Ss x , and the asymptotic representation 120\) is differen- 
tiate. 

The functions F m are uniquely defined by \2U\) . the requirement of 
analyticity at £ = 0, and Fq(0) = ei . 

(ii) The following Gevrey-like estimates hold in Ss 1 for some con- 
stants K\ t 2,B\ > 0: 



\F m (£(x))\ <K x m\B? 



(21) 



m—l 



y(x) - x~ m F m (£(x)) 

j=0 



< K 2 m\B^x~ 



[m G N + , x G S Sl ] 



(22) 
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Comments: 

1. It is interesting to remark that the constant beyond all orders 
C\ is now classically definable in terms of the expansion (|20|) because 
this expansion is unique with its the range of validity, and with the 
given analyticity properties. This is in a sense a generalization of 
Watson's Lemma in the context of transexpansions. 

2. While the classical expansion Q is valid only in any proper 
subsector of S$ 1 n {x : arg(x) < vr/2}, the representation (|20j) holds 
down to a distance going to zero as x becomes large from the (finite- 
plane) singularities of Fo, near which y(x) also develops singularities 
(see Theorem |21 and SJ5J). 

3. A similar picture holds near the lower edge of Strans- The con- 
stant Cj used in which determines the position of singularities 
(see (HHj)) °f y{ x ) related to that direction, is C m . If n\ = 1 then the 
value of C\ is the one in the fourth quadrant (which may differ from 
the one in the first quadrant due to the Stokes phenomenon on 



3.2. Singularity analysis 

We now focus on singularities of y(x) and their connection with sin- 
gularities of Fo- 

3.2.1. Definitions (cf. Figure 1) By © and (O we have F (0) = 0. 
Both Fo and y turn out to be analytic in Ss 1 (Theorems ^i) and 
Eli)); the interesting region is then £ \ Ss-^ (containing the light grey 
region in Figure 1). 

Denote by V a polydisk 

V = {{x-\y) : \x- x \ <pi,|y| < p 2 } (23) 

where g is analytic and continuous up to the boundary. 

Let S be a finite set (possibly empty) of points in the £-plane. This 
set will consist of singular points of Fo thus we assume dist(S', 0) > 8\. 

Denote by 1Z~ the Riemann surface above C \ S. More precisely, 
we assume that 1Z~ is realized as equivalence classes of simple curves 
r : [0, 1] i— > C with -T(O) = modulo homotopies in C\S. 

Let T> C TZe be open, relatively compact, and connected, with the 
following properties: 

(1) Fo(£) is analytic in an ex>-neighborhood of T> with et> > 0, 

(2) sup-p |F (OI : = P3 with p 3 < p 2 

(3) V contains : |f | < 5i}. 7 

7 Conditions (2), (3) can be typically satisfied since Fo(£) = £, + 0(£ 2 ) and 
Si < p2 (see also the examples in JU; borderline cases may be treated after 
choosing a smaller S\. 
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It is assumed that there is an upper bound on the length of the 
curves joining points in V: dx> = sup a mfr c x> ; a,6er length(.T) < 
oo. 

We also need the x-plane counterpart of this domain. 

Let R > (large) and let X = n {x G £ : \x\ > R}. 

Let r be a curve in T>. There is a countable family of curves 7iv 
in the x-plane with £(7jv) = r. The curves are smooth for |x| large 
enough and satisfy 

7jv(t) =2iV7ri + ailn(27riiV)-lnr(t) + lnCi + o(l) (N — > oo) 

(24) 

(For a proof see Appendix Hfi.MI ) 

To preserve smoothness, we will restrict to \x\ > R with R large 
enough, so that along (a smooth representative of) each r € T>, the 
branches of £ _1 are analytic. 

If the curve r is a smooth representative in T> we then have 
^~ 1 (r) = Ujvgn In where 7jv are smooth curves in {x : \x\ > 2R}\X. 

We define V x as the equivalence classes modulo homotopies in 
{x G £ : \x\ > R} \ X (with oo fixed point) of those curves which 
are completely contained in £ n {x : \x\ > 2R}. 

Theorem 2 (i) The functions F m (£); m > 1, are analytic in T> 
( note that by construction Fo is analytic inV) and for some positive 
B,K we have 

|F m (0| <W, £gX> (25) 

(m,) For i? large enough the solution y(x) is analytic in T> x and 
has the asymptotic representation 

oo 

y(x) ~ ^ x~ m F m (£(x)) (x G P,, |x| -> oo) (26) 
m=0 

In /ac£ ; £/ie following Gevrey-like estimates hold 



rn—l 

y(x) - ^ x^F^(x)) 

3=0 



< ^ml-B^x] 



(m G N + , x G Pj;) 



(27) 



(mj Assume Fo /ias an isolated singularity at £ s G S and i/iai i/ie 
projection of T> on C contains a punctured neighborhood of ( or an 
annulus of inner radius r around) £ s . 
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Then, if C\ 7^ 0, y{x) is singular at a distance at most o(l) 
(r + o(l), respectively) of x n £ £ -1 ({£ s }) n "D x , as x n — > 00. 
The collection {x n } n ^ forms a nearly periodic array 

x n = 2nni + a x \n{2nm) + In C\ - In £ s + o(l) (28) 



as n — > 00. 

Some of the conclusions of the theorem hold with T> noncompact, 
under some natural restrictions, see Proposition |HJ 

Comments. 1. The singularities x n satisfy C\e~ Xn x" 1 = £ s (l + 
o(l)) (for n — * 00). Therefore, the singularity array lies slightly to 
the left of the antistokes line iM+ if 3?(ai) < (this case is depicted 
in Figure 1) and slightly to the right of iM+ if 3i(ai) > 0. 

2. In practice it is useful to normalize the system © so that ct\ is 
as small as possible (see the Comment 1 . in § 15.21 and § lfi.7|) . 

3. By (|27|) a truncation of the two-scale series (|26|) at an m de- 
pendent on x (m ~ \x\/B) is seen to produce exponential accuracy 
( e -N/- B l), see e.g. 02|. 

4. Theorem|51can also be used to determine precisely the nature of 
the singularities of y(x). In effect, for any n, the representation (|26|) 
provides o(e~ K ^ Xn ^) estimates on y down to an o(e~ K ^ Xn ^) distance of 
an actual singularity x n . In most instances this is more than sufficient 
to match to a suitable local integral equation, contractive in a tiny 
neighborhood of x n , providing rigorous control of the singularity. See 
also Sand §5\ 



3.3. Singularities for weakly nonlinear systems 

In this section we take g meromorphic in a small enough neighbor- 
hood of (0,0), but nevertheless analytic at (0,0), and only weakly 
nonlinear. Such could be the case if in a sufficiently large neighbor- 
hood of zero only one component of g is singular, and the singular 
manifold of g is approximately a hyperplane. 

Let y(x) be as in Theorem ^ denote f(y) = —A + g(0,y). By 
Theorem^(i) we have yi(x) ~ £(x) and yj(x) = 0(£ 2 ) <C yi(x), j = 
2,...,n when 1 3> |£(x)| 3> |x|~ 2 . For definiteness we assume the 
component 52 to be the only one singular in some neighborhood of 
(0,0). The precise assumptions are 



16 



O. Costin, R. D. Costin 



/ 2 (y) 



h 



Aiyi + ei(y) 

_ -A 2 y2 -72y? + e 2 (y) 



/i(y) 



My) 

- A i l/j - 7j 2/i + e j (y ) (j ^ 1 , 2) 



(29) 



My) 



(y) 



fc=i 



where Qj,7j G C with ||a|| < 5(^2) 1 (see (|23j) ). ej are analytic and 
satisfy 



for some small positive e. 

Choose xo > large enough so that the function g(x _1 ,y) is 
analytic in V if p\ = Xq 1 (see (|23|l) and 



Theorem 3 For almost all values of the parameters 7j,aj with j = 
1, ...,n the following holds. Ife is small, Fo(£) andy{x) are not entire. 
Fq has isolated square root branch points on its circle of analyticity 
and correspondingly y has arrays of branch points for large x. 

More precisely, there exists £ s so that Fq is analytic for |£| < |£ s | 
and y is analytic if \£(x)\ < |£ s | — 5(x) where S(x) — > as x — ► 00. 
Furthermore, Fq is analytic on the Riemann surface of the square 
root at £ s and there is an array x n G 5^ to^ = £s + o(l) as 

|sc n | — > 00 so that y(x) is analytic on the Riemann surface generated 
by curves that encircle at most one of the x n . Near £ s and x n we have 



ej{y)\ < e for |y| < p 2 for j = 1, n + 1 



(30) 



Slip < e 



(31) 



F (0 = F A (K - e-) 




y(z) = y^((^ - x n ) 




(32) 



respectively, where yA and Fa are analytic at zero. 



4. Proofs and further results 



4-1. Further results needed 



A possible proof of Theorem ^ using only classical asymptotics con- 
cepts is sketched in ^6,101 The proof given in this section uses some 
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results in exponential asymptotics. We need the following facts proved 
in H3|. 

Denote by d the complex conjugate of d. We have 

Yu(x) = [_Y k (p)e-P x dp (33) 

Jd 

where the functions Y k have the form 

Y k (p)=p- ka '- 1 A k (p) (34) 

(Lemma 20 in [HI]), with A k analytic near zero, and along curves 
towards infinity avoiding the points pj- k defined as 

Pj;k = Aj — k • A , j = 1, m , ke Z™ 1 (35) 

Because in Strans we have 3?(Ajx) > for j < m and since k > 
it follows that (1) the j>j ;k have no accumulation point and (2) only 
finitely many of them are in Strans = {d : d C Strans} ( see |13| for 
more details). In particular, if pj ;k g" d then there exists a\ > and 
an ai-neighborhood d ai of d (i.e. d ai = {x; dist(x, d) < ai}) where all 
A k , k > are analytic. 

There exist positive constants K\, vq such that for all k > we 
have (113,, Prop. 22(ii) for W k = Y k in V' mv and T k . /3 /_i) 

sup \Y k (p)e- M \ < K lkl (36) 

Also, y k is the classical asymptotic expansion of y k ( JSj, Theorem 
3): 

y k (x) ~ y k (x) (x G de t6 , 6 G (-|, ~) ; x -+ oo) (37) 

In the following we need a better estimate of A k (see ("*"""") and 
1(3*3*1) ). than given in [*f*3*| . 

Proposition 4 For some a±, 62 > and all k and Z we have 

\A.l\p)\ < |r(-k-a / )r 1 Z!ar'e iy ° |p|+ai ^ |k| (38) 
uniformly in d ai . 

The proof amounts to minor modifications in a proof in [*""">]. We 
detail these modifications in Wi.6l 
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4-2. Proof of TheoremUl 

(a) Analyticity at 0. From PBl . ()34[). and Watson's lemma |S] we get 
as x — > 00 along d 



^0 m! ^ 



£ Ai m) (0) 



m>0 



m\x m kc *' 



s m-k.«'-l e - 



5 *=E- £rt)(0)r(m " k ' Q,) 



m>0 



m\x m kc *' 

1 



(39) 



m>0 



so that from Proposition 0] it follows that F m (£) = X^fcLo yfcei;m£ fc 
converges if |£| < 5±, where 5\ < 62, in which ball we also have 



|F m (OI < K x B?m\ 



(40) 



for some B\. 

(b) Asymptoticity. We first note that using again Watson's lemma, 
by ()38j) we have 



M ~ 

Yk-^E^ 
z — ' x l 



1=0 



<e ai Ir^k-a') -1 ! 5 2 ^ ai ' M - l {M + 



< 



v-lc-[k| k-cx'+Af+l„ -Af-1 



r(-k-a')- 62 % «1 



< 8~ M {M + l)\ ai - M - 1 
It is then convenient to write, for x € 5^, 



„k-ct'+M+l 



(41) 



y (x) =Y^C{e- kx x kai y kei (x)+ C k e- k ' Aa; x kQ y k (x) (42) 

k=0 k^0;k^fcei 
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and (cf. (PI) . flUl ) 

oo 

= ^C£e- kx x kai y kei (x) + 0(e~ c ^) := y^(x) + 0{ e - c ^) 



k=0 

Now, 



oo / M \ 

- rn 



k=0 k=0 \m=0 / 

oo / M \ 

+ Y^C k e~ kx x k ^ y kei {x) - £ yk ei ; m x~ m (43) 



k=0 \ m=0 



and theorem ^ follows from (|41|) . Differentiability simply follows from 
the fact that (1201) holds in a nontrivial sector. 



4-3. Special Gevrey estimates 

In the proofs of the main theorems we need optimal estimates of high 
order derivatives of functions of the form ip(z, Ylt=i a kZ k ) in terms 
of estimates of a^, when Ylk^o^k^ are Gevrey type series. Because 
of the truncations involved, the estimates do not follow from Gevrey 
theory. 

Let (f be an analytic function in the polydisk V = {\z\ < pi, \y\ < 
ft} C C x C and continuous up to the boundary. 
Assume first that the series 

oo 

k 



k=i 

converges and denote by a^ the truncation 

m 

a M (z) = ^ afe 



, - 

k=0 



Then 4^Lp(z; a(z))| 2= o is a polynomial in ai,...,a m , of degree 1 in 



dz 

am, • 



1 d m 1 d m 

-(p(z;a(z))\ z = = a y <^(0;a o )a m + — - — -tp{z\ ah~ 1 ^ (z))\ z=0 



ml dz m ml dz m 

(44) 
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Relation (|44jl is meaningful even when a(z) is only a formal sum 
(with no convergence conditions) — in the sense that the LHS is the 
coefficient of z m in the formal series expansion of (p(z;a(z)) at z = 
0. We are primarily interested in the Gevrey— 1 character of a(z) 
(meaning that for some c\,C2 > we have |a&| < cic^kl; see also |42p. 
Proposition [5] below is formulated in a way that permits an inductive 
proof of Gevrey type inequalities, when are defined recursively. 

Proposition 5 Assume p2 — |ao| > 0. There exists a positive C so 
that: for any B,K > and any {a.k}k=i...m such that \&k\ < KB k k\, 
k = 1, 2, m we have 



1 

ml 



< if 2 -B m_1 (m- 1)!(1 + Cm- 1 log 2 m) 

2=0 

(45) 



/or some K% (see \55)) for an estimate of K2). 
For the proof we need the following result. 

Lemma 6 Let a,b satisfy 1 < a < b. There exists C = C(a;b) such 
that if Bmz := Z £ (a,b) then 

^ fc=1 1 <Cm^(lnm) 2 (m € N) 



m\m- m Z m [l - Z- 1 ]- 1 + m^Z 



(46) 



Proof of Lemma In this proof (and in the proof of Proposition [SJ) 
we write 0{f{m)) for terms that go to zero not slower than than 
f(m) uniformly in B, Z (and K). 

Let k m = \ra/Z\. For k < k m the terms k\m~ k Z k are decreasing 
in k, and increasing for k > k m . Thus 

km y 27 2 ?s7 3 

Y B k z k k\ < - + —r + m—^- = Bz(l + Oim.- 1 )) (m -» 00) 
k=\ 

(47) 

Denote p m = [2 In m/ In ZJ . For m large enough we have 



m > k m +p m , 1 -p m /rn > 1/2, p m > k v 



and for p < p Tl 



i=o 



'' ^ > e 2 2,n > << 
m 



p(p+l) _ Pm.(pm + 1) 



(48) 
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N2 k\(Bz) k 
Denote a^ 2 = — — - . — ; we have 
Nl ^ m\(Bz) m 
k=Ni v ' 



k=0 j=0 

while using (|4.3|) it follows that 



m—pm — 



Pm(Pm + l) 



k=0 



1 + 0(m 1 (lnm)' 
(1 - (mBz)- 1 ) 



(50) 



For k £ (k m ,m — p m ), because the terms in the sum are increasing 
we get 



a 



m-p m 



< m- 



(Pm+l)/m 
2?2 lnm/ InZ 



m~ L + 0(m~ 2 (lnm) 2 ) (51) 



Combining (jlHJ) (fBTljl. (|5T|) . and (|I7)l. Lemma El follows. 
□ 

Proof of Proposition O We keep the requirement of uniformity 
with respect to B, K in the notation O(-), as in the Proof of LemmaEl 
Let p = min{/9i,p2 — l a o|} (cf- the beginning of H4.3|) . For small s 
and y we have 



\ip(s, y) - <p(0, 0) - d s <p(0, 0)8 - dy<p{0, 0) • y| 

< 2(n+ p 1 / HI (M 2 + llyf) = ^(N 2 + l|yf) 

We choose a circle of radius r m , where 



(52) 



(m-l)!(flr m ) 



171 — 1 



Br r . 



(53) 



For large m we have Bmr m = e + 0(m _1 ), and we also see that 
the assumptions of Lemma El are satisfied. In particular we have for 
Isl = r m that 



m— 1 
fc=l 



< KBr m {\ + (1 - l/e) -1 )(l + C^m^lnm) 2 )) (54) 



Noting that £ s m+1 ^ m ~ 1 \s)ds = we have, using (JHSJ) and <|52|). 
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1 d" 



ml dz 7 





— <£ 


z=0 


ZTTi J r 



ds 



^(s;aN( s )) 



Q m+1 



< Ml + ^ 2 ^ 2 )|r m | 2 - m (l + O^^lnm) 2 )) 

+ K 2 B 2 )B- l {B m ~ l {m - 1)!)(1 + C^m^lnm) 2 )) (55) 



□ 



Remark 7 ^4 direct calculation shows that the expansion in \2b}) is 
a formal solution of for large x iff the functions F m are solutions 
of the system of equations 



_d 
d_ 



Fo 

F m + :vf 



r 1 (iF -g(0,F )) 



(56) 



ai-rrFm-i + Rm-i for m > 1 (57) 
d? 



where N is the matrix 



rH9yS{0,F )-A) (58) 
and the function R m _i(£) depends only on the F& with k < m: 



171— 1 



(m-l)/ + A F m _! 



1 d m 

m! dz m 



m— 1 



2=0 

(59) 



fsee a/so j^^j. 



4-4- Proof of Theorem^ (i) 

By Theorem^ F m , m > are analytic for |£| < Si. Furthermore, 
F m are analytic in the e^-neighborhood of T> since by assumption, 
Fo is analytic there and equations (|57|) are linear for rn > 1. 

We set as initial conditions for (fST|) in 2? the values of F m (£o) 
provided by Theorem^at a point (o£P with |£o| £ (^i)^i)- 

For £ E T> and m > 1, (fST|) can be integrated, yielding the recursive 
system 
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F m (0 = M(£;£ )F m (£ )+a 1 [F m ^(0 - M(£;£ )F m -i(fo)J 
-a x t M(Z;s)N(s)F m _x(s)ds+ f M^;s)It rn - 1 (s)ds (60) 

ho ho 

where M(£; C) is the fundamental solution of 



^- + NM = with M(£;Q\e=c = I (61) 
Direct estimates in (|6ll)l and (|ST?|) using (flU)) give 



||F m ||2? < MK x m\Bf + (|ai| + M)\\ 

+ 2|a 1 |d I ,M5 2 - 1 (||g|| + ||yi||)|| F m -i \\t> 

+ 25~ 1 M 



(m + A) ||-P m -i|| +d v 



1 H m 



to! dz r ' 



(62) 



where 



M= sup ||M(e;C)ll; ||F||x» = sup|F(0|; l|g||= sup |g(z,y)| 
Uev £sx> | 2 |<pi,|y|<p 2 

A = max ||i^(0|| < ih/2)- 1 {\\g\\ + P||) (63) 



Choosing K, B large enough, the proof of (|25|) is immediate induction 
from (|62|) and Proposition [5J 



4-5. Proof of theorem^ (ii) 

We will prove (|26l) at each point x = x a E (with uniform estimates 
on Pa;). x a is the endpoint of a curve 7at in I? x with £;(7at) = r curve 
in V and satisfying (|24|). 

Denote a = £(x a ). If |a| < <5i then ()26jl follows from Theorem ^ 
so we assume \a\ > 5\. Then we can choose r to go from along a 
direction up to the circle |£| = 5\, not re-entering the circle. 

Let to such that £ = r(t ) G (|5i,5i) and denote -T = r|[ to l ], 
7° = 7Jv|[to,i]i then 7 lies in a bounded region. 

We prove Q26JI in a small, connected, simply connected neighbor- 
hood A/lo of 7 . 
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Denote S(x) = y(x) - F^(x). 

To estimate S(x) we use a contraction argument if m is not too 
large {Case I) and a direct argument for m large (Case IT). 

Let Co > be small and m a be the maximal integer such that 
m\(2B/\a\) m < c . 

Case I: m < m a 

First the differential equation satisfied by F^(x) will be written, 
then the equation of S(x), and finally the independent variable will 
be changed to £, yielding (|65|). 

With the notations introduced in H4.3I and denoting by F(x) the 
formal series ^2^ = QX~ k F m (^(x)) we get from Remark[7] 



ax V x 



+ Z^T AF m(^)) + ^-JrCeW) - -^F m fc(s)) (64) 



Using and @, a direct calculation yields the equation for 
d(x). The map £(x) is a biholomorphism of A/lo onto a neighborhood 
Afpo of .T ; changing the independent variable from x to £ we get 



^ +,M (4> +Ti (^) +T2 (4) (65) 

where 



1 



1 



x£ a/x — 1 
1 1 
£ a/x — 1 
1 1 
£ a/x — 1 



-oul + A + c4^(0,F ) 
g(~,FW+8)-g(-,F* 

X X 

g(-,F [m] ) -g(-,F [m] ) [m 



-d y g(0,F )d 



1 



where To, 1.2 are clearly well defined for small enough £ and S. Fur- 
thermore, they are well defined for £ G A/f° if |<5(£)| < (/?2 — /°3)/2 
and for R large enough (see Appendix t jfi.5|) . 

As in (JfiOj) we obtain for 6 the integral equation 
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S = J (6) where J = Jo + Ji 
with Jo(0 = M(£; £ )<5(£o)+ f Mfc s)T 2 (-i^ ds 



+ J M(^,s)T (^-^jS(s)ds (66) 

Let B be the Banach space of C™-valued analytic functions 5 on 
A/po, continuous up to the boundary, and satisfying <5(£o) = 0, with 
the norm \\S\\ = sup^^ |<5(£)|. 

The integral operator J of ()66() is defined on the ball of radius 
(P2 — P3)/2 in B. We will show that it invariates a ball in B and that 
it is a contraction there. As a consequence, the integral equation (|66|) 
has a solution (5 which is analytic on A/po, therefore y(x) is analytic 
on A/y) ; we will also obtain estimates for which will prove (|26[) in 
Case /. 

We will denote by const a constant independent of a, N, m, B, en, R. 
It will be assumed that B,R > 1, cq < 1. 

Note first that the assumption of Case I implies 

m/|o| < const (67) 

To estimate J{8) note first that 

llm „ const . . 

T < 68 
\a\ 

By (EU when a is large we have |5(£ )| < K 2 {m + lJLBj" + 1 a -m-1 . 
By Theorem [2] (i) , since £ varies in a compact set independent of a, 
and then \M\ < M as in (|63j) . Also estimating derivatives with the 
Cauchy formulas on circles \x— x'\ < pi/2; |y — y'| < ejy/2 and taking 
a large so that |o| < 2|a — du| we get 



|Ti|| < " y — j-r + < const — + ||<5|| ||<5|| (69) 



and 



|9«Ti||< rr ^— + ^P- (70) 
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Also 



M(£;s)T 2 

& 

and using (|BT)l 

so that, 
IIW)II < 



< 



2MKd v (2B) m m\ (2\a 1 \d T > m 



\<i\c p lal 



< consi 



(2B) m m\ 



(71) 
(72) 



M<5(£ )+ / M(£;s)T 2 

& 



ris 



< const 



(2 J B) m m! 



Prom (jggj). PH> we get 



||Ji(5)|| < const ||5|| ^ + ||<5 



1 



< const cq (73) 



(74) 



Also, from (|HH|) . (fHT^ . (fTHj) 



\\J{8 l )-J{8 2 )\\ < COnst -r-r + \\8\ 

a 



(75) 



It is easy to see that for positive constants R, Kq large enough, and 
Co small enough the following holds: if \\S\\ < KqCq then : \\S\\ < P2 ~ P3 ; 
11^(5)11 <# co, also 

(76) 



and ||jr(«5i) - J(S 2 )\\ < X\\Si - S)\\ with A < 1. This shows the 
existence and analyticity of d. 

Finally, to obtain the needed estimate (|26j) note that using 1)73(1 . 
(|75|). we get 

||<5|| = ||J(5)|| < HJbll + \\Ji(S)\\ < const^^ml + h\6\\ 

\a\ m 4 

so that using (|67|) and Lemma QUI 

(2B) m , (2,B) m + 1 / 
o < const — — m! < const — — ; — — — (m + 1)! 

ii H | Q |m | a | m+1 

< const — ; — (m + 1)! 



m+l 



which concludes the proof in Case I. 
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Case II: m > m a 
In this case 



Since 



1511 < 



yO) 



m\(2B/\a\) m > c 
1 



(77) 



Fo Fi 

x 



k=2 



using the result of Case I to estimate the first term 



< const 



(5B) 2 2\ 



m max 



K2\B 2 Km\B r 



and since (2B/\a\) 2 2\ < c < (2B/\a\ 



x 



i m m! 



8B 2 2\ T ^ , (2B) m . ,, (2B^ r " +l 

const— — — +mivm!— < const (m+lj! ( - — - 

|a| m \ |a| 



25 
a 



Fro m (|77|l |a|/(2B) < consi m Cq < const Cq" 1 ; using this, (|77|) 
and (1140(1 we finally get 



<K(c )(m + l)! (4B/co) 



m+l l^l-™- 1 



where -ftT(co) is a constant dependent of Co- □ 

We should stress that while the estimates in this proof clearly 
show the Gevrey character of the expansion, they are very far from 
optimal. In fact the substantial increase in B in the arguments was 
artificially introduced to make the calculations less cumbersome. 

The following is an extension, in some respects, of Theorem |U(ii). 



Proposition 8 Assume T> is not necessarily compact, r is a curve 
of possibly infinite length in T> with the following properties: 

(a) For some e > 0, Ti^^, S) and N(z) are analytic for z in an 
e neighborhood of T and for \S\ < e and in addition Ti^z, 5) = 
0{zd,S 2 ) 

(b) M(£>£i,o) is bounded in an e neighborhood of T and for some K 

and all £ € r we have Jj£ q M(£,£i,o) d|s| < K (where \M\ is some 

Euclidian norm of the matrix M(£,£io),)- 

Then the conclusions of Theorem^ (ii) hold in the x domain V x 
corresponding to T> . 



Noting that 



M(£,a,o; 



d|s| is a finite measure along r, the proof is 



virtually identical to the proof of Theorem |21 
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4-6. Proof of Theorem^ (Hi) 

We need the following result which is in some sense a converse of 
Morera's theorem. 

Lemma 9 Let B r = {£ : |£| < r} and assume that /(£) is analytic 
on the universal covering of B r \{0}. Assume further that for any 
circle around zero C C 5 r \{0} and any g(£) analytic in B r we have 
§c fiOdiO^ = 0- Then f is in fact analytic in B r . 

Proof. Let a G -B r \{0}. It follows that fr f(s)ds is single-valued in 
i? r \{0}. Thus / is single- valued and, by Morera's theorem, analytic 
in £> r \{0}. Since by assumption <f c f(£,)£, n dt; = for all n > 0, there 
are no negative powers of £ in the Laurent series of /(£) about zero: 
/ extends as an analytic function at zero. 

To show Theorem |21 (hi), assume £ s is an isolated singularity of 
Fo (thus £ s 7^ 0) and X = {x : £(2) = £ s }. By lemma 03 there is 
a circle C around £ s and a function analytic in B r {£ — £ s ) such 
that j> c Fo(£)g(£)(i£ = 1- In a neighborhood of x n G I the function 
/(x) = e~ x x ai is a biholomorphism and for large x n 



It follows from lemma^lthat for large enough x n y(x) is not analytic 
inside C either. Since the radius of C can be taken o(l) Theorem |2] 
(hi) follows. 

Note. In many cases the singularity of y is of the same type as 
the singularity of Fq. See ^Slfor further comments. 



4-7. Proof of Theorem^ 

As in £14.51 we can reduce to the study of © in Mpo , where the func- 
tion £(x) is biholomorphic and we can change variables to £. In this 
variable both and (|56[) assume the form (where x = and F 
is F or y) 




(78) 
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at, 

dF 2 _ \ 2 F 2 - l2 F? + e [ 2 l] (x-\F) 



^ = XjFj " ljF ' + e ^ {x ' 1,F) {j + h 2) 



(79) 



where 



hte,F) = l-^a k F k + ( & 1 (x- 1 ,F), 
i=i 

with analytic in V and ||e^||-p < e, % = j, ...,n + 1 (for Fo, we 

havee' 1] (^ 1 ,F ) = e J (F )). 

Generically a 2 ^ 0, then the analytic change of variables to 
Fi, h, F3, F n leads to a system of the form 



+ 



-A 2 + A 2 ajFj - a 2 j 2 F } 



XA 



^ w -,./•?•<; (i> 2) 



(80) 



(81) 



(82) 



The substitution F x = £ + fi,Fj = bj£ 2 + ^ (J > 2), with 
(Xj - 2)~ 1 7 i , in ftt ty -jB5 )l yields 



d£ 



d/i 



= ^ + ai(A 2 -/*)(£ + A) + £ - ^jh)(b m f + fj) 



i>3 



+(/i + 2 E fl A-^-f 

i=3 



[3] 
2 



= A,/, + 2f7i/i + 7,/i 2 - f (j > 2) 



(83) 
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According to the hypothesis of Theorem 01 it is useful to analyze first 
the equation (describing the leading order behavior of h) 



hti = {hC 1 + di + d 2 0h + {-\2C 1 + d 3 + d A Z); h(0) = 1 (84) 

(this Abel type equation cannot be solved in closed form, in general). 
In integral form, 



H0 2 = i + dfk 2 + 4k 

+ / (df ] s + df)h(s)ds + 2\ 2 (h{s)-l)s~ 1 ds (85) 
Jo Jo 



Lemma 10 (i) Equation {8$ has a unique solution ho analytic at 
£ = 0, with ho(0) = 1. 

(ii) For a generic set of d\, ...,d^ the solution ho is not entire and, 
on the boundary of the disk of analyticity, hg has square root branch 
points. 

Proof, (i) It is straightforward to check that, since A2 N (see 
then H85|) has a (unique) formal solution of the form h = 1 + 
EfcLi ~ h kt, k (where h x = (df 1 + df)(2- dfy 1 ). To show h converges 
we take h = 1 + Y,k=l ~ h kt + £ M ^m(0 in ©: 



2i M h M {i) 

= Qm 2M h M (o 2 +e i R(o-4 ] / (t M + x>* M+fe )M*)^ 

J ° k=0 

(86) 

with analytic, or 

2h M (o = Q(m M h M (o 2 +R(o- 

df / (s M ~ 1 + Y j b k s M+k e +1 )h M {is)ds = 2J{h M ) 
Jo k=0 

J is manifestly contractive in the sup norm for small £, if M > 
(ii) The proof, elementary but delicate, is given in ^16.81 □ 



4 
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Lemma 11 Let £o be a branch point singularity on the boundary of 
the disk of analyticity of ho (see Lemma fflJl (ii))- Assume in \8ty 
is small enough and analytic in a (large enough) neighborhood m£, F 

o/(£o,F (&))- Then 

(i) For some < 5\ < 62, F(£) and h(£) are analytic in the cut 
annulus {£ : |£ - £ | G (Si,S 2 ), arg(£ - £ ) / 0}. 

(ii) h and F have a square root branch point at some £ s with 
C s -Co = 0(\\e^\\). 

Proof. We substitute h = ho + fi cf. Lemma ITU1 f satisfies a system 
of the form 

£f = JV(£)f + e(£,f 2 ) 

or 

f = fi + M(0 J M' 1 (s)e{s,f 2 (s))ds 

where the matrices N \M , M _1 and e are analytic in Mp°- Part (i) 
follows now in the same way as Theorem |2] (ii), and \\y — Fo\\j\f r0 = 

0(||e^||). (ii) In a small neighborhood of £0 by part (i) and Proposi- 
tion (generically) j| h 7^ and we may change variables in (|83|) so 
that h is the independent variable (and £ = £(/i). We note that 

d£ £h 

dh~ A(f,Oh + B(f,0 + e(f,h,0 

and /i, e are small while generically S(/(£o)>£o) is n °t small. Then 
(|83|) in the variable h, with initial condition £(/tQ + 0(||e[ 3 l ||)) = £0 + 
0(||e[ 3 ]||), has a solution which is analytic near h = 0. Furthermore 
it is easy to see that (under the same genericity assumptions) we 
have d h Ch=o = but det5 feft (/i,£, .-,f n )h=o + and then Fj(g) = 
Fj((€ ~ £o) 1/2 ) with Fj locally analytic. 

5. Examples 

5.1. Example 1 

We first illustrate how singularities of solutions are found (using 
transasymptotic matching) on a first order Abel equation 8 : 



u = u 3 - z (87) 
The authors are grateful to A. Fokas for pointing out to this example. 
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the first example on which nonintegrability was shown using Kruskal's 
poly-Painleve analysis |34j . 

The study of (|87|) is done in the following steps. Classical asymp- 
totics of differential equations |5J shows (and it also follows from the 
analysis below) that for z — > oo with argz G (^j7r, ^jvr) there is a 
one parameter family of solutions u = u(z; C) such that u(z; C) = 
z 1 / 3 (l + o(l)). Then u ~ u = z 1 / 3 Y^=o z 5fc% ■ The parameter C may 
be chosen to be the constant beyond all orders, see ^5.1.21 

After proper normalization of (|87|) (see H5.1.1|) Theorems ^ an d El 
are applicable and provide a global asymptotic description of u(z; C) 
in a region where the solution is analytic and surrounds its singulari- 
ties for large z (Proposition I13|) . These are algebraic branch points of 
order —1/2 (see 1)102(1 ) and their location, dependent on C, is deter- 
mined asymptotically. Conversely, the C of a particular solution can 
be determined from the asymptotic location of one singularity. 

5.1.1. Normalization Formal solutions provide a good guide in find- 
ing the normalization transformations. A transformation bringing the 
equation to its normal form also brings its transseries solutions to the 
form (|10|) . It is simpler to look for substitutions with this latter prop- 
erty, and then the first step is to find the transseries solutions of 

(EH). 

Power series solutions. Since at this stage we are merely looking 
for useful transformation hints, rigor is naturally not required. Sub- 
stituting of u ~ Az p in (|87|) and looking for maximal balance [S] give 
p = 1/3, A 3 = 1. Then u ~ Az 1 / 3 + Bz q with q < 1/3 determines 
B = I A 2 , q = -4/3. Inductively, one obtains a power series formal 

solution u = Az l / 3 {1 + £Xl uo,kZ~ 5k/3 ). 

General transseries solutions of \8'l\j . In order to determine the 
form of the exponentials in the transseries of u, the method is to look 
for transcendentally small corrections beyond no, by linear perturba- 
tion theory. Substituting u = uq + 5 in (|H7j) yields to leading order in 
S, the equation 



5' = (3A 2 z 2 / 3 + ^\ 5 (88) 

whence 5 oc z 2 ! 3 exp (|^4 2 z 5//3 ) . In Q the exponentials have linear ex- 
ponent, with negative real part. The independent variable should thus 
be x = -(9/5)A 2 z 5 / 3 and 9fc(x) > 0. Then u = x 1 / 5 XXo u 0;kX~ k , 
which compared to © suggests the change of dependent variable 
u(z) = Kx l ^h{x). Choosing for convenience K = A 3 ^{— 135) 1 / 5 
yields 
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ti + — /i + 3/i 3 - - = (89) 
ox 9 

The next step is to achieve leading behavior 0(x~ 2 ). This is easily 
done by subtracting out the leading behavior of h (which can be 
found by maximal balance, as above). With h = y + 1/3 — x _1 /15 we 
get the normal form 



y' = -y + + g{x-\y) (90) 
ox 



where 



rfx-',„) = -Stf +»■) + g - -L, _ _|_ + _L_ (91) 

We see that 

A = l, a = 1/5, and thus £ = Cx 1/B e~ x (92) 

5.1.2. Definition of C for a given solution u(z) After normalization 
l|9~0~|) the results in [2] apply, and the constant C is uniquely associ- 
ated to a u(z) on a direction arg(z) = </> as the limit 

c = & ec*)" 1 («(*)- £ ^Sw) ( 93 ) 

arg(z)=</> y fc<|a:(z)| / 

This limit exists for all € (^j 71 ") ^j 71 ") an d is piecewise constant, 
with one jump discontinuity at the midpoint of this interval. The 
value of C relevant to the singularities of u is the one nearest to the 
edge of Straus where these singularities are calculated, as follows from 
Theorem ^ 

5.1.3. Finding the two-scale expansion \2b]) Having the second scale 
given by ()92|) and all the conditions of Theorem ^ satisfied, the sim- 
plest way to calculate the functions in y = Y^k=o x ~ k Fk{Q is by 
substituting y = y in (JUUJ) and solving the differential equations, as 
in the proof of Theorem |21 (i) ; the equation for Fn(£) is, cf. (Jn7)|) . 

= F (l + 3F + 3F 2 ); i^(0) = 1 (94) 

and, cf. (fo7|). 
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iF' k = (3F + l) 2 F k + R k (F , F k _ x ) 

( for k > 1 and where i?i = ^F 3 ) (95) 

The first term F$ of the expansion of u is then given by 

t = ZoF (0(Fo(0+u>or e (F (0+uJt)-~ e (96) 

with Co = 3- 1 /2 eX p(-i^^/3), too = \ + and = 5 + The 
functions Fk, k > 1 can also be obtained in closed form, order by 
order. 

By Theorem^ the relation y ~ y holds in the sector 

S Sl = {xeC: arg(x) > + 5, \Cx 1/5 e~ x \ < 

for some <5i > and any small 5 > 0. 

Theorem 3 insures that y ~ y holds in fact on a larger region, 
surrounding singularities of Fq (and thus of y) . To apply this result we 
need the surface of analyticity of Fq and an estimate for the location 
of its singularities. 

Lemma 12 (i) The function Fq is analytic on the universal covering 
TZs of C \ S where 

Z = {t P = ("1)^0 exp(p 2 vrV3) : p 1)2 G Z} (97) 

and its singularities are algebraic order —1/2, located at points lying 
above S. 

(ii) (The first Riemann sheet) The function Fq is analytic in C \ 
((-oo,£o]U[£i,oo)). 

(Hi) The Riemann surface associated to Fq is represented in Fig. 

2. 

Proof 

Singularities of Fq. The RHS of Q94|) is analytic except at Fq = 
oo, thus Fq is analytic except at points where Fq — > oo. From (|96|) 
it follows that lirnF _ i . oo £ G S and (i) follows straightforwardly; in 
particular, as £ — > £ p G S we have (£ — £ P ) 1 ^ 2 i ? o(£) — ► \/ — £p/6. 

(ii) We now examine on which sheets in 72.= these singularities 
are located, and start with a study of the first Riemann sheet (where 
^b(C) = £ + 0(£ 2 ) for small £). Finding which of the points £ p are 
singularities of Fq on the first sheet can be rephrased in the following 
way. On which constant phase (equivalently, steepest ascent /descent) 
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paths of £(Fo), which extend to \Fq\ = oo in the plane Fq, is £(Fq) 
uniformly bounded? 

Constant phase paths are governed by the equation ^(dln^) = 0. 
Thus, denoting F = X + iY, since f'/f = (F Q + 3F 2 + 3F 3 ) -1 one 
is led to the real differential equation 9f(£'/£)dX + 3t(£'/£)dY = 0, or 

Y{1 + 6X + 9X 2 - 3Y 2 )dX 

- {X + 3X 2 - 3Y 2 + 3X 3 - 9XY 2 )dY = (98) 

We are interested in the field lines of Q98|) which extend to infinity. 
Noting that the singularities of the field are (0, 0) (unstable node, in 
a natural parameterization) and P± = (— 1/2, ±\/3/6) (stable foci, 
corresponding to — ZJq and — ujq), the phase portrait is easy to draw 
(see Fig. 15.1.3*)) and there are only two curves starting at (0, 0) so that 
\Fq\ — ► oo, £ bounded, namely ±M + , along which £ — ► £o anci £ — > Ci; 
respectively. 

(hi) Thus Fig. 3 encodes the structure of singularities of Fq on 
1Z- in the following way. A given class 7 G 1Z~ can be represented 
by a curve composed of rays and arcs of circle. In Fig. 3, in the Fq- 
plane, this corresponds to a curve 7' composed of constant phase 
(dark gray) lines or constant modulus (light gray ) lines. Curves in 
IZs terminating at singularities of Fq correspond in Fig 3. to curves 
so that \Fq\ — > 00 (the four dark gray separatrices Si, S4). Thus to 
calculate where, on a particular Riemann sheet of IZs, is Fq singular, 
one needs to find the limit of £ in (|96J) . as -Fo — » 00 along along 
7' followed by Si. This is straightforward, since the branch of the 
complex powers 6, 6, is calculated easily from the index of 7' with 
respect to P±. □ 

Theorem [21 can now be applied on relatively compact sub domains 
of IZs an d used to determine a uniform asymptotic representation 
y ~ y in domains surrounding singularities of y(x), and to obtain 
their asymptotic location. Going back to the original variables, similar 
information on u{z) follows. For example, using Theorem [21 for the 
first Riemann sheet (cf. Lemma IT""" (ii)) 

V = {\i\ <K I ^(-00,6) U (£ , +00) , |£-£o| >e,|£-6l >e,} 

(for any small e > and large positive K) the corresponding domain 
in the z-plane is shown in Fig. 4. 

In general, we fix e > small, and some K > and define 
Ak = {z : argz G (^vr - 0, ^vr + 0) , |£(z)| < K} and Kk,~ be 
the universal covering of S Pi Ak and TZ z -,K,e the corresponding Rie- 
mann surface in the z plane, with e- neighborhoods of the points 
projecting on z(x(E)) deleted. 
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Proposition 13 (i) The solutions u = u(z;C) described in the be- 
ginning of 33 have the asymptotic expansion 



u{z) ~ z 1 ^ 



+ + f 

V y k=0 



z 5k/3 I 

(as z — ► oo; # £ TZ z -K,e) (99) 



where 



Mz) = x(z) 1 / 5 e- a '( 2 \ and x(z) = --z 5 / 3 (100) 

5 

(mJ In the "steep ascent" strips arg(£) G (01,02)1 |«2 — oi| < 
7r starting in Ak and crossing the boundary of Ak, the function 
u has at most one singularity, when £(z) = £0 or £i> an d u(z) = 
z i/3 e ±27rr/3^ _|_ (i)) as z _ > qq (t ne s ig n i s determined by arg (£),). 

(Hi) The singularities of u(z;C), for C / 0, are located within 
0(e) of the punctures oflZ z -K,o- 

Applying Theorem |21 to (|90j) it follows that for n — > 00, a given 
solution y is singular at points x P;n such that £(x P)n )/£ p = 1 + o(l) 
large). 

Now, y can only be singular if \y\ — ► 00 (otherwise the r.h.s. of (|90j) 
is analytic). If x p>n is a point where y is unbounded, with 5 = x — x Pin 
and v = l/y we have 



— = vF s (w,o) (101) 
di> 

where F s is analytic near (0,0). It is easy to see that this differential 
equation has a unique solution with 5(0) = and that S'(0) = as 
well. 

The result is then that the singularities of u are also algebraic of 
order —1/2. 

Proposition 14 If zq is a singularity ofu(z;C) then in a neighbor- 
hood of zq we have 



u = ±y/=T/2(z - zq)- x ' 2 M{* ~ z o) 1/2 ) 
where Aq is analytic at zero and Aq(0) = 1. 



(102) 
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Notes. 1. The local behavior near a singularity could have been 
guessed by local Painleve analysis and the method of dominant bal- 
ance, with the standard ansatz near a singularity, u ~ Const. (z—zo) p . 
Our results however are global: Proposition 1131 gives the behavior of 
a fixed solution at infinitely many singularities, and gives the posi- 
tion of these singularities as soon as C\ (or the position of only one 
of these singularities) is known (and in addition show that the power 
behavior ansatz is correct in this case). 

2. Eq. (|90|) can be brought to a form similar to that in TheoremOl 
by the substitution y = v/{\ + v) in ()90|) . The result has the form 



v' = -v - 27 — 10 v 2 + — v + g [1] (t-\v) (103) 

where is a now an 0(t~ 2 ,v~ 2 ) polynomial of total degree 5. The 
singularities of v are at the points where v(t) = —1, and are square 
root branch points, as in Theorem|31 whose technique of proof would 
have also applied, if the more explicit formula (|96|) was unavailable. 



5.2. Example 2: The Painleve equation P\. 

The Painleve functions were studied asymptotically in terms of dou- 
bly periodic functions by Boutroux (see, for example, [21] )• Solutions 



:i<S 
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Fig. 4. Singularities on the boundary of Strans for 1871 . The gray region lies in 
the projection on C of the Riemann surface where 199H holds. The short dotted 
line is a generic cut delimiting a first Riemann sheet. 

of the Pi equation turn out to have arrays of poles and they can 
be asymptotically represented by elliptic functions whose parameters 
change with the direction in the complex plane. Joshi and Kruskal 
carried out this type of expansions for generic solutions, which have 
poles throughout a neighborhood of infinity, to sufficiently many or- 
ders to determine how the parameters of the elliptic functions vary 
|29| . |3Uj . and applied this method to solve the connection problem. 
However, there exist special, one-parameter, families of solutions of Pi 
(the truncated solutions, important in applications) that are free of 
poles in some sectors. These solutions have the same classical asymp- 
totic expansion in the pole free sector to all orders and cannot be 
distinguished by classical asymptotics there. They differ by a con- 
stant C beyond all orders which can be determined by exponential 
asymptotic methods. The results of the present paper apply to this 
special family of solutions and give an asymptotic representation uni- 
formly valid at the ultimate array of poles (neighboring the pole free 
sector) and make the link between the position of these poles position 
and the value of C. 

We note that the behavior of the triply truncated 9 solutions, which 
in a sector have C = and, consequently are pole-free in larger sec- 
tors, does not follow immediately from our analysis. But this case can 
be treated by a similar methodology since after continuation across 

9 They are also known as " doublement tronquees" . 
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a Stokes line the value of C becomes equal to a Stokes multiplier, 
generically nonzero. 

This example extends the asymptotic expansions of |16j to larger 
regions of the complex plane, and also to all orders. 

We consider solutions of the Painleve Pi equation (in the form of 
|25| . which by rescaling gives the form in 24) 



^ = % 2 + * (!0 4 ) 

in a region centered on a Stokes line, say d = {z : argz = ir}. 

To bring 1)104(1 to a normal form the transformations are suggested 
by the general methodology explained in ^5.1.11 There is a one pa- 



rameter family of solutions for each of the behaviors y ~ ±< for 



large z along d. We will study the family with y ~ + w , since the 

other can be treated similarly. Its transseries can be obtained as in 
the previous example, namely determining first the asymptotic series 
yo, then by linear perturbation theory around it one finds the form 
of the small exponential, and notices the exponential is determined 
up to one multiplicative parameter. We get the transseries solution 

— oo 

/Ee*y* (105) 

k=0 

where 

£ = £( z ) = Cx- 1/2 e~ x ; with x = x(z) = ^^/^ (106) 
and yt are power series, in particular 




- , 1 '_]__ 1/0;^ 

V0 8V6(-Z)V2 28-3 z5 - (-*)5*/2 - 

We note that in the sector | arg(z) — ir\ < the constant C of 
a particular solution y (see (|109[) ) changes only once, on the Stokes 
line arg(z) = 7r |13j . 

As in Example 1, the form of the transseries solution (|105[l . (|106|> 
suggests the transformation 

(~24z) 5 / 4 /=J 

x = : viz) = \ Y(x) 

30 y ' V 6 V ' 

which, in fact, coincides with Boutroux's (cf. 24 j); Pi becomes 

Y"(x) - \ y 2 (x) + ~ = -- Y'(z) + ^ \ Y(x) (107) 
2 2 x 25 x z 
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For the present techniques to apply equation (|1()7|) needs to be 
fully normalized and to this end we to subtract the 0(1) and 0(x _1 ) 
terms of the asymptotic behavior of Y(x) for large x. It is conve- 
nient to subtract also the (9(x~ 2 ) term (since the resulting equation 
becomes simpler). Then the substitution 



transforms Pi to 



Y{X) = 1 " 25^ + HX) 



1., . 1.2 392 



y + _*_ fc _ 5fc >__ = o (io 8 ) 

Written as a system, with y = (h, hi) this equation satisfies the as- 
sumptions in ^21 with Ai 5 2 = dbl, a% t 2 = —1/2, and then = 
Ce~ x x~ 1 / 2 . The results of the present paper apply to the normal 
form (|108[) of Pi and we will prove Proposition 1151 below which shows 
in (i) how the constant C beyond all orders is associated to a trun- 
cated solution y(z) of Pi for arg(z) = ir (formula (|109j0 and gives 
the position of one array of poles z n of the solution associated to C 
(formula (jllUjO . and in (ii) provides uniform asymptotic expansion 
to all orders of this solution in a sector centered on arg(z) = ir and 
one array of poles (except for small neighborhoods of these poles) in 
formula (fTH|) . 

Proposition 15 (i) Let y be a solution of j!U4\) such that y(z) ~ 
for large z with arg(z) = ir. For any £ (ir,ir + |tt) the 
following limit determines the constant C (which does not depend on 
4> in this range) in the transseries y of y: 




i^^'IV^'- £ l^ l=c (109) 

arg(z)=</> \ k<\x[z)\ 

(Note that the constants yo-k do not depend on C). With this def- 
inition, if C 7^ 0, the function y has poles near the antistokes line 
arg(z) = 7r + g7r at all points z n , where, for large n 



24 V n V 8 4vr 600vr 2 j 



+ <>Wl .11.., 

n s 



with L n = jrjj In [jj^-Tij, or, more compactly, 
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e(^) = 12+ ( _ 24 ^ )5/4 +0(z- 5/i ) (z n ^oo) (111) 

(ii) Let e £ M + and define 

Z = {z: arg(z) > ^; |£(z)| < 1/e; - 12| > e} 

5 

fi/ie region starts at the antistokes line arg(z) = |7r and extending 

slightly beyond the next antistokes line, arg(z) = \tt). If y ~ y/—z/6 
as \z\ — ► oo, arg(z) = 7r, then for z £ Z we have 



( x 1 f 30*g t (Q \ 

^ 8V6(-z)5/2 ^ (-24z) 5k /*J 

(\z\ oo, z£2) (112) 

Tne functions are rational, and Ho(£) = ^(^/12 — 1)~ 2 . T/ie ex- 
pansion holds uniformly in the sector arg(z) £ (3/5,7/5) 
and a/so on one of its sides, where Hq becomes dominant, down to 
an o(l) if z is large distance of the actual poles of y. 

Proof. We prove the corresponding statements for the normal 
form (|1U8|) . To go back to the variables of ()1U4[) mere substitutions 
are needed, which we omit. 

Most of Proposition El is a direct consequence of Theorems 1 and 
2. For the one-parameter family of solutions which are small in the 
right half plane we then have 




oo 

/i- J>- fe F fc (£(x)) (113) 

k=0 

As in the first example we find H^ by substituting 1)1 13j) in l|108j) . 
The equation of Hq is 

£ 2 H'o + £ H 'o = H o + 

The general solution of this equation are the Weierstrass elliptic func- 
tions of ln£, as expected from the general knowledge of the asymp- 
totic behavior of the Painleve solutions (see [H])- For our special 
initial condition, Hq analytic at zero and Hq{£) = £(1 + o(l)), the 
solution is a degenerate elliptic function, namely, 

m) = (£/J- l)^ 
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Next, one of the two free constants in the general solution H\ is de- 
termined by the condition of analyticity at zero of H± (this constant 
multiplies terms in ln£). It is interesting to note that the remain- 
ing constant is only determined in the next step, when solving the 
equation for H2I This pattern is typical (see ^6.9|) . Continuing this 
procedure we obtain successively: 



Hi = (216 £ + 210 e + 3 e - ^ e 4 J (£ - 12)- 3 (114) 

h 2 - (uwe + 52 38{2 - f e - ^ + ±? + 5^5) « - ur* 

(115) 

We omit the straightforward but quite lengthy inductive proof 
that all Hk are rational functions of £. The reason the calculation is 
tedious is that this property holds for (|1U8|) but not for its generic 
perturbations, and the last potential obstruction to rationality, suc- 
cessfully overcome by (jlU8|) . is at k = 6. On the positive side, these 
calculations are algorithmic and are very easy to carry out with the 
aid of a symbolic language program. 

In the same way as in Example 1 one can show that the cor- 
responding singularities of h are double poles: all the terms of the 
corresponding asymptotic expansion of 1/h are analytic near the sin- 
gularity of h\ All this is again straightforward, and lengthy because 
of the potential obstruction at k = 6. We prefer to rely on an existing 
direct proof, see [TS| . 

Let £ s correspond to a zero of 1/h. To leading order, £ s = 12, by 
Theorem |2] (hi) . To find the next order in the expansion of £ s one 
substitutes £ s = 12 + A/x + 0(x~ 2 ), to obtain 

„ , . (A - 109/10) 2 nl , . 
Vfr(e.) = - 12 3 X 2 +0(1/X S ) 

whence A = 109/10 (because 1/h is analytic at £ s ) and we have 

£ s = 12 + ^ + 0(*- 2 ) (116) 
lUx 

Given a solution h, its constant C in £ for which ()113j) holds can 
be calculated from asymptotic information in any direction above the 
real line by near least term truncation, namely 



C = Jim exp(^)x 1 /2 h{x) < 1 1 7) 

nrg(x)=cf> \ k<\x\ 
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Fig. 3. Poles of igQH for C = -12 (o) and C = 12 (+), 
light circles are on the second line of poles for to C = - 



calculated via l|116jl . The 
12. 



(this is a particular case of much more general formulas Jl]) where 

2fc>o ho,k%~ k is the common asymptotic series of all solutions of 
H1(J8|) which are small in the right half plane. 

□ 

General comments. 1. The expansion scales, x and x~ l / 2 e~ x 
are crucial. Only for this choice one obtains an expansion which is 
valid both in Strans and near poles of (|l()8|h For instance, the more 
general second scale x a e~ x introduces logarithmic singularities in Hj, 
except when a E — ^ + 7L. With these logarithmic terms, the two 
scale expansion would only be valid in an 0(1) region in x, what 
is sometimes called a "patch at infinity", instead of more than a 
sector. Also, a E — ^ — N introduces obligatory singularities at £ = 
precluding the validity of the expansion in Strans- The case a E 
— \ + N produces instead an expansion valid in Strans but not near 
poles. Indeed, the substitution h(x) = g(x)/x n , n E N has the effect 
of changing a to a + n in the normal form. This in turn amounts 
to restricting the analysis to a region far away from the poles, and 
then all Hj will be entire. In general it is useful thus to make (by 
substitutions in ©) a = a minimal compatible with the assumptions 
(al) and (a2), as this ensures the widest region of analysis. 

2. The pole structure can be explored beyond the first array, 
in much of the same way: For large £ induction shows that H n ~ 
Const n .^ n , suggesting a reexpansion for large £ in the form 
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h „ £ £ ^F iZ ; £2 = C^kx- 1 = CCWi- 3 / 2 e- (118) 

By the same techniques it can be shown that 1)118)1 holds and, by 

matching with 1)113(1 at £2 ~ x ~ 2 / 3 , we get H [ q ] = H with CM = 
— 1/60. Hence, if x s belongs to the first line of poles, i.e. £(x s ) = £ s 
cf. (|116|) . the second line of poles is given by the condition 

x~ 3/2 e~ xi = -60 • 12C 
i.e., it is situated at a logarithmic distance of the first one: 

x\ — x a = — liix s + (2n + l)m — ln(60) + o(l) 

(see Fig. 15.2)1 . Similarly, on finds x s ^ and in general x s>n . The second 
scale for the n— th array is x~ n ~ 1 ^ 2 e~ x . 

The expansion (|113J) can be however matched directly to an adi- 
abatic invariant-like expansion valid throughout the sector where h 
has poles, similar to the one in [SO]- I n this language, the successive 
expansions of the form 1)118)1 pertain to the separatrix crossing region. 
We will not pursue this issue here. 



5.3. Example 3: The Painleve equation P2 
This equation reads: 

y " = 2y 3 + xy + a (119) 

(Incidentally, this example also shows that for a given equation dis- 
tinct solution manifolds associated to distinct asymptotic behaviors 
may lead to different normalizations.) After the change of variables 

x=(3t/2) 2/3 ; y(x) = x~ l (th(t) - a) 
one obtains the normal form equation 

„ ti ( 24a 2 + 1\, 8,0 8a, 9 8(a 3 - a) , . 
h" + 1 + = — I h h H h 2 + — 5 — - = (120 

and 

Ax = 1, a x = -1/2; £ = ^; + & = F + ^F 3 

The initial condition is (always): Fq analytic at and Fq{0) = 1. This 
implies 
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Distinct normalizations (and sets of solutions) are provided by 

x = (At) 2 / 3 ; y(x) = (At) 1 / 3 (w(t) - B + JjL) 
if A 2 = -9/8, B 2 = -1/2. In this case, 



„ w' ( Wa l-6a 2 ^ 
w H \-w [ 1 + —. — — | w 



t V tA 9t 2 



w ~^a) w2 + w * + W + 6 ° 2) " rla{a2 ~ 4)) (121) 



so that 



implying 



1 3Ba 
Ai = 1, ai = — 

2 2i 



e 2 ^' + - F = 35F 2 - F 3 



and, with the same initial condition as above, we now have 

2j(l + BQ 

The first normalization applies for the manifold of solutions such 
that y ~ — — (for a = y is exponentially small and behaves like an 

Airy function) while the second one corresponds to y ~ — B — ^x~ 3 / 2 . 
6. Appendix 

6.1. Some results in classical asymptotics 

The notations and assumptions are those of 32 

Theorem 16 Let y(x) be a solution of (0) satisfying |3j) on a direc- 
tion d which is not an antistokes line. 

Let S be the open sector bounded by two consecutive antistokes 
lines which contains d. 

Then 

(i) for any d' C S the solution y(x) is analytic on d' for x large 
enough, and tends to along d' . Also 

(ii) 0) holds on d' . 
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These facts follow from the proof of Theorem 12.1 of (f° r more 
general contexts see also the proofs of |26| . |27j ) and from the proof 
(of a similar theorem) presented in jSJ]. Unfortunately, (i) and (ii) 
of Theorem ^] were not formulated in these references as results in 
their own right. They also follow from the more general results of 
jl.Sj . but their essence is of a classical asymptotics nature and the 
ideas of exponential asymptotics are not really needed. (To compare 
the results obtained using classical versus exponential asymptotics 
approaches see Theorem 1191 and the Remark following it, H6.2I ) We 
therefore include here a self-contained proof of Theorem ll6l The iter- 
ation argument of [SJ is set up as iterations of contractive operators 
on appropriate Banach spaces of analytic functions. 

Proof of (i) 

Setting. Fix rj > and let S v C S be the open subsector whose 
bounding directions form an angle r\ with the boundary of S. We 
assume r\ is small enough, so that d C S v . 

Let xo £ d and let D = xq + £„. It will be shown that if \xq\ is 
large enough, then the solution y(x) satisfying (jHJ) is analytic in D, 
and tends to as |x| — ► oo, x G D. 

Denote yj(x ) = xp-. 

Note that for each index j, K(Ajx) has the same sign for all x E D 
(there are no antistokes lines in S^). Divide the coordinates of y into 
the two sets I + , J_ 



Integral equations. Equation © can be written in the integral form 



where the paths of integration IJj(x) C D are: the segment [xo,x] if 
j £ /+ and the half- line from oo to x, along the direction of x — Xq 
for j e 

Since the solution y(x) goes to along d we see that in (|123(l its 
constants of integration a,j are 

a,j = Q for j E I- , a,j = yj(xo)x" aj e X: > xo for j e 1+ (124) 

By assumption g(a; _1 ,y) is analytic at (0,0), say for {x^ 1 < r 
and |y| < p2, and satisfies |g(a? — 1 , y) | < const (|x| -2 + |y| 2 ) (see ©. 

Let Bo be the Banach space of functions y(x) analytic on D and 
continuous on D (with the sup norm) . Let T be the closed subset of 



I ± = {j = l... n ; n(Xjx) £l± , x € D} 



(122) 
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functions y G Bo with ||y|| < p (where p > will be chosen small 
enough) and satisfying Uj(xo) = y® for j G I + and Uj{xo) = for 

Relations ()123j) . ()124|) can be viewed as an equation y = ip + J{y) 
on T (if \xq\ > r _1 and p < p-i). For p small we show that if y G T 
then if) + J{y) G T and the fact that J is a contraction on T . It will 
follow that the integral equation has a unique solution thus proving 

(0- 

Lemma. Let 

|xo| > max \$t-Otj\ (|A,-| sinn) -1 max{l + \/2 , (V2 sinr/)^ 1 } 

j=l...n 

If : (i) j £ 1+ an d x (t) = x + — ^o) ; t £ [0, 1], or 
(mJ j G 2_ and x(t) = x + i(x — xo) , t > 
then 

1 x 



x(i) 



e £R[A*(*(t)-sO] < x 



(125) 



The proof of this lemma is straightforward (the left side of l|125jl 
is increasing in i in case (i) and decreasing in case (h)). The following 
estimates can be used: | cos arg[Aj(x — xo)]| > sin 7/ > 0, cos[arg(x — 
xo) — argxo] > — cos(2n) > —1, and for (i) R[Xj(x(t) — x)] = —(1 — 
t)\Xj(x — xo)\ cosarg[Aj(x — xo)] < —(1 — t)\Xj(x — xq)\ sin?? , while 
for (ii) $l[\j(x(t) — x)] = t| Aj(x — xo)| cos arg[Aj(x — xo)] < —t\Xj(x — 
xo)| sin 77 

The set T is invariant under iterations 
Let y G T . For j G I + 



\i>j{x) + Jj(y){x)\ < const \y 



X 

x 



jtlXjixo-x)] 



+ const |x — xq| 





X 


10 


x(t) 



JtMxW-x)] (^(t)!^ + || y ||2) dt 



and using (|125|) and that |x(i)| > |xo|(l — cos 2 (2n)) the last term is 
bounded by 



const |y°| + const |x— xq| (|xq| 2 



+ l|y|| 2 ) f 

Jo 



e 2 



i(l-t)|A f (x-x )\siav fa 



< const \yA + const (|xo| 2 + ||y|| 2 ) 
Similar estimates hold for j 6 L, hence (for some K > 0) 

U + J{y)\\ <^(|y°| + ko|- 2 + ||y|| 2 ) 
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Let p be small, such that p < (3K) . Then if |y°| < p(3K) , 
and |x |~ 2 < p{2>K)' 1 we have if> + J(y) £ T. 

Contraction. Let y,y' € T. Writing gj(x ,y) = gj i o(x~ 1 ) + 
Y^k=i...n9j,k( x ~ 1 ^y)yk witn 9j,k analytic for < r and |y| < p 2 

and gj ; k = 0(x~ 2 ) + 0(|y|) for k > 1, and g^o = 0(x -2 ) then 
\gj(x^ 1 ,y) — gj(x^ 1 ,y')\ < const(\x\~ 2 + p)\y — y' \ so that, with esti- 
mates similar to the above, we get || ^7(y) — ^(y')!! < const(|xo|~ 2 + 
p)\y — y'\. For small p and large xq the operator J is a contraction 
on JF, and part (i) of Theorem El is proved. 

Remark. In the estimates above the smaller r\ (i.e. the closer x to 
an antistokes line) the larger xq must be. This is closely related to the 
fact (which is the object of the present paper) that solutions (which 
are analytic in a "sector" — more precisely, in a region described in 
Theorem Hd]| develop (generically) singularities on the edges of this 
"sector" . 

Proof of (ii) 

Let 4>{x) be a solution of © satisfying @ — which is known to 
exist jSJ. Let y(x) be a solution satisfying (|HJ). Let u(x) = y(x) — 
4>{x). It is enough to show that u(x) = 0(|x| _r ) for all r € NU {0}. 

The function u(x) has limit along d and satisfies 

u = -Au + iiu + h(x _1 ,u) (126) 

x 

where 

h (x-\u) = g (x-\ u + (j>(x)) - g (x-\ <f>{x)) (127) 
As in the proof of (i) we write (|126|) in integral form (similar to 

Uj(x) = x a: > e~ x :i x aj+x a i e~ x i x I x^ aj e X: > Xl pj u(xi)) dx\ 

J n 3 {x) 

= ^ j (x)+J j (u)(x) , j = l...n (128) 

where the paths of integration are those of (|123|> and the constants 
aj satisfy the analogue of (|124j) 

o,- = for j'GL , Oj- = Uj(x )x~ aj e X J X0 for j G J + (129) 



Formation of singularities in nonlinear differential systems 



49 



where u = y(x ) - 4>(x ). 

Let D be as in (i), where cj) is analytic. Consider the Banach space 
B r of functions u(x) analytic on D, continuous on D, with the norm 
||u|| = sup X £ D |x r u(x)| (see also (50]). 

Let T be the closed subset of functions u 6 B r with ||u|| < p 
(where p > will be chosen small enough) and satisfying u(xo) = Un. 

Note that <j>(x) = 0(x" 2 ) (cf. ©) hence \</)(x)\ < M\x\~ 2 (for 
x £ D and xo large enough). Then since g(x _1 ,y) was assumed 
0{x~ 2 ) + 0(|y| 2 ) it follows that for x £ D (cf. (tTTTlO 

|h (a; - , u) | < const (|x|~ 2 |u| + |u| 2 ) 

(for zo large enough, so that |</>(x)| < p/2 and for |u| < p/2). 

The same estimates as in the proof of (i) (the only difference 
being that the otj of (i) should be replaced here by r + ctj) show 
that equation u = ifi + ^7( u ) has a unique solution in T if |xo| is 
large enough (depending on r, as expected). Hence ||u|| r < oo which 
concludes the proof of (ii). 



6. 2. Summary of some results in liffi 

This subsection contains details on results of ^3] cited, referred to, 
or relevant for the present paper. A simple consequence of a Lemma 
in JH] (needed for the present paper) is formulated and proved at the 
end of this section (Theorem IT§|) . 

Since the Theorems, Lemmas and some formulas cited in this sec- 
tion are from ^3]) to avoid repetition we will follow by a * sign any 
result cited from [T3] . 

The setting is the same as in the present paper: the equation 
studied is 

y , = -yiy-- J By + g( 2 r 1 ,y) (130) 

x 

(same as with B = —A) having transseries solutions 

y(s) = C k e~ A - ka; x-^ k s k (x) 

ke(Nu{o}) n 

^ £ C k e- A - k ^- M - k y k (x) (131) 

ke(Nu{o})» 

(same as (@J), ©, d} for 0j = -<Xj, M j = + 1, j = l...ri). 
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Since the association between actual and formal solutions depends 
on directions (Stokes phenomena) a sector in the complex a>plane 
is chosen as follows. Fix some non-empty open sector S' C C and 
consider those transseries (jl31|) valid in S' (as explained in ^2.1(1 . 
Some constants Ci, C n may be required to be zero in S', say Cj = 
for j = rti+1, n (with m > 0). Let Strans be the (non-empty, open) 
maximal sector of validity of any transseries (jlHlj) with Cj = for 
j = n\ + 1, n (see ©). 

To simplify the notations it can be assumed (after trivial changes 
of coordinates) that Ai = 1 (see also iJ2.2|) . 

I. The construction of actual solutions associated to transseries 
solutions valid in Strans is done in ^H] using a generalized Borel 
summation as follows. 

Denote by Y(p) the formal inverse Laplace transform of y(x) (i.e. 
Y(p) = J a + lo ° e px y(x) dx, its convergence following from sub- 

sequent analysis). Using the usual properties of the inverse Laplace 
transform (e.g. the transform of y'(x) is —pY(p), multiplication is 
transformed into convolution, etc.) the differential equation (|13Uj) is 
transformed into a convolution equation (eq. (1.13)*). 

The Stokes lines in the Borel plane (p-plane) are defined as the 
complex conjugates of the Stokes lines in the direct (x) space. For 
linear equations the Stokes lines in the p-plane are dj t o = A,-R+, 
j = l...n. For nonlinear equations, there also are other Stokes lines 
(which play a role only in the higher order terms of the transseries, 
with |k| > 2, hence they are "transparent" in the linear case) namely 
d j;k = (Aj - k • \)R + with j = l...n, k e (N U {0}) n (note that 
Pj,k G c?j,k cf.JSSI))- 

Once the sector Strans is fixed, there are only finitely many lines 
dj^r in this sector. 

In any proper subsector of any of the n sectors formed by the 
Stokes lines dj ; o the convolution equation has a unique solution Yo(p) 
which is analytic at p = (Lemma 16*). Yo(p) is in fact the Borel 
transform of the asymptotic series yo(a^) (see (j3J)) (i.e. 
Yo(p) = Ylr yo,rP r / r 0- Yo has singularities on the Stokes lines 
at p G AjZ + , j = l...n (hence the classical Laplace transform cannot 
be taken on M + and the classical Borel sum of yo does not exist). 

Denote by Yq the analytic continuation of Yo on directions above 
d\ t o = M + (but below the neighboring Stokes line), respectively by 
Yq for the continuation below M+; they exist see Lemma 16*. It is 
shown that as p approaches R+ from above (or below) Y^(p) (respec- 
tively, Yq(p)) tends to a distribution on IR+ in an adequate space of 
distributions — the staircase distributions, introduced in ^1 (Lemma 
16*). In general, the two distributions are different. (Of course, a 
similar picture holds at any other Stokes line.) 
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Higher order functions Yk, |k| > 1 are then constructed (Lemma 
20*) by solving the convolution equation on the Stokes lines. 

Consider for example the line M+ = di,o- Fix a solution Yo of the 
convolution equation in the space of staircase distributions on K+. 
The construction of the higher order Yk's with kj = if j > n\ + l 10 
is done in the proof of Lemma 20* as follows. In view of the sought- 
for expansion (|10[). after introducing it in H130j) and identifying the 
coefficients one obtains (a recursive system of) differential equations 
for y^; formal inverse Laplace transform yields (a recursive system of) 
convolution equations for Yk- Once a staircase distribution solution 
Yo on R + is chosen the general solution of this system with regularity 
(|34"|) depends on n\ free constants C±, C ni , in the form C k Yk- 
Outside the Stokes line M+ the solutions Yk are, in fact, analytic 
up to the nearest direction (of positive argument ip + , respectively 
negative argument tp-) which is either a Stokes line d,k which lies in 
the right half-plane (i.e. the half-plane orthogonal to d\ t o = M+) — 
where some of the constants Cj may change, or is an antistokes line 
associated to Ai = 1, i.e. iM+ or iM_ — where the transseries is no 
longer defined 11 — (Lemma 20* (i)-(iv)). 

It is interesting to note that the Yk multiplied by Stokes con- 
stants are generated as differences between different branches of Yo 
(Theorem 4*, Proposition 23*). 

Then (a generalized) Laplace transform is applied to Yk (in the 
space of staircase distributions on the Stokes line under consid- 
eration) yielding (analytic) functions yk = £Yk- 

The last step in the summation of transseries is showing that the 
sum IjlOjl converges (Lemma 20* (v); more details are found in the 
proof of Lemma IT71 of this section). 

The reconstruction of a solution from a transseries is concluded 
showing that the function y(x) obtained as the sum of (jlOj) is a solu- 
tion of the differential equation (jl.SOj) — which follows easily because of 
appropriate convergence and since all functions have been constructed 
from formal objects satisfying the equation (Lemma 20* (v)). 

The correspondence between transseries (i.e. the constants C) and 
actual solutions given by the summation of Lemma 20* is not unique. 
This is due to the non-uniqueness of staircase distribution solutions 
Yo on R + (on which the higher order Yk depend): there is a one 
parameter family of solutions Yq (among which Y^ for a = 0, 1) — 
they are special averages of analytic continuations of the germ of 
analytic function Yq at p = 0. 



Other Yk's are not needed since the corresponding cannot be present in a 
transseries on the fixed Strans- 

11 Yk may be analytic beyond the antistokes lines ±iR+ but the analysis in |13| 
stops there. 
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2. Conversely, any solution of (|130|) satisfying on a direction 
d in the right half-plane is the sum of a transseries, which is unique 
once q is fixed (Theorem 3*(iii)). 

The proof is done in the following steps (Propositions 24* and 
25*). There exists a solution yo a (x) asymptotic to the transseries 
yo (i.e. y with all Cj = 0) as constructed at part 1. above. It is 
also shown that any two solutions y 1,2 (x) with the same asymptotic 
expansion y 1,2 ~ yo on d differ by exponentially small terms: y 1 (x) — 
y 2 (x) = J2j=i n^i e ~ XjXx ~^ j ( e j + on d- Thus the difference 

y (x) — yg (x) fixes the constants Cj = Cj (a, d) , which are then used to 
construct a solution yt(x) from the transseries y with these constants 
(using part 1. above). The last step is showing that y(x) = yt(x). 

3. It has been thus established that (given a direction d, and a 
parameter a) there is a one-to-one correspondence between trans- 
series solutions and actual solutions of (|130j) . The correspondence is 
built using a (family of) generalized Borel summation(s) CB a on d. 
It is shown that the operator CB a is compatible with all algebraic 
operations (performed on transseries, respectively functions on d). 

The Stokes phenomenon is analyzed in Theorems 4* and 5*. 

4. We now state the maximal domain of analyticity 12 of y(x) 
implied by formula (2.41)* of Lemma 20*. 

Lemma 17 Let y(x) be a solution of U3(J\) satisfying on a di- 
rection d above (but close enough to) M + . 

Let CJ , j = 1, ...,n be the constants such that y(x) is represented 
on d as 

y = £Y + ^ (C~) k e- k - Xx x' k ^CY k - (132) 
|k|>0 

(see Theorem 3* (Hi)). 

Then for any e, 5 > there is x\ > such that y(x) is analytic 
on the domain 

— ( _ TT _ TT 

D an = 1 x ; M > x i , arg{x) G [-tp -- + e,-ip + - - e] , 

<,5- 1 , j = l,...,n} (133) 

The constants e, 5 are the same for all solutions of \13CJ\) with 
transseries valid in the same sector St ran s as y- 

The Proof follows immediately from Lemma 20* (v), but we pro- 
vide the details. 

All functions Yi c ~(-e 4(?i ) are Laplace transformable in the space of 
staircase distributions on M + with exponential weight e~ vp (for v > 
large enough) if <f> G (ip_,ip + ) (see Lemma 20* (v)). 

12 This domain is improperly stated in 
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For <j) £ (ip-,0) the Yk _ (pe 1 ^), (p > 0) are analytic (Lemma 
20* (iii) ) , so the Laplace transform in the space of staircase distri- 
butions coincides with the classical Laplace transform (see Lemma 
6* and the classical properties of the Laplace transform and of the 
spaces L\^) where the Laplace transform is defined as 

CF{x) = [ e~ px F(p) dp , for x £ d (134) 
Jd 

(Note that in (fr31|) p £ d and x £ d so that 5i(px) > 0.) 

Remark 18 If F is analytic on a direction d = e~ tv M + (rj £ ~S.) and 

\\ F \\u = [ e~ up \F(p)\ dp < oo 
Jd 

then its Laplace transform \1'34\) ^ s analytic for x £ d = e ir, M + , |x| > 
v . 

Furthermore, for any e > CF has analytic continuation on the 
domain argx £ [i] — ^ + e,rj + ^ — e] = I and \x\ > ^(sine) -1 and 
CF satisfies 

\CF(x)\ < \\F\\ U (135) 

The proof is immediate, noting that the path of integration of 
()134|) can be rotated to any other direction in the interval / — since 
/ is so that $i(xp) > 0. 

To prove Lemma^llet e be small, positive and let eft £ [tp- +e, — e]. 
For any 5 > there is v > so that the functions (pe^), (p > 0) 
satisfy ||Yi c ~(-e l<?!> )|| i y < <5' k ' for all multi-indices k considered (i.e. 
|k| > 0, kj = for j > n\) see Proposition 22* (ii). Using Remark UHl 
and (|132j) the result is immediate. □ 

Remark 

Similarly to Lemma El there is x\ > such that y(x) is analytic 
on D+ n = D+ n (e,5,C+) where 

D tn = \ x \ \x\> xi, arg(x) £ [-ip + - - + e, -ip + + - - e] and 

<6- 1 j = l,...,n] (136) 

so that the domain of analyticity of a solution y(x) includes domains 
of the form D~ n U D+ n . 

Theorem 19 Let y(x) be a solution of hl'JU]) satisfying (0) on d = 
M + . Let e > be small. 

There exists 5, R > such that y(x) is analytic (at least) on 

San = S an (y(x); e) = S+ U S~ (137) 



54 



O. Costin, R. D. Costin 



where 



Sf = ^x;\x\>R, arg(x) G [— — =F e, — =F e] and 

< r 1 for j = l,...,n} (138) 



The constant 5 is the same for all solutions of \13U\) with trans- 
series valid in the same sector Strans as y- (However, R does depend 
on C.J 



Proof 

From Lemma IT71 using the expansion of Y in terms of Yk~, re- 
spectively Yk + it follows that y(x) is analytic for \x\ > x±, argx G 



r+e,— ip +f — e] = / and 



Cfe- 



X i x X 



Mi 



< 5- 



• J 



1, 



■ n. 



Since Mj = — l$t(3j\ + 1 and the sector / is larger than the sector 
where all exponentials in the transseries of y(x) are bounded, the 
result follows. 



Remark 

Fix a direction d (not an antistokes line) and consider solutions 
satisfying (JHJ). It is interesting to compare the result on the domain 
of analyticity of these solutions as given by Theorem El (obtained 
using results of exponential asymptotics) to the result of Theorem 
(obtained in a classical setting). There are special families of 
solutions for which the sector of analyticity given by Theorem 
is, in fact, larger than the sector between two consecutive antistokes 
lines (i.e. solutions having the corresponding C,- zero). For Painleve 
PI equation these special solutions are called triply truncated. 

Convention 

The Borel summation used in the present paper is CB = CBi . 



6.3. Proof of {££p 

Proof. This follows from the definition £ = C\x ai e~ x and from the 
asymptotic behavior of the functional inverse W of se s (see e.g. |12|). 
For large t > the branch of W which is real has the expansion 
(convergent, as it is not difficult to show) 



ttt, xi i . lnln * A(lnlnt) 2 -lnlnt 

W (t) = hit- In lilt + + ^ TT-rr= + 

mi (lntj^ 
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6.4- Points on 7^ 



have the same magnitude 



Let R be large, so that 



P3 + 



3KB 
R 



< P2 



(139) 



Lemma 20 There is a small enough neighborhood MjO 0/7^ so that 
any x' , x" E A/lo satisfy 



7N(to)\ - \j N (t)\ = 9 (ln(r(t)) - ln(r(t ))) + o(l) (N -> 00) 



Therefore lirri2v->oo TaKO/t^KO = 1 uniformly for G [£o>l]- So 
for TV large \x/x'\ < 3/2 for all x, x' on 7^ between xq and a, so in a 
small enough neighborhood |x/a;'| < 1/2 which proves the Lemma. 

6.5. Special estimates 

We show that the second argument of g in (|65|) has absolute value 
less than p 2 (cf. (J23J)). 

We need the following lemma 

Lemma 21 Let A, en > 0. There exist Aq, k > sue/t i/iai 



if A > Aq and m\A~ m < c . 

Before giving the proof of the lemma, we show how (j!41|) is used. 
Let R be large, and cq small, so that 




(140) 



Proof 

Using Q24|) for t G [to, 1] we get the uniform estimate 




(141) 




(142) 



In view of Theorem for s G A/lo 
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and in view of Lemma QUI the last term is bounded by 

fc=i v 1 1 ' 



Using 1)143|) , Lemma |^ and the bound p% on Fo and ()142|) we 

have 

|FH ( X )| < |F (*(*)) \ + K K ( 2 ^ + c ) < < P2 



Finally, if \S(x)\ < on Jfy then 



Proof of Lemma 1211 

Estimates like (|141j) are common in proofs using least term trun- 
cation of factorially divergent series (see e.g. Jl]). The proof of ()141j) 
is included here for completeness. 

The series 

k>l 

is divergent. Its terms decrease for k < A and increase for k > A; the 
term with k = [A\ (or k = [A\ + 1) is called the least term (see e.g. 
0)- 

Case I: m < A 

In this case the terms in the l.h.s. of (|141|) are decreasing, hence 
kl 1 2 , ,3 

Ei^i+i2(--i)<i 

fc=l 

Case //: A < m < eA/2 

The terms in the l.h.s. of 1)141)1 are increasing for A < k < m, not 
exceeding the second term: m\/A m < 21/ A 2 . 

Indeed, this is a simple estimate using Stirling's formula and the 
fact that the function F(A) = A^ 2 a aeA is decreasing for A > Aq (if 
Aq is large enough). 

Then as in Case I 

^ k\ 1 2 . . e + 1 

Y.T^A + l¥ {m - l)< — 

k=l 



Case II: eA/2 < m 

Denote p = [A\ and q = [^-p - \\ . 
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Write 



fc=i 



where 



Jfc! 



Sin + Sn + Si 



k\ 



k\ 



Sni = ^2jk ' Yl -jk > Si= E # 



k=l 



k=p+l 



A k 



k=q+l 



A k 



and estimate each sum separately. 
To estimate Si 



S < ^ < m ' 

k=q+l F F 



A;=g+1 



(5 + 2) m ~ fc 



and since + 2) < 2/(^4 + 1) < 1 this is less than 

m! yl + 1 m! e + 2 
< p™A^T < p™e-2 

To estimate 5// note that 



^4 \ m ~<? 

< (q - p) -7- < -r^O? - p) ( — , 



q\ to! 



and since < po < 1 if A > Aq (for An large enough) 



< j^{q-p)pq 9 < '—r>\. 



ml 



Finally 

The result of Lemma |2*T1 follows. 



5m< i + Jj (p-1)< i 



6. 6. Proof of Proposition ^] 

A consolidation of one of the norms in Example (3a) in ^3] is first 
needed. For convenience we repeat that part. The notations are those 
in pig. _ 

"(3a) For > and fa ^ fa, let T^UV) = {/ : f(p) = 

p@F(p)}, where F is analytic in the interior of £ U V and continuous 
in its closure. We use the family of (equivalent) norms 
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u ,p = \r(P + l)\K sup |e-^/(p)| 
sesuv 



(145) 



It is clear that convergence of / in \\\\ v fi implies uniform convergence 
of F on compact sets in £■ U V (for p near zero, this follows from 
Cauchy's formula). Tp are thus Banach spaces and focusing spaces 
in || 11^/3 by (j!45j) . The spaces {Tp}p are isomorphic to each-other. 
Convolution is defined as 



~01~02 



~\fi * h)(p) = C t^Ftipt^l - tf 2 F 2 (p(l - t))dt = F(p) 
Jo 

(146) 



where F is manifestly analytic, and the application 
(•*•): Tp 1 xT ft « T Pl+ p 2+ i 

is continuous: 



(147) 



||/l * f2\\u,0 1 +/3 2 +l 



\K sup 


z~ up [ 


V 


Jo 



s Pl Fi(s){p - sf 2 F 2 {p - s)ds 



r(A + fa + 2) 

< =z sup 



A'F 1 (s)e- I/s s /31 KF 2 (p - s)e- u( - p ~ s \p - s] 



02 



r(Pi + 1) 



d|s| 



m + 1) 

^WhUpAhWv^ (148) 



Estimating the norm of Y k exactly as in ^3] but using this inequality 
instead of (2.8) of ^3] we get that 

|Y k (p)| < |r(-k-a / ) _1 |^' k| e 1/obl 

and thus, using straightforward Cauchy estimates in d ai of deriva- 
tives, is proved. □ 
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6. 7. Note on normalization of the on 

The reference |13| uses a transformation that makes f3j < (ay > 
in the present notation). To determine the singularities of y it is now 
important to make ctj as small as possible, as explained in Comment 
1, N5.21 In some cases we must then allow for m < in |13| . Eq. 
(2.43). This does not affect the estimates (2.44) through (2.46) in the 
space 7k.pi -i, the only one that relevant to the present paper. Minor 
modifications of the proof following Lemma 20 in ^3] are needed. For 
completeness we redo redo here the whole proof. 

For |k| > 1 with Wk := Yk and Rk := Tk, the functions Wk 
satisfy the equations 

(1 + Jk)W k = Q^Rk (149) 

with Qk := ( — A + P + k • A) (notice that for |k| > 1 and p E S' we 
have detQk(p) 7^ 0). 



(J k W)(p) := Q,, 1 ((s + m-k) jTw(«)ds 



v 



£ / W^-Djip - s)ds ] (150) 



o 



Proposition 22 i) For large v and constants K\ and K.2{v) inde- 
pendent o/k, with K^iy) = 0{y~ l ) we have \\Q%. || < t& and 

|| J k || < K 2 H (151) 

iij For Zarge v, the operators (1 + Jk) defined in 7>' mv , and also in 
7 ]s _pi_ 1 for |k| > 1 and in 7\ for |k| = 1 are simultaneously invertible. 
Given Yo and C, the Wk, |k| > 1 are uniquely determined. For any 
5 > there is a large enough v, so that 

||W k || < 5 |k| , k = 0,1,.. (152) 

(in the V mv topology, \15ty) hold uniformly in (j) G [ip- + e,0] and 
(f> 6 [0, -0+ — e] /or any smaZZ e > 0). 

□ 

Proof. 

(i) follows immediately from ()148|) . 

(ii) From Q149|) and (i) we get, for some K and j > 1 ||Wk|| < 
if||Rk||- We first show inductively that the Wk are bounded. Choos- 
ing a suitably large v{e) we can make max| k |<i ||Wk|| y < e for any 
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positive e (uniformly in (p). We show by induction that ||Wk||^ < e 
for all k. In the same way as in |13j we get 

HWkii, < ^iiR k |i, < y, k i^ 1 1 

l<k (i mp ) 
|k| 

< el k l ^2 ^2 ni d k l +s )2 s+ni < (Cie)l k l (153) 

s=0 

where C\ does not depend on e, k. Choosing e so that e < C-f 2 we 
have, for |k| > 2 (Cie)l k l < e completing the induction step. But as 
we now know that HW^H^ < e, the same inequalities (j!53j) show that 
in fact HWklli/ < (Cie)' k L Choosing e small enough, the first part of 
Proposition 1221 (ii) follows. □ 

6. 8. Proof of Lemma MIA (ii) 

1. Generically ho is not entire. Assume /in is analytic in a neighbor- 
hood of the disk Br = {|£| < R} and let Mr = sup| z | =R /in(z). We 
have 

sup {z~ l {h {z) - 1) = sup {z^ihoiz) - 1) < R~ 1 {M + 1) 

\z\<R \z\=R 
and thus, for some constants Cj we have 

M 2 < CiR 2 + C 2 + {C 3 R 2 + C 4 )M (154) 

whence 



M < C 5 R 2 + C 6 (155) 

If /in is entire it then follows that /in is a quadratic polynomial in 
£. But it is straightforward to check that does not, generically, 
admit quadratic solutions. Thus the radius of analyticity of /in is 
finite, say R , and 13 

sup {aT 1 (/*)(*) - 1} < RqHMr, + 1) and Mr, < C 5 R 2 + C 6 
\z\<R 

(156) 

13 An upper bound for Ro can be found by comparing /io"(0) with its estimate 
from Cauchy's formula and l|154|l 



Formation of singularities in nonlinear differential systems 



61 



2. /iq is uniformly continuous on Br . Indeed, if € £>_r we have 
by (jHSJ) and (11561 that /i 2 , is in fact Lipschitz: 

1^(0-^(01 < Connie (157) 

3. If £o £ 5^i? an d ^o(Cq) 7^ then £o is a regular point of eq. (|51)l 
and thus ho is analytic at £o- 

4. J/£ s G 0Br is a singular point of ho and — A^.7 1 + ^3 + d^ s ^ 
then £ s is a square root branch point of ho, i.e. ho(£) = — 4s) 1 / 2 ) 
where h± is analytic at zero. From parts 2 and 3 above, ho(£ s ) = 0. 
It is convenient to look at the equation for £(7in) derived from ((Ml) : 

= 2 -h . m = e 

dh (A 2 £~ 1 + d 1 + d 2 £)ho + (-As^- 1 + d 3 + d 4 

(158) 

whose unique solution is analytic near zero. The claim now follows 
by noting that f'(0) = and £"(0) = 2(-A 2 £: 1 + d ?J + d£ s ). 

5. We now restrict the analysis to a smaller but generic set of coeffi- 
cients. We denote by K s the following subset of parameters (see (|79j) 
and (|HH) ) 



K s = {(d,j) = (dj, 7j)j=i,..., n G C 2n : ho not entire and 

P(£) = -2A 2 + \2diC + (f 3 4 2 has distinct roots } (159) 

We show in parts 6 through 8 that if £ s G 9Bji is a singular point 
of ho and P(£ s ) = 0, then a generic small variation of (0,7) in i^T s 
makes P(£ s ) 0) and by part 4, £ s becomes a square root branch 
point of ho- 

We thus assume that P(£ s ) = 0. The substitution £ — £ s = i in 
(JBU) gives 

2/jo/i'o = B x h + tB 2 (160) 

where 



Bi = 2\ 2 {t + 6)" 1 - d x + (t + 6)c/ 2 ; # 2 = *tf% g rf 3 

and the roots of P(£) are £1 7^ £ s . 

We now study /io in the following smaller region. Let to be small 
on the segment [0, — £ s ]. Choose Z\ = {t : \t — to\ < \to\} C Br q — £ s 
to be a disk tangent at t = to i?/j — £ s which does not contain 
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the points £ s — £ and — £ s . Then /io is analytic in Z\ and continuous 

in Z\\{0) (while hn is continuous in Z\), and lim hn(t) = 0. We 

A3t-*0 

assume t = is a singularity of ho- 

6. There exists a sequence {t n } n in A with t n — > suc/i that 
lim^oo h (t n )/t n = L with 2L 2 - LBi(O) - B 2 (0) = 0. 

We first show that ho — * 0, then prove ho/t is bounded below and 
above, and finally that ho/t has a limit. 

(a) . We estimate 

.Mt = max |/in(s)l 
l--*l<l*l 

for t £ (0, to) from (|16U|) written as 

/»§(*) = f B 3 (s)h (s)ds + t 2 B 3 (t) 
Jo 

where B 3 ^ are analytic in A and continuous on A. Then 

M 2 < \t\M t max |B 3 | + \t 2 \ max |B 4 | 
A A 

and thus Ait < | i | , for some ifi > and all t £ [0, to]. 

(b) Cauchy's formula on the circle \s — t\ = \t\ implies |/i (t)| < K\. 

(c) Equation IjlbUj) written in the form t/ho = (2h' — B^B^ 1 
implies now \t/ho\ < < oo. In conclusion, 



t 

h 



£ [K 2 ,Ki] for all t E [0,t ] 



To conclude the proof of the statement at the beginning of this part, 
the function y = ho/t, which is analytic in A and continuous on 
A \ {0} satisfies the equation 

2ty' = -2y + B x + B 2 /y 

which can be written as 

i y' 



2t P (y)+tf(t,y) 



(161) 



with Po(y) = —2y + Bi(0) + B 2 (0)/y. Assume, to get a contradiction, 
that for some e > we had |P (y(*))l > e for t £ (0,ti] C (0,t }. 
Since 



+ th(t,y) (162) 



Po(y) + tf(t,y) P (y) 



Formation of singularities in nonlinear differential systems 



63 



with fi(t,y) bounded on (0, t{\ if t\ is small, we get by integrating 
(Q2TJ) on [t,ti] C (0,ii] 

i ln(iAi) = FoMt)) - Fo(y(t x )) + jf ai/00/l(*, 2/00)^ 

where i*o is a primitive of l/Po> anci thus Fo(y(t)) is bounded for £ E 
(0,ii). Hence the r.h.s. of (j!62|) is uniformly bounded for t E (0, t±), 
which is a contradiction, given the l.h.s. of ([162)) . 
7. With S(t) defined by h (t) = Lt(l + <5(t)) we have S(t) -> as 
i — > 0, i E (0, to]. We show this by a contraction argument. 5 satisfies 
the equation 



tf'= (-2+^-\s + B(t,S) (163) 



where 



#(M) 



gi(0) 
2L 



1 + 5 



+ t 



Bi(t) 



2Lf 



B 1 (0) + B 2 (t) 



B 2 (oy 



2LH{\ + 5) 



or in integral form, 



6(t) = (JS)(t) = (l-Y S(t n )+t b ' t s-^BfaSisfids (164) 

\ 71 / J t n 

where {t n } n is the sequence found in 6, such that S(t n ) — ► 0. We 
have two cases, according to whether is positive or negative 

(3i(£>i) = is nongeneric). 

(a) In the case M(bi) > we define J on the space A of analytic 
functions in B n and continuous in B n , where B n is the ball having 
[0,t n ] as a diameter. Let A r = {5 E A : ||<5||oo < r} and r n = 2|5(t n )|. 

(b) For 3?(&i) < we instead define J on the space A' of analytic 
functions in B' n and continuous in B' n , where B' n is the ball having 
[t n ,t n -%] as a diameter. Let „4.' r = {5 E .4' : ||£||oo — r } an d r n = 
2\5{t n )\+2(t n _ l -t n ). 

Lemma 23 ^aj If $t{b{) > 0, /or n /arye enough, J : A Tn i— ► A Tn 
is contractive. Therefore, as n is large, we have \5{t)\ < 2\5(t n )\ on 
[0,t n ] so that 8(t) ast^O in (0,t ]. 

(b) 7/5ft(6i) < 0, forn large enough, J : A', i— > A', is contractive. 
Therefore, as n is large, we have \8(t)\ < 2\5(t n )\ + £ n _i — t n on 
[t n ,t n —i] so that, again, 5(t) — ► as t — ► in (0, to]- 



64 



O. Costin, R. D. Costin 



Proof. A straightforward calculation. 

8. Now we bootstrap the information that 6(t) —* to sharpen the 
characterization of 5(t) for small t. 

Lemma 24 (a) If dt(bi) < then 5(t) is analytic in t at zero, and 
<5(0) = 0, thus 5{t) = 0{t). 

(b) If Sft(6i) > then 5(t) = 0(t bl ) + 0(t) for small t. 

Proof, (a) Since 6(t) — > we have 



S(t)=J (tzb 3 (S(tz),tz) +b 2 5 2 {zt)^z- bl - l dz (165) 

which for small t is manifestly contractive in a small sup ball in a 
space of analytic functions in a neighborhood of t = 0. 

(b) Fixing some n, from (|164j) and Lemma l28l we see that 

W)\< l -^1* + (t n - *)*P[|oo 

□ 

(In fact it is not difficult to show that in case (b), 5 can be written 
as an analytic function in the two variables t and t bl .) 
With u = we get the equation in variations 



uti + h u' = (2A 2 r X + a 2 )u + 1 (166) 

whence u is analytic as long as ho ^ 0. This equation being linear we 
can use classical Frobenius theory and, in the only interesting case 
R(6i) > we have u(£ ) = (Ai j2 - 2X 2 ^ 1 - aa)" 1 ^ 0. Thus an 
arbitrarily small variation of A 2 makes /io(£o) 7^ 0. 



6. 9. The recursive system for F m 

In applications it is usually more convenient to determine the func- 
tions F m recursively, from their differential equation. Formally the 
calculation is the following. 

The series F = J2m>o x ~~ m ^m{€) is a formal solution of @; substi- 
tution in the equation and identification of coefficients of x~ m yields 
the recursive system (JSHJ, (jBTj) . To determine the F m 's associated 
to y we first note that these functions are analytic at £ = (cf. 
Theorem P). Denoting by Fmj, j = l,..,n the components of F m , 
a simple calculation shows (|56*]l has a unique analytic solution satis- 
fying F 0)1 (0 = £ + 0(O and F 0ii (£) = C>(£ 2 ) for j = 2,...,n. For 
m = 1, there is a one parameter family of solutions of (|57|) having a 
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Taylor series at £ = 0, and they have the form Fi t i(£) = ci£ + 0(£ 2 ) 
and jFi 3 - (£) = 0(£ 2 ) for i = 2, ...,n. The parameter c% is determined 
from the condition that (|57[) has an analytic solution for m = 2. For 
this value of ci there is a one parameter family of solutions F2 ana- 
lytic at £ = and this new parameter is determined by analyzing the 
equation of F3. The procedure can be continued to any order in m, 
in the same way; in particular, the constant c m is only determined at 
step m + 1 from the condition of analyticity of F m +i. 



6.10. Sketch of a classical proof of Theorem^ 

It is also interesting to mention a direct, classical proof (i.e. not in- 
volving results of exponential asymptotics) of Theorem ^ (Since we 
do not rely on this more involved approach, we only give a brief out- 
line of this proof.) 

Having determined the initial conditions for F m as above, equa- 
tions Q56|) . (|57|) can be transformed to integral equations possessing 
unique analytic solutions F m for small £. 

To show (J2U|) let y(x) be a solution of (j2J such that y(x) — » in 
Strans- Denote 

/ N 
R N (x) = x N+1 y(x) - x- m F m (C(x)) 

V m=0 

(the remainder of y(x) with respect to the truncated expansion). 
Then R^r satisfies the differential equation 

R N = F N (x,R N ) (167) 
where 

F N (x, Rjy) = d y g(0, F )R N + F(x)R N + u N {x) + F^ N {x, R N ) 

(168) 

with F^ N {x,R N ) = 0(R 2 N ). 

Let RnJi j = 1, denote the components of Rn and Enj be 
the components of E^v- 

Equation (|167|) can be written in the integral form 



art 

dx 



N 



+ 



A + x~ l (A - (N + 1)1) 



R NJ (x) = e-^ x x N+1 -^ e^ s s- N - 1+ ^E N (s,R N (s))ds (169) 

J 

(for j = l,...,n). For an appropriate (rather delicate) choice of the 
initial points x® and of the contours of integration in ()169j) the in- 
tegral operators defined by the r.h.s. of (|169|) are contractive for N 
sufficiently large, hence (j!69|) has a unique analytic solution R^v- 
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1. Introduction 



In generic analytic nonlinear differential systems in the complex plane, 
we study the position and the type of singularities formed by solu- 
tions when an irregular singular point of the system is approached 
, along an antistokes direction 1 . Placing the singularity of the system 

{Sj ' at infinity we look at equations of the form y' = f (x , y) with f ana- 

lytic in a neighborhood of (0, 0), with genericity assumptions; x = oo 
is then a rank one singular point. We analyze the singularities of 
those solutions y(x) which tend to zero for x — * oo in some sectorial 
4^ ' region, on the edges of the maximal region (also described) with this 

property. 

After standard normalization of the differential system, it is shown 
that singularities occuring in antistokes directions are grouped in 
nearly periodical arrays of similar singularities as x — » oo, the loca- 
tion of the array depending on the solution while the (near-) period 
and type of singularity are determined by the form of the differential 
system. 

This regularity in type and position of movable singularities has 
been observed previously in various examples of nonlinear systems: 



^1" 



X 



Send offprint requests to: 

1 In the sense stemming from Stokes original papers and the one favored in ex- 
ponential asymptotics literature, Stokes lines are those where a small exponential 
is purely real; on an antistokes line the exponential becomes purely oscillatory. In 
some references these definitions are interchanged. 
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Painleve equations ((BUI, ESj- ^El) third order nonlinear equations 
(HH, |22j) nonintegrable Abel equations ([22], EH) 

among others. 

We show these features are rather universal and find a formalism to 
calculate them (asymptotically). 

When f is meromorphic and satisfies some estimates the singular- 
ities in the arrays are generically square root branch points. 

The mechanism of singularity formation is elucidated by expo- 
nential asymptotic analysis, which also provides a general and effec- 
tive calculation tool for determining the type and position of singu- 
larities. The present method generalizes that of ^Bl- The analysis 
yields two-scale asymptotic expansions of solutions, valid in a re- 
gion which includes the directions along which y — ► and extend- 
ing, on appropriate Riemann surfaces, into regions where the solu- 
tions typically develop singularities. The expansions have the form 
y ~ F(s;£(aO) = E^o^^m^W), where £(x) = Ce~ Xx x a - A, a 
and the functions F m are uniquely determined, modulo trivial trans- 
formations, by f; the constant C depends on the solution y. F m 
satisfy a recursive system of equations, typically simpler than the 
original system. In particular, for all first order equations and for 
Painleve's PI and P2 equations, the solution of the recursive system 
is expressible by quadratures. 

The method can be interpreted as a trans asymptotic matching 
technique in that the expansion F of y matches (and is fully deter- 
mined by) its Ecalle transseries 2 in a sector where y is asymptotic 
to a power series. The constant C in the definition of £ is one of the 
constants beyond all orders in the transseries of y. In some instances, 
the technique provides a connection method even in nonintegrable 
systems (in which case, the connection data are path-dependent). 
The constant C becomes thus accessible by classical asymptotics and 
determines the position of singularities of y. 

The expansion y ~ F(x;S t (x)) satisfies Gevrey-type inequalities, 
and thus produces exponentially accurate estimates of y (see 

Some examples are outlined. In the first one, a nonintegrable Abel 
equation, the method provides a description of the exact type of all 
but finitely many singularities of solutions in a sector, and of the as- 
sociated the Riemann surfaces. The connection between these com- 
plicated Riemann surfaces and the numerically observed chaoticity of 
solutions [2S1 is briefly discussed. 

As other examples we consider the Painleve equations Pi and Pn 
for which we use the technique to express the asymptotic distribution 
of poles near the antistokes lines of the so called truncated solutions. 



2 In our context these are algebraically determinable formal combinations of se- 
ries in x~ x and small exponentials, and generalize classical asymptotic expansions 

EE] 
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The position of the poles only depends on an exponential asymptotics 
quantity, the constant beyond all orders. 

2. Setting 

We adopt, with few exceptions that we mention, the same conditions, 
notations and terminology as ^Hj; the results on formal solutions and 
their generalized Borel summability are also taken from |13j . 
The differential system considered has the form 

y' = f(x-\y) yeC\ xgC (1) 

where 

(i) f is analytic in a neighborhood V x x V y of (0, 0), under the gener- 
icity conditions that: 

(ii) the eigenvalues Xj of the matrix A = — < f^-(0, 0) > are 

I % J i,j=l,2,...n 

linearly independent over Z (in particular Xj ^ 0) and such that 
(hi) argAj are all different. 

(In fact somewhat less restrictive conditions are used, namely 
those of HI §1.1.2.) 

By elementary changes of variables, the system can be brought 
to the normalized form ^5] , [IS] 

y' = -iy + iiy + g(x-\y) (2) 

where A = diag{Aj}, A = diag{<x,-} are constant matrices, g is ana- 
lytic at (0, 0) and g(x _1 , y) = 0(x~ 2 )+0(\y\ 2 ) as x — > oo and y — > 0. 
Performing a further transformation of the type y i— > y — a k x ~ k 
(which takes out M terms of the formal asymptotic series solutions 
of the equation) , makes 

g(|x|-\y) = 0(x- M - 1 ;|y| 2 ;|x- 2 y|) (x -+ oo; y - 0) 

where 

M > max3f?(aj) 

3 

and 0(a; b; c) means (at most) of the order of the largest among a, b, c. 

Our analysis applies to solutions y(x) such that y(x) — > as 
x — ► oo along some arbitrary direction d = {x £ C : arg(x) = eft}. A 
movable singularity of y(x) is a point x E C with G where 
y(x) is not analytic. The point at infinity is an irregular singular 
point of rank 1; it is a fixed singular point of the system since, after 
the substitution x = z _1 the r.h.s of the transformed system, ^| = 
— z~ 2 i{z,y) has, under the given assumptions, a pole at z = 0. 
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2.1. Classical versus exponential asymptotics 

In order to understand the properties of solutions of (j2j) for large 
x, one way is to find formal asymptotic solutions, then use asymp- 
toticity relations to deduce information about the true solutions from 
the formal ones. It is easy to see that equation © admits a unique 
asymptotic formal power series solution |51| 

oo 
r=2 

The coefficients yo, r of yo can be computed recursively by substitu- 
tion in (j2J) and identification of the coefficients of x~ r ; the series yo 
is a formal solution, and is usually divergent. Its Borel summability 
was shown, in a more general setting, by Braaksma [§]. 

Given an open sector of the complex x-plane, of angle less than ir, 
there exists a true solution of ^ which is asymptotic to © in that 
sector (as |x| — ► oo) [HJ- This solution is not unique in general. 

To illustrate the way different solutions with the same asymptotic 
series (in a sector) can be distinguished consider the simple linear 
equation f'(x) = —f(x) + x _1 with the general solution f(x; C) = 
e~ x Ei(x) + Ce~ x where Ei(x) = P f*^ t~ x e l dt. Any solution f(x; C) 
has the same (divergent) power series asymptotic expansion in the 

right-half plane: f(x;C) ~ fo(x) = Yl r >o r] - x ~ r l ^ or x ~~ * °°> > 
0. The parameter C which distinguishes different solutions multiplies 
the term e~ x which is much smaller than all the terms of the asymp- 
totic series /o : C is a constant beyond all orders. 

The theory of linear equations with an irregular singular point is 
well developed and there are comprehensive results; we mention the 
works of Babbit and Varadarajan 1 , Balser, Braaksma, Jurkat, Lutz, 
0, [H], |5) EH) Balser, Braaksma, Ramis and Sibuya [Hj, Deligne [T7| . 
Jurkat [SJ , Katz [HS] , Levelt [SHI , Levelt and Van den Essen j37j , Lutz 
and Schafke gU, Manin [HH], Olver 00], Malgrange EH, Ramis gl], 
Ramis and Martinet 13 , Ramis and Sibuya jH] , Sibuya |46j , Turritin 
|49j . and others — see also (HOI and the references therein. 

For linear equations there exist fundamental systems of solutions 
in terms of which one can speak of exponentially small terms. A 
formal analogue for nonlinear equations is represented by formal ex- 
ponential series. 

An n-parameter formal solution of @ (under the assumptions 
mentioned) as a combination of powers and exponentials is found in 
the form 

y(x) = Yl C k e- A - k V k s k (x) (4) 

ke(NU{0})" 
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where §k are (usually factorially divergent) formal power series: so = 
yo (see (J3J) and in general 

oo ~ 
r=0 

that can be determined by formal substitution of (J3J) in (J2J); C € C n 
is a vector of parameters 3 (we use the notations C k = YYj=i ; 
A = (Ai, A n ), a = (a 1} a n ), |k| = fcj. + ... + k n ). 

Note the structure of Q : an infinite sum of (generically) divergent 
series multiplying exponentials. They are called formal exponential 
power series [5T] . 

Formal solutions (Jlj) of differential equations ((21) were introduced 
by Fabry [2^ and studied extensively by Cope [TT] . 

From the point of view of correspondence of these formal solutions 
to actual solutions it was recognized that not all expansions (J1J should 
be considered meaningful; also they are defined relative to a sector 
(or a direction). 

Given a direction d in the complex x-plane the transseries (on 
d), introduced by Ecalle ^HJ, are, in our context, those exponen- 
tial series (Jlj) which are formally asymptotic on d, i.e. the terms 
c k e -A-kx x a,k x -r ( with k G ( N u { }) n , r € N U {0}) form a well 
ordered set with respect to 3> on d (see also |H3)- (For example, this 
is the case when the terms of the formal expansion become (much) 
smaller when k becomes larger.) 

For linear systems any exponential power series solution is also 
a transseries: it consists of n power series multiplying exponentials 
since §k = for |k| > 2. 

In the nonlinear case if a formal exponential power series (Jlj) sat- 
isfies the condition Cj = if e~ x i x -f* as x — > oo, x € d then Q is 
a transseries on d. In fact, it is clear that (jlj) is a transseries on (any 
direction of) the open sector St ra ns defined by 

Straus = {x € C ; if Cj^ then X(\ jX ) > , j = 1, n } (6) 

This sector may be empty; it may be the whole C if all Cj = 0; 
otherwise it lies between two antistokes lines, and has opening at 
most TT. 

If yo is divergent (which is generic) then the terms containing 
exponentials in (j3J (i.e. terms with |k| > 1) are much smaller than 

3 In the general case when some assumptions made here do not hold, the general 
formal solution may also involve compositions of exponentials, logs and powers 
|19|. The present paper only discusses equations in the setting explained at the 
beginning of S|3 

4 We note here a slight difference between our transseries and those of Ecalle, 
in that we are allowing complex constants. 
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all powers of x in yo and cannot be denned by classical asymptotic 
inequalities in the Poincare sense. Hence their designation: terms be- 
yond all orders. 

Transseries and their correspondence with functions are the sub- 
ject of exponential asymptotics, which developed substantially in 
the eighties with the work of M. Berry (hyperasymptotics), J. Ecalle 
(the theory of analyzable functions), and M. Kruskal (theory of tower 
representations and nice functions) (see references). 

From a historical point of view we must stress the importance 
of the fundamental work of Iwano. Generalizing earlier results of 
Malmquist he proved in wide generality that locally meromorphic 
systems of differential equations have expansions, which for the class 
(J2J) discussed in the present paper have the form 

y(s) = 0o(*) + E C k e- A - k ^ a ' k k (x) (7) 

keN" 

convergent, and with 0k analytic, in appropriate sectors [2S], |27| - 

Later, Ecalle introduced a very large space of expansions, relevant 
to differential, difference, integral and other equations. The funda- 
mental space where formal solutions are sought is purely algebraic — 
the transseries (for example, see the expansion © for solutions 
of ©)• Then a general procedure (based on Borel summation), inde- 
pendent of the equation where the expansions originate, is outlined 
to associate functions to formal expansions. As a consequence the 
summation procedure can be used in a broad class of problems and 
yields a complete isomorphism between formal expansions and a class 
of functions (analyzable functions). For differential equations this pro- 
cedure shows that all the terms 0k in (jZj) can be in fact obtained from 
00 by a form of analytic continuation (Ecalle's resurgence relations). 
Also, Stokes phenomena can be described in detail. The paper [P3*] 
proves this procedure in the context of differential systems with a 
rank 1 singularity. 

Braaksma has recently extended the theory to nonlinear difference 
equations [TO] . 

The present paper studies the solutions in a region where the ex- 
pansions ® and those of .13; diverge. It is shown that in this region 
the solutions of generic systems actually do have singularities (see 
Theorems El and |3J) grouped in regular arrays. Thus, a posteriori we 
know that expansions of the form (jJJ) cannot converge there. Never- 
theless, asymptotic representations, in terms of functions themselves 
singular derived in the present paper hold in this region (see l)26jl ) and 
they enable finding the singularities of y. The region where the se- 
ries (ffi))) is asymptotic to solutions and the region where (J7J) , or (jTU)) , 
converge do intersect so the two expansions of the same solution can 
be matched in this region (see Theorem 
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2.2. Further notations and results referred to 

We recall that the antistokes lines of (jSJ are the 2n directions of the 
x-plane i\j M+, — i\j M+, j = 1, n, i.e. the directions along which 
some exponential e~^ x of the general formal solution (jlj) is purely 
oscillatory. 

In the context of differential systems with an irregular singular 
point, asymptoticity should be (generically) discussed relative to a 
direction towards the singular point; in fact, under the present as- 
sumptions (of nondegeneracy) asymptoticity can be defined on sec- 
tors. 

A first question is to determine which are the solutions asymptotic 
to the power series solution @, and to find their regularity. 

Let d be a direction in the x-plane which is not an antistokes line. 
The solutions y(x) of ((2J) which satisfy 



are analytic for large x in a sector containing d, between two neigh- 
boring antistokes lines and have the same asymptotic series 



(see Add endix 16 . 1 1 for more precise statements and details). 

A sweeping correspondence between general transseries and the 
class of analyzable functions has been introduced in the monumental 
work of Ecalle HBJ-EOJ- 

In the context of © , a generalized Borel summation CB of trans- 
series (jlj is defined in • The rest of this section states some results 
of ^H! needed in the present paper; more details are included in Sec- 
tion ^16.21 of the Appendix. 

The formal solutions Q are determined by the equation ((2j) that 
they satisfy, except for the parameters C. Then a correspondence 
between actual and formal solutions of the equation is an association 
between solutions and constants C. This is done using a generalized 
Borel summation CB. 

The operator CB constructed in can be applied to any trans- 
series solution (@J of (J2J) (valid on its open sector Straus, assumed 
non-empty) on any direction d C Straus and yields an actual solution 
y = CBy of (J2J), analytic in a domain S an (see (|137j0 . Conversely, any 
solution y{x) satisfying Q on a direction d is represented as CBy(x), 
on d, for some unique y(x): 



y(x) ^0 (x G d; \x\ — > oo) 



(8) 



y(x) ~ yo {x G d; \x\ -> oo) 



(9) 
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y(x) = C k e- A ' k ^ M - k y k (^) 

k>0 

= C k e~ A - k:c x M - k £fiy k (x) = CBy(x) (10) 

k>0 

for some constants C G C n , where Ma = L^QyJ + 1 (L'J * s the integer 
part), and 

oo 

y k (*) = £^kfe (a' = a-M) (11) 

r=0 

(for technical reasons the Borel summation procedure is applied to 
the series 

y k (x) = x W s k (x) (12) 

rather than to s k (x) cf. ©,©). 

In any direction d, CB is a one-to-one map between the transseries 
solutions on d and actual solutions satisfying Q, see |I3J, Theorem 
3. 

The map y ^ CB(y) depends on the direction d, and (typically) 
is discontinuous at the finitely many Stokes lines, see [T3J, Theorem 
4. _ 

For linear equations only the directions AjlR+, j = l,...,n are 
Stokes lines, but for nonlinear equations there are also other Stokes 
lines, recognized first by Ecalle (the complex conjugate directions to 
Pj ;k M + cf. (|35|): see ^Hl)- CB is only discontinuous because of the jump 
discontinuity of the vector of "constants" C across Stokes directions 
(Stokes' phenomenon); between Stokes lines CB does not vary with 
d. 

The function series in (jlUj) is uniformly convergent and the func- 
tions y k are analytic on domains S an defined in (|137|) (for some 6 > 0, 
R = R(y(x),5) > — see Theorem d of flOl . 

2. 3. Heuristic discussion of transasymptotic matching 

There is a sharp distinction between linear and nonlinear systems 
with respect to the behavior beyond Straus- 

In the linear case there are only finitely many (at most n) nonzero 
y k in ([TO]) , and (flT)|) holds in a full (possibly ramified) neighborhood 
of infinity, except for jumps in the components of C, one at each 
Stokes line (see jlHlj also jlj, ^3], [SO] and the references therein). 
The map CB is continuous at the antistokes lines, and thus the trans- 
series y of y is the same on both sides of an antistokes line. What 
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changes at such a direction is the classical asymptotic expansion of 
y, because classical asymptotics only retains the dominant series in 
y, and exponentials in the transseries y exchange dominance at antis- 
tokes lines. From the point of view of exponential asymptotics, where 
transseries are considered rather than just the dominant series, be- 
havior of solutions of linear equations at antistokes lines is relatively 
simple. 

In the nonlinear case however, generically all components yk are 
nonzero and, beyond Strans (for example, in the notations of § El 
below, for arg(x) > vr/2), (|TU)) will typically blow up because of a 
growing exponential (e~ x in this example). 

The divergence of turns out to mark an actual change in the 
behavior of y(x), which usually develops singularities in this region. 
The information about the singularities is contained in H10|). 

The key to understanding the behavior of y(x) for x beyond S an 
is to look carefully at the borderline region where converges 
but barely so. Because of nonresonance, for arg(x) = ir/2 we have 
?R.(\jx) > 0,j = 2,...,ni. 5 By (JHZl) all terms in (@J) with k not a 
multiple of ei = (1,0, ...,0) are subdominant (small). Thus, for x 
near iR + we only need to look at 

yW(x) = Y J C k l z- kx x kM ^y kei {x) (13) 

The region of convergence of (|13|) (thus of (|10|l ) is then determined by 
the effective variable £ = Cie~ x x ai (since yfc ei ~ yke 1 -,o/ x k ^ ai ~ Ml ^) . 
Convergence is marginal along curves such that £ is small enough 
but, as |x| — ► oo, is nevertheless larger than all negative powers of x. 
In this case, any term of the form C k e~ kx x kMl y kei .Q is much larger 
than the terms C[e~ lx x lMl yi ei]r x~ r if k, I > and r > 0. Hence the 
leading behavior of yM is expected to be 

yW(x) ~ ^(C ie -^ Ql ) fc s fcei; o = F (0 (14) 

fc>0 

(cf. (fTTj) ): moreover, taking into account all terms in Sfc ei we get 

yW^-^-^W-E^f 1 (15) 

r=0 k=0 j=0 

Expansion (|15|) has a two-scale structure, with the scales £ and x. 

It may come as a surprise that each Fj is a convergent series in 
^ (though the whole expansion (|15j) is still divergent). 

5 In the notations explained below in 321 =0 for j > n 1 . 
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It turns out that the reshuffling 1)15)1 is meaningful and yields 
the correct asymptotic representation of yW, and therefore of y, be- 
yond the upper edge of S an . In fact, (fT5|) extends © right into the 
regions in C where y is singular, as near as (under mild assump- 
tions) 0( e - co ^-N) of these singularities. Once these two scales are 
known and once the validity of ()15)1 is proved for our class of systems 
(Theorems 1 and 3 below), it is easier to calculate the Fj by direct 
substitution of (|T5|) in ((21) and identification of the powers of x (see 
Remark and ffi).9JI . The exact form of the second scale £ is decisive 
for the domain of validity of the expansion, see H5.2I 

To leading order we have y ~ Fo (see also (fTl)) ) where Fo satisfies 
the autonomous (after a substitution £ = ) equation 

£F(,(£) = iF (£)-g(0,Fo) 

which can be solved in closed form for first order equations (n = 1) 
(the equation for Fq is separable, and for k > 1 the equations are 
linear), as well as in other interesting cases (see e.g. $5.21 ^5.3(1 . 

Assume that Fo(£) has an isolated singularity at £ = £ s . Then 
y(x) must also be singular near x s , if £(x s ) = £ s - Indeed, it is not 
difficult to see (see H4.6j) that there must exist some #(£) analytic at 
£ s so that that / g(t)F (t)dt = 1 on a small circle around £ s . Taking 
x s large we must have by {HJ <f(l + o(l))g(£(t))y(t)dt = 1 + o(l) 
on a small circle around x s . In many instances one can refine these 
arguments to see that the singularities of y(x) and Fo(£(x)) must be 
exactly of the same type. It is clear, on the other hand, that x s form 
a nearly periodic array of points as |a? s | — > oo (see Theorem |2J). 

In the following we will make rigorous these intuitive arguments 
and then proceed to explore further properties and consequences. 

3. Main results 

Let d be a direction in the x-plane (not an antistokes line). Consider 
a solution y(x) of (j2J satisfying Q hence Q. Let (fTfl|) be its rep- 
resentation as summation of a transseries y(x) (see @) on d. Let 
Straus be the sector of validity of y(x) see (jSJ). 

For simplicity we assume, what is generically the case, that no 
pj^k (see 1)35)1) lies on the antistokes lines bounding S trans- 

We assume that not all parameters Cj are zero, say C\ ^ 0. Then 
Strans is bounded by two antistokes lines and its opening is at most 

7T. 

Notations. It can be assumed without loss of generality (possibly 
after a linear transformation in x and renumbering the coordinates 
of y) that 

(a) Ai = 1, and 
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(b) Cj = for j > m (where n\ G {2, ...,n}) 

(c) arg(Ai) < arg(A 2 ) < ... < arg(A ni ) 

(d) Strans is bounded by iM+ and the direction arg(A ni x) = — ir/2 
(which are antistokes lines associated to Ai and A n J, and Strans is 
contained in the right half-plane. 

The solution y(x) is then analytic in a region S an (see ()137)) ). 

The singularities of y(x) that we find are related to the two an- 
tistokes directions bounding Strans- We will formulate the results for 
the direction iM + (and similar results hold for the other direction 
arg(A ni x) = -it/2). 

The locations of singularities of y(x) depend on the constant C\ 
(constant which may change when d crosses Stokes lines). We need 
its value in the sector between iM + and the neighboring Stokes line 
in Strans- Let d! C Strans be a direction in the first quadrant and 
consider the representation (|1U|) of y(x) on d'. 6 From here on we will 
rename d! as d. 

Fix some small, positive 8 and c. Denote 



Z = £(x) = Cie- 



(16) 



and 



£ = < x ; arg(x) G 



Also let 



7T . 7T 

-2 +5 <2 + * 



and 



»(Ajs/|x|) > c for all j with 2 < j < m} (17) 



S Sl ={xe£ ■ <*i} 



(18) 



The sector £ contains Strans > except for a thin sector at the lower 
edge of Strans (excluded by the conditions 5R(Aj-a;/|x|) > c for 2 < 
j < tii, or, if ni = 1, by the condition arg(x) > — f + 5), and may 
extend beyond iM+ since there is no condition on K(Aix) — hence 
3?(Aix) = 3?(x) may change sign in £ and 5^. 

Figure 1 is drawn for n\ = 1; £ contains the gray regions and 
extends beyond the curved boundary. 



3.1. Asymptotic behavior of y(x) in Ss 1 

Theorem 1 (i) There exists 5\ > so that for |£| < 5\ the power 
series 

oo 

F m (0=J2^ykev,m, rn = 0,1,2,... (19) 

fc=0 

6 Ci does not change at the possible secondary Stokes lines dj,k, |k| > 1 lying 
between R+ and iR+. 
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I x I =x 




arg (x) =-jt/2+8 

Fig. 1. Region T> x where l)2ri|l holds, when n\ = 1. The dark gray subregion is Ss 1 . 
Curves like the spiraling gray curves surround points in X (close to singularities 
of y) generate the region T> x . The picture is drawn with ni — 1,A = jQ,a = 
— h,5i = 3 • 10 6 ,a;o = 40. In this case Strans is a sector where | arg(a;)| < 5 — 0. 



converge (for notations see ^jj, |2J) ; \lb\) and for an estimate of 5\ 
see Proposition^. 
Furthermore 



y(x) ~ x~ m F m (t(x)) {xeS Sl , x^ oo) 



(20) 



m=0 



uniformly in Ss x , and the asymptotic representation 120\) is differen- 
tiate. 

The functions F m are uniquely defined by \2U\) . the requirement of 
analyticity at £ = 0, and Fq(0) = ei . 

(ii) The following Gevrey-like estimates hold in Ss 1 for some con- 
stants K\ t 2,B\ > 0: 



\F m (£(x))\ <K x m\B? 



(21) 



m—l 



y(x) - x~ m F m (£(x)) 

j=0 



< K 2 m\B^x~ 



[m G N + , x G S Sl ] 



(22) 
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Comments: 

1. It is interesting to remark that the constant beyond all orders 
C\ is now classically definable in terms of the expansion (|20|) because 
this expansion is unique with its the range of validity, and with the 
given analyticity properties. This is in a sense a generalization of 
Watson's Lemma in the context of transexpansions. 

2. While the classical expansion Q is valid only in any proper 
subsector of S$ 1 n {x : arg(x) < vr/2}, the representation (|20j) holds 
down to a distance going to zero as x becomes large from the (finite- 
plane) singularities of Fo, near which y(x) also develops singularities 
(see Theorem |21 and SJ5J). 

3. A similar picture holds near the lower edge of Strans- The con- 
stant Cj used in which determines the position of singularities 
(see (HHj)) °f y{ x ) related to that direction, is C m . If n\ = 1 then the 
value of C\ is the one in the fourth quadrant (which may differ from 
the one in the first quadrant due to the Stokes phenomenon on 



3.2. Singularity analysis 

We now focus on singularities of y(x) and their connection with sin- 
gularities of Fo- 

3.2.1. Definitions (cf. Figure 1) By © and (O we have F (0) = 0. 
Both Fo and y turn out to be analytic in Ss 1 (Theorems ^i) and 
Eli)); the interesting region is then £ \ Ss-^ (containing the light grey 
region in Figure 1). 

Denote by V a polydisk 

V = {{x-\y) : \x- x \ <pi,|y| < p 2 } (23) 

where g is analytic and continuous up to the boundary. 

Let S be a finite set (possibly empty) of points in the £-plane. This 
set will consist of singular points of Fo thus we assume dist(S', 0) > 8\. 

Denote by 1Z~ the Riemann surface above C \ S. More precisely, 
we assume that 1Z~ is realized as equivalence classes of simple curves 
r : [0, 1] i— > C with -T(O) = modulo homotopies in C\S. 

Let T> C TZe be open, relatively compact, and connected, with the 
following properties: 

(1) Fo(£) is analytic in an ex>-neighborhood of T> with et> > 0, 

(2) sup-p |F (OI : = P3 with p 3 < p 2 

(3) V contains : |f | < 5i}. 7 

7 Conditions (2), (3) can be typically satisfied since Fo(£) = £, + 0(£ 2 ) and 
Si < p2 (see also the examples in JU; borderline cases may be treated after 
choosing a smaller S\. 
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It is assumed that there is an upper bound on the length of the 
curves joining points in V: dx> = sup a mfr c x> ;a ,&er length(.T) < 
oo. 

We also need the x-plane counterpart of this domain. 

Let R > (large) and let X = n {x G £ : \x\ > R}. 

Let r be a curve in T>. There is a countable family of curves 7iv 
in the x-plane with £(7jv) = r. The curves are smooth for |x| large 
enough and satisfy 

7jv(t) =2iV7ri + ailn(27riiV)-lnr(t) + lnCi + o(l) (N — > oo) 

(24) 

(For a proof see Appendix Hfi.Ml ) 

To preserve smoothness, we will restrict to \x\ > R with R large 
enough, so that along (a smooth representative of) each r € T>, the 
branches of £ _1 are analytic. 

If the curve r is a smooth representative in T> we then have 
^~ 1 (r) = Ujvgn In where 7jv are smooth curves in {x : \x\ > 2R}\X. 

We define V x as the equivalence classes modulo homotopies in 
{x G £ : \x\ > R} \ X (with oo fixed point) of those curves which 
are completely contained in £ n {x : \x\ > 2R}. 

Theorem 2 (i) The functions F m (£); m > 1, are analytic in T> 
( note that by construction Fq is analytic inV) and for some positive 
B,K we have 

|F m (0| <KmlB m , £eV (25) 

(ii) For R large enough the solution y(x) is analytic in T> x and 
has the asymptotic representation 

oo 

y(x) ~ x- m F m (£(x)) (x G V x , \x\ -> oo) (26) 
m=0 

In /ac£ ; £/ie following Gevrey-like estimates hold 



rn—l 

y(x) - ^ x^F^(x)) 

3=0 



< ^ml-B^x] 



(m G N + , x e Pjj) 



(27) 



(mj Assume Fq has an isolated singularity at £ s G S and i/iai i/ie 
projection of T> on C contains a punctured neighborhood of ( or an 
annulus of inner radius r around) £ s . 
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Then, if C\ 7^ 0, y{x) is singular at a distance at most o(l) 
(r + o(l), respectively) of x n £ £ -1 ({£ s }) n "D x , as x n — > 00. 
The collection {x n } n ^ forms a nearly periodic array 

x n = 2nni + a x \n{2nm) + In C\ - In £ s + o(l) (28) 



as n — > 00. 

Some of the conclusions of the theorem hold with T> noncompact, 
under some natural restrictions, see Proposition |HJ 

Comments. 1. The singularities x n satisfy C\e~ Xn x" 1 = £ s (l + 
o(l)) (for n — * 00). Therefore, the singularity array lies slightly to 
the left of the antistokes line iM+ if 3?(ai) < (this case is depicted 
in Figure 1) and slightly to the right of iM+ if 3i(ai) > 0. 

2. In practice it is useful to normalize the system © so that ct\ is 
as small as possible (see the Comment 1 . in § 15.21 and § lfi.7|) . 

3. By (|27|) a truncation of the two-scale series (|26|) at an m de- 
pendent on x (m ~ \x\/B) is seen to produce exponential accuracy 
( e -N/- B l), see e.g. 02|. 

4. Theorem|51can also be used to determine precisely the nature of 
the singularities of y(x). In effect, for any n, the representation (|26|) 
provides o(e~ K ^ Xn ^) estimates on y down to an o(e~ K ^ Xn ^) distance of 
an actual singularity x n . In most instances this is more than sufficient 
to match to a suitable local integral equation, contractive in a tiny 
neighborhood of x n , providing rigorous control of the singularity. See 
also Sand §5\ 



3.3. Singularities for weakly nonlinear systems 

In this section we take g meromorphic in a small enough neighbor- 
hood of (0,0), but nevertheless analytic at (0,0), and only weakly 
nonlinear. Such could be the case if in a sufficiently large neighbor- 
hood of zero only one component of g is singular, and the singular 
manifold of g is approximately a hyperplane. 

Let y(x) be as in Theorem ^ denote f(y) = —A + g(0,y). By 
Theorem^(i) we have yi(x) ~ £(x) and yj(x) = 0(£ 2 ) <C yi(x), j = 
2,...,n when 1 3> |£(x)| 3> |x|~ 2 . For definiteness we assume the 
component 52 to be the only one singular in some neighborhood of 
(0,0). The precise assumptions are 
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/ 2 (y) 



h 



Aiyi + ei(y) 

_ -A 2 y2 -72y? + e 2 (y) 



/i(y) 



My) 

- A i l/j - 7j 2/i + e j (y ) (j ^ 1 , 2) 



(29) 



My) 



(y) 



fc=i 



where Qj,7j G C with ||a|| < 5(^2) 1 (see (|23j) ). ej are analytic and 
satisfy 



for some small positive e. 

Choose xo > large enough so that the function g(x _1 ,y) is 
analytic in V if p\ = Xq 1 (see (|23|l) and 



Theorem 3 For almost all values of the parameters 7j,aj with j = 
1, ...,n the following holds. Ife is small, Fo(£) andy{x) are not entire. 
Fq has isolated square root branch points on its circle of analyticity 
and correspondingly y has arrays of branch points for large x. 

More precisely, there exists £ s so that Fq is analytic for |£| < |£ s | 
and y is analytic if \£(x)\ < |£ s | — 5(x) where S(x) — > as x — ► 00. 
Furthermore, Fq is analytic on the Riemann surface of the square 
root at £ s and there is an array x n G 5^ to^ = £s + o(l) as 

|sc n | — > 00 so that y(x) is analytic on the Riemann surface generated 
by curves that encircle at most one of the x n . Near £ s and x n we have 



ej{y)\ < e for |y| < p 2 for j = 1, n + 1 



(30) 



Slip < e 



(31) 



F (0 = F A (K - e-) 




y(z) = y^((^ - x n ) 




(32) 



respectively, where yA and Fa are analytic at zero. 



4. Proofs and further results 



4-1. Further results needed 



A possible proof of Theorem ^ using only classical asymptotics con- 
cepts is sketched in ^6,101 The proof given in this section uses some 
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results in exponential asymptotics. We need the following facts proved 
in H3|. 

Denote by d the complex conjugate of d. We have 

Yu(x) = [_Y k (p)e-P x dp (33) 

Jd 

where the functions Y k have the form 

Y k (p)=p- ka '- 1 A k (p) (34) 

(Lemma 20 in [HI]), with A k analytic near zero, and along curves 
towards infinity avoiding the points pj- k defined as 

Pj;k = Aj — k • A , j = 1, m , ke Z™ 1 (35) 

Because in Strans we have 3?(Ajx) > for j < m and since k > 
it follows that (1) the j>j ;k have no accumulation point and (2) only 
finitely many of them are in Strans = {d : d C Strans} ( see |13| for 
more details). In particular, if pj ;k g" d then there exists a\ > and 
an ai-neighborhood d ai of d (i.e. d ai = {x; dist(x, d) < ai}) where all 
A k , k > are analytic. 

There exist positive constants K\, vq such that for all k > we 
have (113,, Prop. 22(ii) for W k = Y k in V' mv and T k . /3 /_i) 

sup \Y k (p)e- M \ < K lkl (36) 

Also, y k is the classical asymptotic expansion of y k ( JSj, Theorem 
3): 

y k (x) ~ y k (x) (x G de t6 , 6 G (-|, ~) ; x -+ oo) (37) 

In the following we need a better estimate of A k (see ("*"""") and 
1(3*3*1) ). than given in [*f*3*| . 

Proposition 4 For some a±, 62 > and all k and Z we have 

\A.l\p)\ < |r(-k-a / )r 1 Z!ar'e iy ° |p|+ai ^ |k| (38) 
uniformly in d ai . 

The proof amounts to minor modifications in a proof in [*""">]. We 
detail these modifications in Wi.6l 
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4-2. Proof of TheoremUl 

(a) Analyticity at 0. From PBl . ()34[). and Watson's lemma |S] we get 
as x — > 00 along d 



^0 m! ^ 



£ Ai m) (0) 



m>0 



m\x m kc *' 



s m-k.«'-l e - 



5 *=E- £rt)(0)r(m " k ' Q,) 



m>0 



m\x m kc *' 

1 



(39) 



m>0 



so that from Proposition 0] it follows that F m (£) = X^fcLo yfcei;m£ fc 
converges if |£| < 5±, where 5\ < 62, in which ball we also have 



|F m (OI < K x B?m\ 



(40) 



for some B\. 

(b) Asymptoticity. We first note that using again Watson's lemma, 
by ()38j) we have 



M ~ 

Yk-^E^ 
z — ' x l 



1=0 



<e ai Ir^k-a') -1 ! 5 2 ^ ai ' M - l {M + 



< 



v-lc-[k| k-cx'+Af+l„ -Af-1 



r(-k-a')- 62 % «1 



< 8~ M {M + l)\ ai - M - 1 
It is then convenient to write, for x € 5^, 



„k-ct'+M+l 



(41) 



y (x) =Y^C{e- kx x kai y kei (x)+ C k e- k ' Aa; x kQ y k (x) (42) 

k=0 k^0;k^fcei 
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and (cf. (PI) . flUl ) 

oo 

= ^C£e- kx x kai y kei (x) + 0(e~ c ^) := y^(x) + 0{ e - c ^) 



k=0 

Now, 



oo / M \ 

- rn 



k=0 k=0 \m=0 / 

oo / M \ 

+ Y^C k e~ kx x k ^ y kei {x) - £ yk ei ; m x~ m (43) 



k=0 \ m=0 



and theorem ^ follows from (|41|) . Differentiability simply follows from 
the fact that (1201) holds in a nontrivial sector. 



4-3. Special Gevrey estimates 

In the proofs of the main theorems we need optimal estimates of high 
order derivatives of functions of the form ip(z, Ylt=i a kZ k ) in terms 
of estimates of a^, when Ylk^o^k^ are Gevrey type series. Because 
of the truncations involved, the estimates do not follow from Gevrey 
theory. 

Let (f be an analytic function in the polydisk V = {\z\ < pi, \y\ < 
ft} C C x C and continuous up to the boundary. 
Assume first that the series 

oo 

k 



k=i 

converges and denote by a^ the truncation 

m 

a M (z) = ^ afe 



, - 

k=0 



Then 4^Lp(z; a(z))| 2= o is a polynomial in ai,...,a m , of degree 1 in 



dz 

am, • 



1 d m 1 d m 

-(p(z;a(z))\ z = = a y <^(0;a o )a m + — - — -tp{z\ ah~ 1 ^ (z))\ z=0 



ml dz m ml dz m 

(44) 
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Relation (|44jl is meaningful even when a(z) is only a formal sum 
(with no convergence conditions) — in the sense that the LHS is the 
coefficient of z m in the formal series expansion of (p(z;a(z)) at z = 
0. We are primarily interested in the Gevrey— 1 character of a(z) 
(meaning that for some c\,C2 > we have |a&| < cic^kl; see also |42p. 
Proposition [5] below is formulated in a way that permits an inductive 
proof of Gevrey type inequalities, when are defined recursively. 

Proposition 5 Assume p 2 — |ao| > 0. There exists a positive C so 
that: for any B,K > and any {a-k}k=l,...,m such that \&k\ < KB k k\, 
k = 1, 2, m we have 



1 

ml 



< K 2 B m - 1 (m - 1)1(1 + Cm- 1 log 2 m) 

z=0 

(45) 



/or some K% (see \55)) for an estimate of K 2 ). 
For the proof we need the following result. 

Lemma 6 Let a,b satisfy 1 < a < b. There exists C = C(a;b) such 
that if Bmz := Z £ (a,b) then 

^ fc=1 1 <Cm^(lnm) 2 (m € N) 



m!m~ m Z m [l - Z- 1 ]- 1 + m^Z 



(46) 



Proof of Lemma In this proof (and in the proof of Proposition [SJ) 
we write 0(f(m)) for terms that go to zero not slower than than 
f(m) uniformly in B, Z (and K). 

Let k m = \_m/Z\. For k < k m the terms k\m~ k Z k are decreasing 
in k, and increasing for > k m . Thus 

V £ fc z fe A;! < - + —r + m—^- = Bz(l + O^ 1 )) (m -» oo) 
^ ~ m m 2 m 3 ' 

k=\ 

(47) 

Denote p m = [2 In m/ In Z\ . For m large enough we have 



m > k m +p m , 1 -p m /rn > 1/2, p m > k v 



and for p < p Tl 



i=0 



'' ^ > e 2 2,n > << 
m 



p(p+l) _ Pm.(pm + 1) 



(48) 
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N2 k\(Bz) k 
Denote a^ 2 = — — - . — ; we have 
Nl ^ m\(Bz) m 
k=Ni v ' 



k=0 j=0 

while using (|4.3|) it follows that 



m—pm — 



Pm(Pm + l) 



k=0 



1 + 0(m 1 (lnm)' 
(1 - (mBz)- 1 ) 



(50) 



For k £ (k m ,m — p m ), because the terms in the sum are increasing 
we get 



a 



m-p m 



< m- 



(Pm+l)/m 
2?2 lnm/ InZ 



m~ L + 0(m~ 2 (lnm) 2 ) (51) 



Combining (jlHJ) (fBTljl. (|5T|) . and (|I7)l. Lemma El follows. 
□ 

Proof of Proposition O We keep the requirement of uniformity 
with respect to B, K in the notation O(-), as in the Proof of LemmaEl 
Let p = min{/9i,p2 — l a o|} (cf- the beginning of H4.3|) . For small s 
and y we have 



\ip(s, y) - <p(0, 0) - d s <p(0, 0)8 - dy<p{0, 0) • y| 

< 2(n+ p 1 / HI (M 2 + llyf) = ^(N 2 + l|yf) 

We choose a circle of radius r m , where 



(52) 



(m-l)!(flr m ) 



171 — 1 



Br r . 



(53) 



For large m we have Bmr m = e + 0(m _1 ), and we also see that 
the assumptions of Lemma El are satisfied. In particular we have for 
Isl = r m that 



m— 1 
fc=l 



< KBr m {\ + (1 - l/e) -1 )(l + C^m^lnm) 2 )) (54) 



Noting that £ s m+1 ^ m ~ 1 \s)ds = we have, using (JHSJ) and <|52|). 
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1 d" 



ml dz 7 





— <£ 


z=0 


ZTTi J r 



ds 



^(s;aN( s )) 



Q m+1 



< Ml + ^ 2 ^ 2 )|r m | 2 - m (l + O^^lnm) 2 )) 

+ K 2 B 2 )B- l {B m ~ l {m - 1)!)(1 + C^m^lnm) 2 )) (55) 



□ 



Remark 7 ^4 direct calculation shows that the expansion in \2b}) is 
a formal solution of for large x iff the functions F m are solutions 
of the system of equations 



_d 
d_ 



F = r 1 (iF -g(0,F )) 
F m + jVF„ 



(56) 



d 



ai— F m _i + Rm-l form>\ (57) 
d? 



where N is the matrix 



(T\d y %(Q,¥ Q )-A) (58) 
and i/ie function R m _i(£) depends only on the F& A; < m: 



£R 



m— 1 



(m -1)7 + 1 



m— 1 



1 d m 

ml dz m 



m—l 
3=0 



2=0 

(59) 



(see also \44\))- 



4-4- Proof of Theorem^ (i) 

By Theorem^ F m , m > are analytic for |£| < Si. Furthermore, 
F m are analytic in the eu-neighborhood of T> since by assumption, 
Fo is analytic there and equations (|57|) are linear for rn > 1. 

We set as initial conditions for (fST|) in V the values of F m (£o) 
provided by Theorem^at a point (o£P with |£o| G 

For £ E 2? and m > 1, (fST|) can be integrated, yielding the recursive 
system 
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F m (0 = M(£;£ )F m (£ )+a 1 [F m ^(0 - M(£;£ )F m -i(fo)J 
-a x t M(Z;s)N(s)F m _x(s)ds+ f M^;s)It rn - 1 (s)ds (60) 

ho ho 

where M(£; C) is the fundamental solution of 



^- + NM = with M(£;Q\e=c = I (61) 
Direct estimates in (|6ll)l and (|ST?|) using (flU)) give 



||F m ||2? < MK x m\Bf + (|ai| + M)\\ 

+ 2|a 1 |d I ,M5 2 - 1 (||g|| + ||yi||)|| F m -i \\t> 

+ 25~ 1 M 



(m + A) ||-P m -i|| +d v 



1 H m 



to! dz r ' 



(62) 



where 



M= sup ||M(e;C)ll; ||F||x» = sup|F(0|; l|g||= sup |g(z,y)| 
Uev £sx> | 2 |<pi,|y|<p 2 

A = max ||i^(0|| < ih/2)- 1 {\\g\\ + P||) (63) 



Choosing K, B large enough, the proof of (|25|) is immediate induction 
from (|62|) and Proposition [5J 



4-5. Proof of theorem^ (ii) 

We will prove (|26l) at each point x = x a E (with uniform estimates 
on Pa;). x a is the endpoint of a curve 7at in I? x with £;(7at) = r curve 
in V and satisfying (|24|). 

Denote a = £(x a ). If |a| < <5i then ()26jl follows from Theorem ^ 
so we assume \a\ > 5\. Then we can choose r to go from along a 
direction up to the circle |£| = 5\, not re-entering the circle. 

Let to such that £ = r(t ) G (|5i,5i) and denote -T = r|[ to l ], 
7° = 7Jv|[to,i]i then 7 lies in a bounded region. 

We prove Q26JI in a small, connected, simply connected neighbor- 
hood A/lo of 7 . 
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Denote 6{x) = y(x) - F^(x). 

To estimate S{x) we use a contraction argument if m is not too 
large {Case I) and a direct argument for m large {Case II). 

Let Co be positive and small and let m a be the maximal integer 
such that 
m\{2B/\a\) m < c . 

Case I: m < m a 

First the differential equation satisfied by F^{x) will be written, 
then the equation of 6{x), and finally the independent variable will 
be changed to £, yielding 1)65(1 . 

With the notations introduced in M,3I and denoting by F{x) the 
formal series X^fcLo x_fc F m (£(x)) we get from Remark[7] 



■fpW( x ) = (-A + -A)F^{x)+g{x- 1 ,F^{x))^ 
ax \ x ) 

+ JtT Af ^(^ + «(*)) - ^F m fc(x)) (64) 

Using (|64j) and ©, a direct calculation yields the equation for 
d{x). The map £(x) is a biholomorphism of A/Lo onto a neighborhood 
A/po of -T ; changing the independent variable from x to £ we get 



k +m = T °{^) s+T iwr s )M^) (65) 

where 



\x J 

\x , 



1 



1 



x£ a/x — 1 
1 1 
£ a/x — 1 
1 1 
£ a/x — 1 



-aTl + A + a9 y Sf(0,F 



x x 



3(-,F 

x 



m] \ [m] 



-a y 0(o,F o )5 

s 7^ 17 m 



r m+l 



+ - 



.m+1 



(ml + A)F r 



where To, 1,2 are clearly well defined for small enough £ and (5. Fur- 
thermore, they are well defined for £ G A/po if |<5(£)| < (p2 — P3)/2 
and for large enough (see Appendix H6.5|) . 

As in (|6U|) we obtain for 6 the integral equation 
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S = J (6) where J = Jo + Ji 
with Jo(0 = M(£; £ )<5(£o)+ f Mfc s)T 2 (-i^ ds 



+ J M(^,s)T (^-^jS(s)ds (66) 

Let B be the Banach space of C™-valued analytic functions 5 on 
A/po, continuous up to the boundary, and satisfying <5(£o) = 0, with 
the norm \\S\\ = sup^^ |<5(£)|. 

The integral operator J of ()66() is defined on the ball of radius 
(P2 — P3)/2 in B. We will show that it invariates a ball in B and that 
it is a contraction there. As a consequence, the integral equation (|66|) 
has a solution (5 which is analytic on A/po, therefore y(x) is analytic 
on A/y) ; we will also obtain estimates for which will prove (|26[) in 
Case /. 

We will denote by const a constant independent of a, N, m, B, en, R. 
It will be assumed that B,R > 1, cq < 1. 

Note first that the assumption of Case I implies 

m/|o| < const (67) 

To estimate J{8) note first that 

llm „ const . . 

T < 68 
\a\ 

By (EU when a is large we have |5(£ )| < K 2 {m + lJLBj" + 1 a -m-1 . 
By Theorem [2] (i) , since £ varies in a compact set independent of a, 
and then \M\ < M as in (|63j) . Also estimating derivatives with the 
Cauchy formulas on circles \x— x'\ < pi/2; |y — y'| < ejy/2 and taking 
a large so that |o| < 2|a — du| we get 



|Ti|| < " y — j-r + < const — + ||<5|| ||<5|| (69) 



and 



|9«Ti||< rr ^— + ^P- (70) 
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Also 



M(£;s)T 2 

& 

and using (|BT)l 

so that, 
IIW)II < 



< 



2MKd v (2B) m m\ (2\a 1 \d T > m 



\<i\c p lal 



< consi 



(2B) m m\ 



(71) 
(72) 



M<5(£ )+ / M(£;s)T 2 

& 



ris 



< const 



(2 J B) m m! 



Prom (jggj). PH> we get 



||Ji(5)|| < const ||5|| ^ + ||<5 



1 



< const cq (73) 



(74) 



Also, from (|HH|) . (fHT^ . (fTHj) 



\\J{8 l )-J{8 2 )\\ < COnst -r-r + \\8\ 

a 



(75) 



It is easy to see that for positive constants R, Kq large enough, and 
Co small enough the following holds: if \\S\\ < KqCq then : \\S\\ < P2 ~ P3 ; 
11^(5)11 <# co, also 

(76) 



and ||jr(«5i) - J(S 2 )\\ < X\\Si - S)\\ with A < 1. This shows the 
existence and analyticity of d. 

Finally, to obtain the needed estimate (|26j) note that using 1)73(1 . 
(|75|). we get 

||<5|| = ||J(5)|| < HJbll + \\Ji(S)\\ < const^^ml + h\6\\ 

\a\ m 4 

so that using (|67|) and Lemma QUI 

(2B) m , (2,B) m + 1 / 
o < const — — m! < const — — ; — — — (m + 1)! 

ii H | Q |m | a | m+1 

< const — ; — (m + 1)! 



m+l 



which concludes the proof in Case I. 
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Case II: m > m a 
In this case 



Since 



1511 < 



yO) 



m\(2B/\a\) m > c 
1 



(77) 



Fo Fi 

x 



k=2 



using the result of Case I to estimate the first term 



< const 



(5B) 2 2\ 



m max 



K2\B 2 Km\B r 



and since (2B/\a\) 2 2\ < c < (2B/\a\ 



x 



i m m! 



8B 2 2\ T ^ , (2B) m . ,, (2B^ r " +l 

const— — — +mivm!— < const (m+lj! ( - — - 

|a| m \ |a| 



25 
a 



Fro m (|77|l |a|/(2B) < consi m Cq < const Cq" 1 ; using this, (|77|) 
and (1140(1 we finally get 



<K(c )(m + l)! (4B/co) 



m+l l^l-™- 1 



where -ftT(co) is a constant dependent of Co- □ 

We should stress that while the estimates in this proof clearly 
show the Gevrey character of the expansion, they are very far from 
optimal. In fact the substantial increase in B in the arguments was 
artificially introduced to make the calculations less cumbersome. 

The following is an extension, in some respects, of Theorem |U(ii). 



Proposition 8 Assume T> is not necessarily compact, r is a curve 
of possibly infinite length in T> with the following properties: 

(a) For some e > 0, Ti^^, S) and N(z) are analytic for z in an 
e neighborhood of T and for \S\ < e and in addition Ti^z, 5) = 
0{zd,S 2 ) 

(b) M(£>£io) is bounded in an e neighborhood of T and for some K 

and all £ € r we have Jj5 q M(£,£i,o) d|s| < K (where \M\ is some 

Euclidian norm of the matrix M(£,£io),)- 

Then the conclusions of Theorem^ (ii) hold in the x domain V x 
corresponding to T> . 



Noting that 



M(£,a,o; 



d|s| is a finite measure along r, the proof is 



virtually identical to the proof of Theorem |21 
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4-6. Proof of Theorem^ (Hi) 

We need the following result which is in some sense a converse of 
Morera's theorem. 

Lemma 9 Let B r = {£ : |£| < r} and assume that /(£) is analytic 
on the universal covering of B r \{0}. Assume further that for any 
circle around zero C C 5 r \{0} and any g(£) analytic in B r we have 
§c fiOdiO^ = 0- Then f is in fact analytic in B r . 

Proof. Let a G -B r \{0}. It follows that fr f(s)ds is single-valued in 
i? r \{0}. Thus / is single- valued and, by Morera's theorem, analytic 
in £> r \{0}. Since by assumption <f c f(£)£ n d£ = for all n > 0, there 
are no negative powers of £ in the Laurent series of /(£) about zero: 
/ extends as an analytic function at zero. 

To show Theorem |21 (hi), assume £ s is an isolated singularity of 
Fo (thus £ s 7^ 0) and X = {x : £(2) = £ s }. By lemma 03 there is 
a circle C around £ s and a function analytic in B r (£ — £ s ) such 
that j> c Fo(£)g(£)(i£ = 1- In a neighborhood of x n G I the function 
/(x) = e~ x x ai is a biholomorphism and for large x n 



It follows from lemma^lthat for large enough x n y(x) is not analytic 
inside C either. Since the radius of C can be taken o(l) Theorem |2] 
(hi) follows. 

Note. In many cases the singularity of y is of the same type as 
the singularity of Fq. See ^Slfor further comments. 



4-7. Proof of Theorem^ 

As in £14.51 we can reduce to the study of © in Mpo , where the func- 
tion £(x) is biholomorphic and we can change variables to £. In this 
variable both and (|56[) assume the form (where x = and F 
is F or y) 




(78) 
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at, 

dF 2 _ \ 2 F 2 - l2 F? + e [ 2 l] (x-\F) 



^ = XjFj " ljF ' + e ^ {x ' 1,F) {j + h 2) 



(79) 



where 



hte,F) = l-^a k F k + ( & 1 (x- 1 ,F), 
i=i 

with analytic in V and ||e^||-p < e, % = j, ...,n + 1 (for Fo, we 

havee' 1] (x~ 1 ,F ) = e j (F )). 

Generically a 2 is nonzero. Then the analytic change of variables 
to Fx, h, F3, F n leads to a system of the form 



+ 



JV2 j>3 



-A 2 + A 2 ajFj - a 2 j 2 F } 



XA 



^ V-) -,./•?•<; (i> 2) 



(80) 



(81) 



(82) 



The substitution F x = f + = bj£ 2 + ^ (J > 2), with 

(Xj - 2)~ 1 7 i , in ftt ty -jB5 )l yields 



d£ 



d/i 



£/* ^7 = ^ + ai(A 2 -/*)(£ + /1) + £ - ^jh)(b m f + fj) 



i>3 



+(/i + 2 E fl A-^-f 

i=3 



[3] 
2 



= A,/, + 2f7i/i + 7,/i 2 - f (j > 2) 



(83) 
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According to the hypothesis of Theorem 01 it is useful to analyze first 
the equation (describing the leading order behavior of h) 



hti = {hC 1 + di + d 2 0h + {-\2C 1 + d 3 + d A Z); h(0) = 1 (84) 

(this Abel type equation cannot be solved in closed form, in general). 
In integral form, 



H0 2 = i + dfk 2 + 4k 

+ / (df ] s + df)h(s)ds + 2\ 2 (h{s)-l)s~ 1 ds (85) 
Jo Jo 



Lemma 10 (i) Equation {8$ has a unique solution ho analytic at 
£ = 0, with ho(0) = 1. 

(ii) For a generic set of d\, ...,d^ the solution ho is not entire and, 
on the boundary of the disk of analyticity, hg has square root branch 
points. 

Proof, (i) It is straightforward to check that, since A2 N (see 
then H85|) has a (unique) formal solution of the form h = 1 + 
EfcLi ~ h kt, k (where h x = (df 1 + df)(2- dfy 1 ). To show h converges 
we take h = 1 + Y,k=l ~ h kt + £ M ^m(0 in ©: 



2i M h M {i) 

= Qm 2M h M (o 2 +e i R(o-4 ] / (t M + x>* M+fe )M*)^ 

J ° k=0 

(86) 

with analytic, or 

2h M (o = Q(m M h M (o 2 +R(o- 

df / (s M ~ 1 + Y j b k s M+k e +1 )h M {is)ds = 2J{h M ) 
Jo k=0 

J is manifestly contractive in the sup norm for small £, if M > 
(ii) The proof, elementary but delicate, is given in ^16.81 □ 



4 
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Lemma 11 Let £o be a branch point singularity on the boundary of 
the disk of analyticity of ho (see Lemma fflJl (ii))- Assume in \8ty 
is small enough and analytic in a (large enough) neighborhood m£, F 

o/(£o,F (&))- Then 

(i) For some < 5\ < 62, F(£) and h(£) are analytic in the cut 
annulus {£ : |£ - £ | G (Si,S 2 ), arg(£ - £ ) / 0}. 

(ii) h and F have a square root branch point at some £ s with 
C s -Co = 0(\\e^\\). 

Proof. We substitute h = ho + fi cf. Lemma ITU1 f satisfies a system 
of the form 

£f = JV(£)f + e(£,f 2 ) 

or 

f = fi + M(0 J M' 1 (s)e{s,f 2 (s))ds 

where the matrices N \M , M _1 and e are analytic in Mp°- Part (i) 
follows now in the same way as Theorem |2] (ii), and \\y — Fo\\j\f r0 = 

0(||e^||). (ii) In a small neighborhood of £0 by part (i) and Proposi- 
tion (generically) j| h 7^ and we may change variables in (|83|) so 
that h is the independent variable (and £ = £(/i). We note that 

d£ £h 

dh~ A(f,Oh + B(f,0 + e(f,h,0 

and /i, e are small while generically S(/(£o)>£o) is n °t small. Then 
(|83|) in the variable h, with initial condition £(/tQ + 0(||e[ 3 l ||)) = £0 + 
0(||e[ 3 ]||), has a solution which is analytic near h = 0. Furthermore 
it is easy to see that (under the same genericity assumptions) we 
have d h Ch=o = but det5 feft (/i,£, .-,f n )h=o + and then Fj(g) = 
Fj((€ ~ £o) 1/2 ) with Fj locally analytic. 

5. Examples 

5.1. Example 1 

We first illustrate how singularities of solutions are found (using 
transasymptotic matching) on a first order Abel equation 8 : 



u = u 3 - z (87) 
The authors are grateful to A. Fokas for pointing out to this example. 



32 



O. Costin, R. D. Costin 



the first example on which nonintegrability was shown using Kruskal's 
poly-Painleve analysis |34j . 

The study of (|87|) is done in the following steps. Classical asymp- 
totics of differential equations |5J shows (and it also follows from the 
analysis below) that for z — > oo with argz G (^j7r, ^jvr) there is a 
one parameter family of solutions u = u(z; C) such that u(z; C) = 
z 1 / 3 (l + o(l)). Then u ~ u = z 1 / 3 Y^=o z 5fc% ■ The parameter C may 
be chosen to be the constant beyond all orders, see ^5.1.21 

After proper normalization of (|87|) (see H5.1.1|) Theorems ^ an d El 
are applicable and provide a global asymptotic description of u(z; C) 
in a region where the solution is analytic and surrounds its singulari- 
ties for large z (Proposition I13|) . These are algebraic branch points of 
order —1/2 (see 1)102(1 ) and their location, dependent on C, is deter- 
mined asymptotically. Conversely, the C of a particular solution can 
be determined from the asymptotic location of one singularity. 

5.1.1. Normalization Formal solutions provide a good guide in find- 
ing the normalization transformations. A transformation bringing the 
equation to its normal form also brings its transseries solutions to the 
form (|10|) . It is simpler to look for substitutions with this latter prop- 
erty, and then the first step is to find the transseries solutions of 

(EH). 

Power series solutions. Since at this stage we are merely looking 
for useful transformation hints, rigor is naturally not required. Sub- 
stituting of u ~ Az p in (|87|) and looking for maximal balance [S] give 
p = 1/3, A 3 = 1. Then u ~ Az 1 / 3 + Bz q with q < 1/3 determines 
B = I A 2 , q = -4/3. Inductively, one obtains a power series formal 

solution u = Az l / 3 {1 + £Xl uo,kZ~ 5k/3 ). 

General transseries solutions of \8'l\j . In order to determine the 
form of the exponentials in the transseries of u, the method is to look 
for transcendentally small corrections beyond no, by linear perturba- 
tion theory. Substituting u = uq + 5 in (|H7j) yields to leading order in 
S, the equation 



5' = (3A 2 z 2 / 3 + ^\ 5 (88) 

whence 5 oc z 2 ! 3 exp (|^4 2 z 5//3 ) . In Q the exponentials have linear ex- 
ponent, with negative real part. The independent variable should thus 
be x = -(9/5)A 2 z 5 / 3 and 9fc(x) > 0. Then u = x 1 / 5 XXo u 0;kX~ k , 
which compared to © suggests the change of dependent variable 
u(z) = Kx l ^h{x). Choosing for convenience K = A 3 ^{— 135) 1 / 5 
yields 
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ti + — /i + 3/i 3 - - = (89) 
ox 9 

The next step is to achieve leading behavior 0(x~ 2 ). This is easily 
done by subtracting out the leading behavior of h (which can be 
found by maximal balance, as above). With h = y + 1/3 — x _1 /15 we 
get the normal form 



y' = -y + + g{x-\y) (90) 
ox 



where 



rfx-',„) = -Stf +»■) + g - -L, _ _|_ + _L_ (91) 

We see that 

A = l, a = 1/5, and thus £ = Cx 1/B e~ x (92) 

5.1.2. Definition of C for a given solution u(z) After normalization 
l|9~0~|) the results in [2] apply, and the constant C is uniquely associ- 
ated to a u(z) on a direction arg(z) = </> as the limit 

c = & ec*)" 1 («(*)- £ ^Sw) ( 93 ) 

arg(z)=</> y fc<|a:(z)| / 

This limit exists for all € (^j 71 ") ^j 71 ") an d is piecewise constant, 
with one jump discontinuity at the midpoint of this interval. The 
value of C relevant to the singularities of u is the one nearest to the 
edge of Straus where these singularities are calculated, as follows from 
Theorem ^ 

5.1.3. Finding the two-scale expansion \2b]) Having the second scale 
given by ()92|) and all the conditions of Theorem ^ satisfied, the sim- 
plest way to calculate the functions in y = Y^k=o x ~ k Fk{Q is by 
substituting y = y in (JUUJ) and solving the differential equations, as 
in the proof of Theorem |21 (i) ; the equation for Fn(£) is, cf. (Jn7)|) . 

= F (l + 3F + 3F 2 ); i^(0) = 1 (94) 

and, cf. (fo7|). 
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iF' k = (3F + l) 2 F k + R k (F , F k _ x ) 

( for k > 1 and where i?i = ^F 3 ) (95) 

The first term F$ of the expansion of u is then given by 

t = ZoF (0(Fo(0+u>or e (F (0+uJt)-~ e (96) 

with Co = 3- 1 /2 eX p(-i^^/3), too = \ + and = 5 + The 
functions Fk, k > 1 can also be obtained in closed form, order by 
order. 

By Theorem^ the relation y ~ y holds in the sector 

S Sl = {xeC: arg(x) > + 5, \Cx 1/5 e~ x \ < 

for some <5i > and any small 5 > 0. 

Theorem 3 insures that y ~ y holds in fact on a larger region, 
surrounding singularities of Fq (and thus of y) . To apply this result we 
need the surface of analyticity of Fq and an estimate for the location 
of its singularities. 

Lemma 12 (i) The function Fq is analytic on the universal covering 
TZs of C \ S where 

Z = {t P = ("1)^0 exp(p 2 vrV3) : p 1)2 £ Z} (97) 

and its singularities are algebraic of order —1/2, located at points 
lying above S. 

(ii) (The first Riemann sheet) The function Fq is analytic in C\ 
((-oo,£o]U[£i,oo)). 

(Hi) The Riemann surface associated to Fq is represented in Fig. 

2. 

Proof 

Singularities of Fq. The RHS of Q94|) is analytic except at Fq = 
oo, thus Fq is analytic except at points where Fq — > oo. From (|96|) 
it follows that lirnF _ i . oo £ G S and (i) follows straightforwardly; in 

particular, as £ — > £ p G S we have (£ — C P ) 1 ^ 2 i ? o('?) — ► \J — £p/6. 

(ii) We now examine on which sheets in 72.= these singularities 
are located, and start with a study of the first Riemann sheet (where 
^b(C) = £ + 0(£ 2 ) for small £). Finding which of the points £ p are 
singularities of Fq on the first sheet can be rephrased in the following 
way. On which constant phase (equivalently, steepest ascent /descent) 



Formation of singularities in nonlinear differential systems 



35 



paths of £(Fo), which extend to \Fq\ = oo in the plane Fq, is £(Fq) 
uniformly bounded? 

Constant phase paths are governed by the equation ^(dln^) = 0. 
Thus, denoting F = X + iY, since f'/f = (F Q + 3F 2 + 3F 3 ) -1 one 
is led to the real differential equation 9f(£'/£)dX + 3t(£'/£)dY = 0, or 

Y{1 + 6X + 9X 2 - 3Y 2 )dX 

- {X + 3X 2 - 3Y 2 + 3X 3 - 9XY 2 )dY = (98) 

We are interested in the field lines of Q98|) which extend to infinity. 
Noting that the singularities of the field are (0, 0) (unstable node, in 
a natural parameterization) and P± = (— 1/2, ±\/3/6) (stable foci, 
corresponding to — ZJq and — ujq), the phase portrait is easy to draw 
(see Fig. 2) and there are only two curves starting at (0, 0) so that 
\Fq\ — ► oo, £ bounded, namely ±M + , along which £ — ► £o and £ — ► £i, 
respectively. 

(hi) Thus Fig. 2 encodes the structure of singularities of Fq on 
1Z- in the following way. A given class 7 G 7£,~ can be represented 
by a curve composed of rays and arcs of circle. In Fig. 2, in the Fq- 
plane, this corresponds to a curve 7' composed of constant phase 
(dark gray) lines or constant modulus (light gray ) lines. Curves in 
IZs terminating at singularities of Fq correspond in Fig 2. to curves 
so that \Fq\ — > 00 (the four dark gray separatrices Si, S4). Thus to 
calculate where, on a particular Riemann sheet of IZs, is Fq singular, 
one needs to find the limit of £ in (|96|) . as -Fo — » 00 along along 
7' followed by S 1 ,. This is straightforward, since the branch of the 
complex powers 6, 6, is calculated easily from the index of 7' with 
respect to P±. □ 

Theorem |21 can now be applied on relatively compact sub domains 
of IZs an d used to determine a uniform asymptotic representation 
y ~ y in domains surrounding singularities of y(x), and to obtain 
their asymptotic location. Going back to the original variables, similar 
information on u{z) follows. For example, using Theorem [21 for the 
first Riemann sheet (cf. Lemma IT2*1 (ii)) 

V = {\i\ <K I ^(-00,6) U (£ , +00) , |e — eo| >e,|£-6l >e,} 

(for any small e > and large positive K) the corresponding domain 
in the z-plane is shown in Fig. 3. 

In general, we fix e > small, and some K > and define Ak = 
{z : avgz G (^vr - 0, ^vr + 0) , \$(z)\ < K} and let K K ,s be the 
universal covering of S H Ak and H z -K e the corresponding Riemann 
surface in the z plane, with e- neighborhoods of the points projecting 
on z(x(E)) deleted. 
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Proposition 13 (i) The solutions u = u(z;C) described in the be- 
ginning of 33 have the asymptotic expansion 



u{z) ~ z 1 ^ 



+ + f 

V y k=0 



z 5k/3 I 

(as z — ► oo; # £ TZ z -K,e) (99) 



where 



Mz) = x(z) 1 / 5 e- a '( 2 \ and x(z) = --z 5 / 3 (100) 

5 

(mJ In the "steep ascent" strips arg(£) G (01,02)1 |«2 — oi| < 
7r starting in Ak and crossing the boundary of Ak, the function 
u has at most one singularity, when £(z) = £0 or £i> an d u(z) = 
z i/3 e ±27rr/3^ _|_ (i)) as z _ > qq (t ne s ig n i s determined by arg (£),). 

(Hi) The singularities of u(z;C), for C / 0, are located within 
0(e) of the punctures oflZ z -K,o- 

Applying Theorem |21 to (|90j) it follows that for n — > 00, a given 
solution y is singular at points x P;n such that £(x P)n )/£ p = 1 + o(l) 
large). 

Now, y can only be singular if \y\ — ► 00 (otherwise the r.h.s. of (|90j) 
is analytic). If x p>n is a point where y is unbounded, with 5 = x — x Pin 
and v = l/y we have 



— = vF s (w,o) (101) 
di> 

where F s is analytic near (0,0). It is easy to see that this differential 
equation has a unique solution with 5(0) = and that S'(0) = as 
well. 

The result is then that the singularities of u are also algebraic of 
order —1/2. 

Proposition 14 If zq is a singularity ofu(z;C) then in a neighbor- 
hood of zq we have 



u = ±y/=T/2(z - zq)- x ' 2 M{* ~ z o) 1/2 ) 
where Aq is analytic at zero and Aq(0) = 1. 



(102) 



Formation of singularities in nonlinear differential systems 



37 




Notes. 1. The local behavior near a singularity could have been 
guessed by local Painleve analysis and the method of dominant bal- 
ance, with the standard ansatz near a singularity, u ~ Const. (z—zo) p . 
Our results however are global: Proposition 1131 gives the behavior of 
a fixed solution at infinitely many singularities, and gives the posi- 
tion of these singularities as soon as C\ (or the position of only one 
of these singularities) is known (and in addition show that the power 
behavior ansatz is correct in this case). 

2. Eq. (|90|) can be brought to a form similar to that in TheoremOl 
by the substitution y = v/{\ + v) in ()90|) . The result has the form 



v' = -v - 27 — 10 v 2 + — v + g [1] (t-\v) (103) 

where is a now an 0(t~ 2 ,v~ 2 ) polynomial of total degree 5. The 
singularities of v are at the points where v(t) = —1, and are square 
root branch points, as in Theorem|31 whose technique of proof would 
have also applied, if the more explicit formula (|96|) was unavailable. 



5.2. Example 2: The Painleve equation P\. 

The Painleve functions were studied asymptotically in terms of dou- 
bly periodic functions by Boutroux (see, for example, [21] )• Solutions 
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Fig. 3. Singularities on the boundary of Strans for 1871 . The gray region lies in 
the projection on C of the Riemann surface where 199H holds. The short dotted 
line is a generic cut delimiting a first Riemann sheet. 

of the Pi equation turn out to have arrays of poles and they can 
be asymptotically represented by elliptic functions whose parameters 
change with the direction in the complex plane. Joshi and Kruskal 
carried out this type of expansions for generic solutions, which have 
poles throughout a neighborhood of infinity, to sufficiently many or- 
ders to determine how the parameters of the elliptic functions vary 
|29| . |3Uj . and applied this method to solve the connection problem. 
However, there exist special, one-parameter, families of solutions of Pi 
(the truncated solutions, important in applications) that are free of 
poles in some sectors. These solutions have the same classical asymp- 
totic expansion in the pole free sector to all orders and cannot be 
distinguished by classical asymptotics there. They differ by a con- 
stant C beyond all orders which can be determined by exponential 
asymptotic methods. The results of the present paper apply to this 
special family of solutions and give an asymptotic representation uni- 
formly valid at the ultimate array of poles (neighboring the pole free 
sector) and make the link between the position of these poles position 
and the value of C. 

We note that the behavior of the triply truncated 9 solutions, which 
in a sector have C = and, consequently are pole-free in larger sec- 
tors, does not follow immediately from our analysis. But this case can 
be treated by a similar methodology since after continuation across 

9 They are also known as " doublement tronquees" . 
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a Stokes line the value of C becomes equal to a Stokes multiplier, 
generically nonzero. 

This example extends the asymptotic expansions of |16j to larger 
regions of the complex plane, and also to all orders. 

We consider solutions of the Painleve Pi equation (in the form of 
|25| . which by rescaling gives the form in 24) 



^ = % 2 + * (!0 4 ) 

in a region centered on a Stokes line, say d = {z : argz = ir}. 

To bring 1)104(1 to a normal form the transformations are suggested 
by the general methodology explained in ^5.1.11 There is a one pa- 



rameter family of solutions for each of the behaviors y ~ ±< for 



large z along d. We will study the family with y ~ + w , since the 

other can be treated similarly. Its transseries can be obtained as in 
the previous example, namely determining first the asymptotic series 
yo, then by linear perturbation theory around it one finds the form 
of the small exponential, and notices the exponential is determined 
up to one multiplicative parameter. We get the transseries solution 

— oo 

/Ee*y* (105) 

k=0 

where 

£ = £( z ) = Cx- 1/2 e~ x ; with x = x(z) = ^^/^ (106) 
and yt are power series, in particular 




- , 1 '_]__ 1/0;^ 

V0 8V6(-Z)V2 28-3 z5 - (-*)5*/2 - 

We note that in the sector | arg(z) — ir\ < the constant C of 
a particular solution y (see (|109[) ) changes only once, on the Stokes 
line arg(z) = 7r |13j . 

As in Example 1, the form of the transseries solution (|105[l . (|106|> 
suggests the transformation 

(~24z) 5 / 4 /=J 

x = : viz) = \ Y(x) 

30 y ' V 6 V ' 

which, in fact, coincides with Boutroux's (cf. 24 j); Pi becomes 

Y"(x) - \ y 2 (x) + ~ = -- Y'(z) + ^ \ Y(x) (107) 
2 2 x 25 x z 
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For the present techniques to apply equation (|1()7|) needs to be 
fully normalized and to this end we subtract the O(l) and 0(x _1 ) 
terms of the asymptotic behavior of Y(x) for large x. It is conve- 
nient to subtract also the 0(x -2 ) term (since the resulting equation 
becomes simpler). Then the substitution 



transforms Pi to 



Y{X) = 1 " 25^ + HX) 



1., . 1.2 392 



y + _*_ fc _ 5fc >__ = o (io 8 ) 

Written as a system, with y = (h, hi) this equation satisfies the as- 
sumptions in ^21 with Ai 5 2 = dbl, a% t 2 = —1/2, and then = 
Ce~ x x~ 1 / 2 . The results of the present paper apply to the normal 
form (|108[) of Pi and we will prove Proposition 1151 below which shows 
in (i) how the constant C beyond all orders is associated to a trun- 
cated solution y(z) of Pi for arg(z) = ir (formula (|109j0 and gives 
the position of one array of poles z n of the solution associated to C 
(formula (jllUjO . and in (ii) provides uniform asymptotic expansion 
to all orders of this solution in a sector centered on arg(z) = ir and 
one array of poles (except for small neighborhoods of these poles) in 
formula (fTT2|) . 

Proposition 15 (i) Let y be a solution of j!U4\) such that y(z) ~ 
for large z with arg(z) = ir. For any £ (ir,ir + |tt) the 
following limit determines the constant C (which does not depend on 
4> in this range) in the transseries y of y: 




i^^'IV^'- £ l^ l=c (109) 

arg(z)=</> \ k<\x[z)\ 

(Note that the constants yo-k do not depend on C). With this def- 
inition, if C 7^ 0, the function y has poles near the antistokes line 
arg(z) = 7r + g7r at all points z n , where, for large n 



24 V n V 8 4vr 600vr 2 j 



n s 



with L n = jrjj In [jj^-Tij, or, more compactly, 
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e(^) = 12+ ( _ 24 ^ )5/4 +0(z- 5/i ) (z n ^oo) (111) 

(ii) Let e £ M + and define 

Z = {z: arg(z) > ^; |£(z)| < 1/e; - 12| > e} 

5 

fi/ie region starts at the antistokes line arg(z) = g7r and extends 

slightly beyond the next antistokes line, arg(z) = If y ~ y^—z/6 
as \z\ — ► oo, arg(z) = 7r, i/ien /or z£ Z ioe /iaue 



A 1 f 30*g t (Q \ 

^ 8V6(-z)5/2 ^ (-24z) 5k /*J 

(\z\ oo, z£2) (112) 

T/ie functions are rational, and Ho(£) = ^(^/12 — 1)~ 2 . T/ie ex- 
pansion holds uniformly in the sector arg(z) £ (3/5,7/5) 
and a/so on one of its sides, where Hq becomes dominant, down to 
an o(l) distance of the actual poles of y if z is large. 

Proof. We prove the corresponding statements for the normal 
form (|1U8|) . To go back to the variables of ()1U4[) mere substitutions 
are needed, which we omit. 

Most of Proposition El is a direct consequence of Theorems 1 and 
2. For the one-parameter family of solutions which are small in the 
right half plane we then have 




oo 

/i- J>- fe F fc (£(x)) (113) 

k=0 

As in the first example we find H^ by substituting 1)1 13j) in l|108j) . 
The equation of Hq is 

£ 2 H'o + £ H 'o = H o + 

The general solution of this equation are the Weierstrass elliptic func- 
tions of ln£, as expected from the general knowledge of the asymp- 
totic behavior of the Painleve solutions (see [H])- For our special 
initial condition, Hq analytic at zero and Hq{£) = £(1 + o(l)), the 
solution is a degenerate elliptic function, namely, 

m) = (£/J- l)^ 
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Next, one of the two free constants in the general solution H\ is de- 
termined by the condition of analyticity at zero of H± (this constant 
multiplies terms in ln£). It is interesting to note that the remain- 
ing constant is only determined in the next step, when solving the 
equation for H2I This pattern is typical (see ^6.9|) . Continuing this 
procedure we obtain successively: 



Hi = (216 £ + 210 e + 3 e - ^ e 4 J (£ - 12)- 3 (114) 

h 2 - (uwe + 52 38{2 - f e - ^ + ±? + 5^5) « - ur* 

(115) 

We omit the straightforward but quite lengthy inductive proof 
that all Hk are rational functions of £. The reason the calculation is 
tedious is that this property holds for (|1U8|) but not for its generic 
perturbations, and the last potential obstruction to rationality, suc- 
cessfully overcome by (jlU8|) . is at k = 6. On the positive side, these 
calculations are algorithmic and are very easy to carry out with the 
aid of a symbolic language program. 

In the same way as in Example 1 one can show that the cor- 
responding singularities of h are double poles: all the terms of the 
corresponding asymptotic expansion of 1/h are analytic near the sin- 
gularity of h\ All this is again straightforward, and lengthy because 
of the potential obstruction at k = 6. We prefer to rely on an existing 
direct proof, see [TS| . 

Let £ s correspond to a zero of 1/h. To leading order, £ s = 12, by 
Theorem |2] (hi) . To find the next order in the expansion of £ s one 
substitutes £ s = 12 + A/x + 0(x~ 2 ), to obtain 

„ , . (A - 109/10) 2 nl , . 
Vfr(e.) = - 12 3 X 2 +0(1/X S ) 

whence A = 109/10 (because 1/h is analytic at £ s ) and we have 

£ s = 12 + ^ + 0(*- 2 ) (116) 
lUx 

Given a solution h, its constant C in £ for which ()113j) holds can 
be calculated from asymptotic information in any direction above the 
real line by near least term truncation, namely 



C = Jim exp(^)x 1 /2 h{x) < 1 1 7) 

nrg(x)=cf> \ k<\x\ 
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Fig. 4. Poles of igQH for C = -12 (o) and C = 12 (+), 
light circles are on the second line of poles for to C = - 



calculated via l|116jl . The 
12. 



(this is a particular case of much more general formulas Jl]) where 

2fc>o ho,k%~ k is the common asymptotic series of all solutions of 
H1(J8|) which are small in the right half plane. 

□ 

General comments. 1. The expansion scales, x and x~ l / 2 e~ x 
are crucial. Only for this choice one obtains an expansion which is 
valid both in Strans and near poles of (|l()8|h For instance, the more 
general second scale x a e~ x introduces logarithmic singularities in Hj, 
except when a E — ^ + 7L. With these logarithmic terms, the two 
scale expansion would only be valid in an 0(1) region in x, what 
is sometimes called a "patch at infinity", instead of more than a 
sector. Also, a E — ^ — N introduces obligatory singularities at £ = 
precluding the validity of the expansion in Strans- The case a E 
— \ + N produces instead an expansion valid in Strans but not near 
poles. Indeed, the substitution h(x) = g(x)/x n , n E N has the effect 
of changing a to a + n in the normal form. This in turn amounts 
to restricting the analysis to a region far away from the poles, and 
then all Hj will be entire. In general it is useful thus to make (by 
substitutions in ©) a = a minimal compatible with the assumptions 
(al) and (a2), as this ensures the widest region of analysis. 

2. The pole structure can be explored beyond the first array, 
in much of the same way: For large £ induction shows that H n ~ 
Const n .^ n , suggesting a reexpansion for large £ in the form 
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h „ £ £ ^F iZ ; £2 = C^kx- 1 = CCWi- 3 / 2 e- (118) 

By the same techniques it can be shown that (|118|) holds and, by 

matching with 1)113(1 at £2 ~ x ~ 2 / 3 , we get H [ q ] = H with CM = 
— 1/60. Hence, if x s belongs to the first line of poles, i.e. £(x s ) = £ s 
cf. (|116|) . the second line of poles is given by the condition 

x~ 3/2 e~ xi = -60 • 12C 
i.e., it is situated at a logarithmic distance of the first one: 

x\ — x a = — liix s + (2n + l)m — ln(60) + o(l) 

(see Fig. 4). Similarly, on finds x Sj 3 and in general x s>n . The second 
scale for the n— th array is x~ n ~ 1 / 2 e~ x . 

The expansion (|113|) can be however matched directly to an adi- 
abatic invariant-like expansion valid throughout the sector where h 
has poles, similar to the one in [SO]- I n this language, the successive 
expansions of the form (|118j) pertain to the separatrix crossing region. 
We will not pursue this issue here. 



5.3. Example 3: The Painleve equation P2 
This equation reads: 

y " = 2y 3 + xy + a (119) 

(Incidentally, this example also shows that for a given equation dis- 
tinct solution manifolds associated to distinct asymptotic behaviors 
may lead to different normalizations.) After the change of variables 

x = (3t/2) 2/3 ; y(x) = x^{t h(t) - a) 
one obtains the normal form equation 

„ ti ( 24a 2 + 1\, 8,0 8a, 9 8(a 3 - a) , . 
h" + 1 + = — I h h H h 2 + — 5 — - = (120 

and 

Ax = 1, a x = -1/2; £ = ^; + & = F + ^F 3 

The initial condition is (always): Fq analytic at and Fq{0) = 1. This 
implies 
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Distinct normalizations (and sets of solutions) are provided by 

x = (At) 2 / 3 ; y(x) = (At) 1 / 3 (w(t) - B + JjL) 
if A 2 = -9/8, B 2 = -1/2. In this case, 



„ w' ( Wa l-6a 2 ^ 
w H \-w [ 1 + —. — — | w 



t V tA 9t 2 



w ~^a) w2 + w * + W + 6 ° 2) " rla{a2 ~ 4)) (121) 



so that 



implying 



1 3Ba 
Ai = 1, ai = — 

2 2i 



e 2 ^' + - F = 35F 2 - F 3 



and, with the same initial condition as above, we now have 

2j(l + BQ 

The first normalization applies for the manifold of solutions such 
that y ~ — — (for a = y is exponentially small and behaves like an 

Airy function) while the second one corresponds to y ~ — B — ^x~ 3 / 2 . 
6. Appendix 

6.1. Some results in classical asymptotics 

The notations and assumptions are those of 32 

Theorem 16 Let y(x) be a solution of (0) satisfying |3j) on a direc- 
tion d which is not an antistokes line. 

Let S be the open sector bounded by two consecutive antistokes 
lines which contains d. 

Then 

(i) for any d' C S the solution y(x) is analytic on d' for x large 
enough, and tends to along d' . Also 

(ii) 0) holds on d' . 
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These facts follow from the proof of Theorem 12.1 of (f° r more 
general contexts see also the proofs of |26| . |27j ) and from the proof 
(of a similar theorem) presented in jSJ]. Unfortunately, (i) and (ii) 
of Theorem ^] were not formulated in these references as results in 
their own right. They also follow from the more general results of 
jl.Sj . but their essence is of a classical asymptotics nature and the 
ideas of exponential asymptotics are not really needed. (To compare 
the results obtained using classical versus exponential asymptotics 
approaches see Theorem 1191 and the Remark following it, H6.2I ) We 
therefore include here a self-contained proof of Theorem ll6l The iter- 
ation argument of [SJ is set up as iterations of contractive operators 
on appropriate Banach spaces of analytic functions. 

Proof of (i) 

Setting. Fix rj > and let S v C S be the open subsector whose 
bounding directions form an angle r\ with the boundary of S. We 
assume r\ is small enough, so that d C S v . 

Let xq £ d and let D = xq + £„. It will be shown that if \xq\ is 
large enough, then the solution y(x) satisfying (jHJ) is analytic in D, 
and tends to as |x| — ► oo, x £ D. 

Denote yj(x ) = xp-. 

Note that for each index j, K(Ajx) has the same sign for all x £ D 
(there are no antistokes lines in S^). Divide the coordinates of y into 
the two sets I + , J_ 

J ± = {j = l, n ; M(\jx) £l± , x £ D} (122) 
Integral equations. Equation © can be written in the integral form 

Vj (x) = .,■'•■< A • .r" ( A [ .r, "■< X -' : g l (xr\y(xi)) dx x 

Jn^x) 

= ^ ] (x)+J 3 (y){x) , j = l,...,n (123) 

where the paths of integration Ilj(x) C D are: the segment [xo,x] if 
j £ /+ and the half- line from oo to x, along the direction of x — Xq 
for j G 2_ . 

Since the solution y(x) goes to along d we see that in (|123(l its 
constants of integration a,j are 

aj = for j E I- , aj = y j (xo)x" a i e x i x ° for j £ I + (124) 

By assumption g(a; _1 ,y) is analytic at (0,0), say for {x^ 1 < r 
and |y| < p2, and satisfies |g(a? — 1 , y) | < const (|x| -2 + |y| 2 ) (see ©. 

Let Bo be the Banach space of functions y(x) analytic on D and 
continuous on D (with the sup norm) . Let T be the closed subset of 
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functions y G Bo with ||y|| < p (where p > will be chosen small 
enough) and satisfying Uj(xo) = y® for j G I + and Uj{xo) = for 

Relations ()123j) . ()124|) can be viewed as an equation y = ip + J{y) 
on T (if \xq\ > r _1 and p < p-i). For p small we show that if y G T 
then if) + J{y) G T and the fact that J is a contraction on T . It will 
follow that the integral equation has a unique solution thus proving 

(0- 

Lemma. Let 

|xo| > max |3f?ttj| (|AJ sinr/) -1 max{l + \/2 , (v^sin??) -1 } 

j=l,...,n 

If : (i) 3 € 1+ an d x (t) = x o + — xo) 7 t £ [0,1], or 
(mJ j G 2_ and x(t) = x + i(x — xq) , t > 

1 x 



x(i) 



e £R[A*(*(t)-sO] < 1 



(125) 



The proof of this lemma is straightforward (the left side of l|125jl 
is increasing in t in case (i) and decreasing in case (h)). The following 
estimates can be used: | cos arg[Aj(x — Xo)]\ > s'mr] > 0, cos[arg(x — 
xo) — argxo] > — cos(2n) > —1, and for (i) R[Xj(x(t) — x)} = —(1 — 
t)\Xj(x — xo)\ cosarg[Aj(x — xo)] < —(1 — t)\Xj(x — xn)| sin?? , while 
for (ii) di[Xj(x(t) — x)] = t\ Xj(x — xq)\ cos arg[Aj(x — xq)] < — t\Xj(x — 
xo)| sin 77 

The set T is invariant under iterations 
Let y G T . For j G I + 



\i>j{x) +Jj(y)(x)\ < const \y 



X 
Xq 



jtlXjixo-x)] 



+ const \x — xq\ 





X 


f 

10 


x(t) 



JtMxW-x)] (^(t)!^ + || y ||2) dt 



and using ()125() and that |x(i)| > |xo|(l — cos 2 (2n)) the last term is 
bounded by 



const |y°| + const |x— xq| (|xq| 2 



+ l|y|| 2 ) f 

Jo 



e 2 



i(l-t)|A f (x-x )\siav fa 



< const \yA + const (|xo| 2 + ||y|| 2 ) 
Similar estimates hold for j 6 L, hence (for some K > 0) 

||^ + J(y)|| <^(|y°| + ko|- 2 + ||y|| 2 ) 
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Let p be small, such that p < (3K) . Then if |y°| < p(3K) , 
and |x |~ 2 < p{2>K)' 1 we have if> + J(y) £ T. 

Contraction. Let y and y' be in T . Writing gj^x" 1 , y) as gjfi(x~ 1 )+ 
^2k=i,...,n9j,k( x ~ 1 ^y)yk witn 9j,k analytic for < r and |y| < p 2 

and gj^k = 0(x~ 2 ) + 0(|y|) for k > 1, and gjp = 0(x~ 2 ) then 
\gj{x~ ,y) — 9j(x~ , y')| < const(|x| -2 + p)|y — y'| so that, with esti- 
mates similar to the above, we get || J{y) — J{y')\\ < const(|xo|~ 2 + 
p)\y — y'|. For small p and large xq the operator J" is a contraction 
on J 7 , and part (i) of Theorem is proved. 

Remark. In the estimates above the smaller 77 (i.e. the closer x to 
an antistokes line) the larger xq must be. This is closely related to the 
fact (which is the object of the present paper) that solutions (which 
are analytic in a "sector" — more precisely, in a region described in 
Theorem lltjf) develop (generically) singularities on the edges of this 
"sector" . 

Proof of (ii) 

Let 4>{x) be a solution of © satisfying Q — which is known to 
exist Let y(x) be a solution satisfying (jHJ). Let u(x) = y(x) — 
4>(x). It is enough to show that u(x) = 0(\x\~ r ) for all r G NU {0}. 

The function u(x) has limit along d and satisfies 

u' = -Au + iiu + h(x _1 ,u) (126) 

x 

where 

h (aT 1 , u) = g (aT 1 , u + <f>(x)) - g (aT 1 , <ft(x)) (127) 
As in the proof of (i) we write (|126|) in integral form (similar to 

Uj-(x) = r A 'a ; r A a / .r ; ' V V > ; (a^ 1 , U^)) dx x 

Jnj(x) 

= ^j(x) + Jj(u)(x) , j = l,...,n (128) 

where the paths of integration are those of (|12cS|) and the constants 
dj satisfy the analogue of (|124f) 

o,- = for j € I_ , Oj- = Uj(xo)x~ aj e X J xo for j G 7+ (129) 
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where u = y(x ) - 4>(x ). 

Let D be as in (i), where cj) is analytic. Consider the Banach space 
B r of functions u(x) analytic on D, continuous on D, with the norm 
||u|| = sup X £ D |x r u(x)| (see also (50]). 

Let T be the closed subset of functions u 6 B r with ||u|| < p 
(where p > will be chosen small enough) and satisfying u(xo) = Un. 

Note that <j>(x) = 0(x" 2 ) (cf. ©) hence \</)(x)\ < M\x\~ 2 (for 
x £ D and xo large enough). Then since g(x _1 ,y) was assumed 
0{x~ 2 ) + 0(|y| 2 ) it follows that for x £ D (cf. (tTTTlO 

|h (a; - , u) | < const (|x|~ 2 |u| + |u| 2 ) 

(for zo large enough, so that |</>(x)| < p/2 and for |u| < p/2). 

The same estimates as in the proof of (i) (the only difference 
being that the otj of (i) should be replaced here by r + aj) show 
that equation u = ifi + ^7(u) has a unique solution in T if |xo| is 
large enough (depending on r, as expected). Hence ||u|| r < oo which 
concludes the proof of (ii). 



6. 2. Summary of some results in liffi 

This subsection contains details on results of ^3] cited, referred to, 
or relevant for the present paper. A simple consequence of a Lemma 
in JH] (needed for the present paper) is formulated and proved at the 
end of this section (Theorem IT§|) . 

Since the Theorems, Lemmas and some formulas cited in this sec- 
tion are from ^3]) to avoid repetition we will follow by a * sign any 
result cited from [T3] . 

The setting is the same as in the present paper: the equation 
studied is 

y , = -yiy-- J By + g( 2 r 1 ,y) (130) 

x 

(same as with B = —A) having transseries solutions 

y(s) = C k e~ A - ka; x-^ k s k (x) 

ke(Nu{o}) n 

^ £ C k e-^* M - k yk (z) (131) 

ke(NU{0}) n 

(same as ©, ©, (HU for 0j = -aj, Mj = - + 1, j = 1, ...,n). 
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Since the association between actual and formal solutions depends 
on directions (Stokes phenomena) a sector in the complex a>plane 
is chosen as follows. Fix some non-empty open sector S' C C and 
consider those transseries (j!31[) valid in S' (as explained in ^2.1(1 . 
Some constants Ci, C n may be required to be zero in S', say Cj = 
for j = rti+1, n (with m > 0). Let Strans be the (non-empty, open) 
maximal sector of validity of any transseries (jlHlj) with Cj = for 
j = n\ + 1, n (see ©). 

To simplify the notations it can be assumed (after trivial changes 
of coordinates) that Ai = 1 (see also t J2.2j) . 

I. The construction of actual solutions associated to transseries 
solutions valid in Strans is done in ^H] using a generalized Borel 
summation as follows. 

Denote by Y(p) the formal inverse Laplace transform of y(x) (i.e. 
Y(p) = J°+. l0 ° e px y(x) dx, its convergence following from sub- 

sequent analysis). Using the usual properties of the inverse Laplace 
transform (e.g. the transform of y'(x) is —pY(p), multiplication is 
transformed into convolution, etc.) the differential equation (|13Uj) is 
transformed into a convolution equation (eq. (1.13)*). 

The Stokes lines in the Borel plane (p-plane) are defined as the 
complex conjugates of the Stokes lines in the direct (x) space. For 
linear equations the Stokes lines in the p-plane are dj t o = AjR+, 
j = 1, n. For nonlinear equations, there also are other Stokes lines 
(which play a role only in the higher order terms of the transseries, 
with |k| > 2, hence they are "transparent" in the linear case) namely 
d jik = (Xj - k • A)R+ with j = 1, ...,n, k G (N U {0}) n (note that 
Pj,k G d j,k of.®). 

Once the sector Strans is fixed, there are only finitely many lines 
djfc in this sector. 

In any proper subsector of any of the n sectors formed by the 
Stokes lines dj t o the convolution equation has a unique solution Yo(p) 
which is analytic at p = (Lemma 16*). Yo(p) is in fact the Borel 
transform of the asymptotic series yo(^) (see (JHJ) (i.e. 
Yo(p) = Ylr yo,rP r /^0- Yo has singularities on the Stokes lines 
at p £ AjZ + , j = l,...,n (hence the classical Laplace transform 
cannot be taken on R + and the classical Borel sum of yo does not 
exist). 

Denote by Yq the analytic continuation of Yo on directions above 
d\ t o = M + (but below the neighboring Stokes line), respectively by 
Yq for the continuation below R + ; they exist see Lemma 16*. It is 
shown that as p approaches R+ from above (or below) Yq [p) (respec- 
tively, Yq(p)) tends to a distribution on R + in an adequate space of 
distributions — the staircase distributions, introduced in (Lemma 
16*). In general, the two distributions are different. (Of course, a 
similar picture holds at any other Stokes line.) 
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Higher order functions Yk, |k| > 1 are then constructed (Lemma 
20*) by solving the convolution equation on the Stokes lines. 

Consider for example the line M+ = di,o- Fix a solution Yo of the 
convolution equation in the space of staircase distributions on K+. 
The construction of the higher order Yk's with kj = if j > n\ + l 10 
is done in the proof of Lemma 20* as follows. In view of the sought- 
for expansion (|10[). after introducing it in H130j) and identifying the 
coefficients one obtains (a recursive system of) differential equations 
for y^; formal inverse Laplace transform yields (a recursive system of) 
convolution equations for Yk- Once a staircase distribution solution 
Yo on R + is chosen the general solution of this system with regularity 
(|34"|) depends on n\ free constants C±, C ni , in the form C k Yk- 
Outside the Stokes line M+ the solutions Yk are, in fact, analytic 
up to the nearest direction (of positive argument ip + , respectively 
negative argument tp-) which is either a Stokes line d,k which lies in 
the right half-plane (i.e. the half-plane orthogonal to d\ t o = M+) — 
where some of the constants Cj may change, or is an antistokes line 
associated to Ai = 1, i.e. iM+ or iM_ — where the transseries is no 
longer defined 11 — (Lemma 20* (i)-(iv)). 

It is interesting to note that the Yk multiplied by Stokes con- 
stants are generated as differences between different branches of Yo 
(Theorem 4*, Proposition 23*). 

Then (a generalized) Laplace transform is applied to Yk (in the 
space of staircase distributions on the Stokes line under consid- 
eration) yielding (analytic) functions yk = £Yk- 

The last step in the summation of transseries is showing that the 
sum IjlOjl converges (Lemma 20* (v); more details are found in the 
proof of Lemma IT71 of this section). 

The reconstruction of a solution from a transseries is concluded 
showing that the function y(x) obtained as the sum of (jlOj) is a solu- 
tion of the differential equation (jl.SOj) — which follows easily because of 
appropriate convergence and since all functions have been constructed 
from formal objects satisfying the equation (Lemma 20* (v)). 

The correspondence between transseries (i.e. the constants C) and 
actual solutions given by the summation of Lemma 20* is not unique. 
This is due to the non-uniqueness of staircase distribution solutions 
Yo on R + (on which the higher order Yk depend): there is a one 
parameter family of solutions Yq (among which Y^ for a = 0, 1) — 
they are special averages of analytic continuations of the germ of 
analytic function Yq at p = 0. 



Other Yk's are not needed since the corresponding cannot be present in a 
transseries on the fixed Strans- 

11 Yk may be analytic beyond the antistokes lines ±iR+ but the analysis in |13| 
stops there. 
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2. Conversely, any solution of (|130|) satisfying on a direction 
d in the right half-plane is the sum of a transseries, which is unique 
once q is fixed (Theorem 3*(iii)). 

The proof is done in the following steps (Propositions 24* and 
25*). There exists a solution yo a (x) asymptotic to the transseries 
yo (i.e. y with all Cj = 0) as constructed at part 1. above. It is 
also shown that any two solutions y 1,2 (x) with the same asymptotic 
expansion y 1,2 ~ yo on d differ by exponentially small terms: y 1 (x) — 
y 2 (x) = J2j=i n Cje~ x i x x~ l3 3 (ej + o(l)) on d. Thus the difference 
y(x)— yg (x) fixes the constants Cj = Cj(a, d), which are then used to 
construct a solution yt(x) from the transseries y with these constants 
(using part 1. above). The last step is showing that y(x) = yt(x). 

3. It has been thus established that (given a direction d, and a 
parameter a) there is a one-to-one correspondence between trans- 
series solutions and actual solutions of (|130j) . The correspondence is 
built using a (family of) generalized Borel summation(s) CB a on d. 
It is shown that the operator CB a is compatible with all algebraic 
operations (performed on transseries, respectively functions on d). 

The Stokes phenomenon is analyzed in Theorems 4* and 5*. 

4. We now state the maximal domain of analyticity 12 of y(x) 
implied by formula (2.41)* of Lemma 20*. 

Lemma 17 Let y(x) be a solution of U3(J\) satisfying on a di- 
rection d above (but close enough to) M + . 

Let CJ , j = 1, ...,n be the constants such that y(x) is represented 
on d as 

y = £Y + ^2 (C-) k e- k Aa; x M - k £Y k - (132) 

[k|>0 

(see Theorem 3* (Hi)). 

Then for any e, 5 > there is x\ > such that y(x) is analytic 
on the domain 

— ( _ TT _ TT 

D an = 1 x ; M > x i , arg{x) G [-tp -- + e,-ip + - - e] , 

<,5- 1 , j = l,...,n} (133) 

The constants e, 5 are the same for all solutions of \13CJ\) with 
transseries valid in the same sector St ran s as y- 

The Proof follows immediately from Lemma 20* (v), but we pro- 
vide the details. 

All functions Yi c ~(-e 4(?i ) are Laplace transformable in the space of 
staircase distributions on M + with exponential weight e~ vp (for v > 
large enough) if <f> G (ip_,ip + ) (see Lemma 20* (v)). 

12 This domain is improperly stated in 
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For <j) £ (ip-,0) the Yk _ (pe 1 ^), (p > 0) are analytic (Lemma 
20* (iii) ) , so the Laplace transform in the space of staircase distri- 
butions coincides with the classical Laplace transform (see Lemma 
6* and the classical properties of the Laplace transform and of the 
spaces L\^) where the Laplace transform is defined as 

CF{x) = [ e~ px F(p) dp , for x £ d (134) 
Jd 

(Note that in (fr31|) p £ d and x £ d so that 5i(px) > 0.) 

Remark 18 If F is analytic on a direction d = e~ tv M + (rj £ ~S.) and 

\\ F \\u = [ e~ up \F(p)\ dp < oo 
Jd 

then its Laplace transform \1'34\) ^ s analytic for x £ d = e ir, M + , |x| > 
v . 

Furthermore, for any e > CF has analytic continuation on the 
domain argx £ [i] — ^ + e,rj + ^ — e] = I and \x\ > ^(sine) -1 and 
CF satisfies 

\CF(x)\ < \\F\\ U (135) 

The proof is immediate, noting that the path of integration of 
()134|) can be rotated to any other direction in the interval / — since 
/ is so that $i(xp) > 0. 

To prove Lemma^llet e be small, positive and let eft £ [tp- +e, — e]. 
For any 5 > there is v > so that the functions (pe^), (p > 0) 
satisfy ||Yi c ~(-e l<?!> )|| i y < <5' k ' for all multi-indices k considered (i.e. 
|k| > 0, kj = for j > n\) see Proposition 22* (ii). Using Remark UHl 
and (|132j) the result is immediate. □ 

Remark 

Similarly to Lemma El there is x\ > such that y(x) is analytic 
on D+ n = D+ n (e,5,C+) where 

D tn = \ x \ \x\> xi, arg(x) £ [-ip + - - + e, -ip + + - - e] and 

<6- 1 j = l,...,n] (136) 

so that the domain of analyticity of a solution y(x) includes domains 
of the form D~ n U D+ n . 

Theorem 19 Let y(x) be a solution of hl'JU]) satisfying (0) on d = 
M + . Let e > be small. 

There exists 5, R > such that y(x) is analytic (at least) on 

San = S an (y(x); e) = S+ U S~ (137) 
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where 



Sf = ^x;\x\>R, arg(x) G [— — =F e, — =F e] and 

< r 1 for j = l,...,n} (138) 



The constant 5 is the same for all solutions of \13U\) with trans- 
series valid in the same sector Strans as y- (However, R does depend 
on C.J 



Proof 

From Lemma IT71 using the expansion of Y in terms of Yk~, re- 
spectively Yk + it follows that y(x) is analytic for \x\ > x±, argx G 



r+e,— ip +f — e] = / and 



Cfe- 



X i x X 



Mi 



< 5- 



• J 



1, 



■ n. 



Since Mj = — l$t(3j\ + 1 and the sector / is larger than the sector 
where all exponentials in the transseries of y(x) are bounded, the 
result follows. 



Remark 

Fix a direction d (not an antistokes line) and consider solutions 
satisfying (JHJ). It is interesting to compare the result on the domain 
of analyticity of these solutions as given by Theorem El (obtained 
using results of exponential asymptotics) to the result of Theorem 
(obtained in a classical setting). There are special families of 
solutions for which the sector of analyticity given by Theorem 
is, in fact, larger than the sector between two consecutive antistokes 
lines (i.e. solutions having the corresponding C,- zero). For Painleve 
PI equation these special solutions are called triply truncated. 

Convention 

The Borel summation used in the present paper is CB = CBi . 



6.3. Proof of {££p 

Proof. This follows from the definition £ = C\x ai e~ x and from the 
asymptotic behavior of the functional inverse W of se s (see e.g. |12|). 
For large t > the branch of W which is real has the expansion 
(convergent, as it is not difficult to show) 



ttt, xi i . lnln * A(lnlnt) 2 -lnlnt 

W (t) = hit- In lilt + + ^ TT-rr= + 

mi (lntj^ 
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6.4- Points on 7^ 



have the same magnitude 



Let R be large, so that 



P3 + 



3KB 
R 



< P2 



(139) 



Lemma 20 There is a small enough neighborhood MjO 0/7^ so that 
any x' , x" E A/lo satisfy 



7N(to)\ - \j N (t)\ = 9 (ln(r(t)) - ln(r(t ))) + o(l) (N -> 00) 



Therefore lirri2v->oo TaKO/t^KO = 1 uniformly for G [£o>l]- So 
for TV large \x/x'\ < 3/2 for all x, x' on 7^ between xq and a, so in a 
small enough neighborhood |x/a;'| < 1/2 which proves the Lemma. 

6.5. Special estimates 

We show that the second argument of g in (|65|) has absolute value 
less than p 2 (cf. (J23J)). 

We need the following lemma. 

Lemma 21 Let A, en > 0. There exist Aq, k > sue/t i/iai 



if A > Aq and m\A~ m < c . 

Before giving the proof of the lemma, we show how (j!41|) is used. 
Let R be large, and cq small, so that 




(140) 



Proof 

Using Q24|) for t G [to, 1] we get the uniform estimate 




(141) 




(142) 



In view of Theorem for s G A/lo 
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and in view of Lemma QUI the last term is bounded by 

fc=i v 1 1 ' 



Using 1)143|) , Lemma |^ and the bound p% on Fo and ()142|) we 

have 

|FH ( X )| < |F (*(*)) \ + K K ( 2 ^ + c ) < < P2 



Finally, if \S(x)\ < on Jfy then 



Proof of Lemma 1211 

Estimates like (|141j) are common in proofs using least term trun- 
cation of factorially divergent series (see e.g. Jl]). The proof of ()141j) 
is included here for completeness. 

The series 

k>l 

is divergent. Its terms decrease for k < A and increase for k > A; the 
term with k = [A\ (or k = [A\ + 1) is called the least term (see e.g. 
0)- 

Case I: m < A 

In this case the terms in the l.h.s. of (|141|) are decreasing, hence 
kl 1 2 , ,3 

Ei^i+i2(--i)<i 

fc=l 

Case //: A < m < eA/2 

The terms in the l.h.s. of 1)141)1 are increasing for A < k < m, not 
exceeding the second term: m\/A m < 21/ A 2 . 

Indeed, this is a simple estimate using Stirling's formula and the 
fact that the function F(A) = A^ 2 a aeA is decreasing for A > Aq (if 
Aq is large enough). 

Then as in Case I 

^ k\ 1 2 . . e + 1 

Y.T^A + l¥ {m - l)< — 

k=l 



Case III: eA/2 < m 

Denote p = [A\ and q = [^-p - \\ . 



Formation of singularities in nonlinear differential systems 



57 



Write 



fc=i 



where 



Jfc! 



Sin + Sn + Si 



k\ 



k\ 



Sni = ^2jk ' Yl -jk > Si= E # 



k=l 



k=p+l 



A k 



k=q+l 



A k 



and estimate each sum separately. 
To estimate Si 



S < ^ < m ' 

k=q+l F F 



A;=g+1 



(5 + 2) m ~ fc 



and since + 2) < 2/(^4 + 1) < 1 this is less than 

m! A + 1 m! e + 2 
p"*A - 1 < p™e - 2 

To estimate 5// note that 



£// < - p) -7- < -r-(<? - p) ( , 

' Ai A mK 1 \ q + 1 



q! to! 



and since < po < 1 if A > Aq (for Aq large enough) 



< j^{q-p)pq 9 < '—r>\. 



ml 



Finally 

The result of Lemma |2*T1 follows. 



SjII< J + ^ (p ~ 1)< J 



6. 6. Proof of Proposition ^] 

A consolidation of one of the norms in Example (3a) in ^3] is first 
needed. For convenience we repeat that part. The notations are those 
in pig. _ 

"(3a) For > and fa ^ fa, let T^UV) = {/ : f(p) = 

p@F(p)}, where F is analytic in the interior of £ U V and continuous 
in its closure. We use the family of (equivalent) norms 
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u ,p = \r(P + l)\K sup |e-^/(p)| 
sesuv 



(145) 



It is clear that convergence of / in \\\\ v fi implies uniform convergence 
of F on compact sets in £■ U V (for p near zero, this follows from 
Cauchy's formula). Tp are thus Banach spaces and focusing spaces 
in || 11^/3 by (j!45j) . The spaces {Tp}p are isomorphic to each-other. 
Convolution is defined as 



~01~02 



~\fi * h)(p) = C t^Ftipt^l - tf 2 F 2 (p(l - t))dt = F(p) 
Jo 

(146) 



where F is manifestly analytic, and the application 
(•*•): Tp 1 xT ft « T Pl+ p 2+ i 

is continuous: 



(147) 



||/l * f2\\u,0 1 +/3 2 +l 



\K sup 


z~ up [ 


V 


Jo 



s Pl Fi(s){p - sf 2 F 2 {p - s)ds 



r(A + fa + 2) 

< =z sup 



A'F 1 (s)e- I/s s /31 KF 2 (p - s)e- u( - p ~ s \p - s] 



02 



r(Pi + 1) 



d|s| 



m + 1) 

^WhUpAhWv^ (148) 



Estimating the norm of Y k exactly as in ^3] but using this inequality 
instead of (2.8) of ^3] we get that 

|Y k (p)| < |r(-k-a / ) _1 |^' k| e 1/obl 

and thus, using straightforward Cauchy estimates in d ai of deriva- 
tives, is proved. □ 
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6. 7. Note on normalization of the on 

The reference |13| uses a transformation that makes f3j < (ay > 
in the present notation). To determine the singularities of y it is now 
important to make ctj as small as possible, as explained in Comment 
1, N5.21 In some cases we must then allow for m < in |13| . Eq. 
(2.43). This does not affect the estimates (2.44) through (2.46) in the 
space 7k.pi -i, the only one that relevant to the present paper. Minor 
modifications of the proof following Lemma 20 in ^3] are needed. For 
completeness we redo redo here the whole proof. 

For |k| > 1 with Wk := Yk and Rk := Tk, the functions Wk 
satisfy the equations 

(1 + Jk)W k = Q^Rk (149) 

with Qk := ( — A + P + k • A) (notice that for |k| > 1 and p E S' we 
have detQk(p) 7^ 0). 



(J k W)(p) := Q,, 1 ((s + m-k) jTw(«)ds 



v 



£ / W^-Djip - s)ds ] (150) 



o 



Proposition 22 i) For large v and constants K\ and K.2{v) inde- 
pendent o/k, with K^iy) = 0{y~ l ) we have \\Q%. || < t& and 

|| J k || < K 2 H (151) 

iij For Zarge v, the operators (1 + Jk) defined in 7>' mv , and also in 
7 ]s _pi_ 1 for |k| > 1 and in 7\ for |k| = 1 are simultaneously invertible. 
Given Yo and C, the Wk, |k| > 1 are uniquely determined. For any 
5 > there is a large enough v, so that 

||W k || < 5 |k| , k = 0,1,.. (152) 

(in the V mv topology, \15ty) hold uniformly in (j) G [ip- + e,0] and 
(f> 6 [0, -0+ — e] /or any smaZZ e > 0). 

□ 

Proof. 

(i) follows immediately from ()148|) . 

(ii) From Q149|) and (i) we get, for some K and j > 1 ||Wk|| < 
if||Rk||- We first show inductively that the Wk are bounded. Choos- 
ing a suitably large v{e) we can make max| k |<i ||Wk|| y < e for any 
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positive e (uniformly in (p). We show by induction that ||Wk||^ < e 
for all k. In the same way as in |13j we get 

HWkii, < ^iiR k |i, < y, k i^ 1 1 

l<k (i mp ) 
|k| 

< el k l ^2 ^2 ni d k l +s )2 s+ni < (Cie)l k l (153) 

s=0 

where C\ does not depend on e, k. Choosing e so that e < C-f 2 we 
have, for |k| > 2 (Cie)l k l < e completing the induction step. But as 
we now know that HW^H^ < e, the same inequalities (j!53j) show that 
in fact HWklli/ < (Cie)' k L Choosing e small enough, the first part of 
Proposition 1221 (ii) follows. □ 

6. 8. Proof of Lemma MIA (ii) 

1. Generically ho is not entire. Assume /in is analytic in a neighbor- 
hood of the disk Br = {|£| < R} and let Mr = sup| z | =R /in(z). We 
have 

sup {z-\Hq{z) - 1)} = sup {aT 1 (/*,(*) - 1)} < R- 1 {M + 1) 

\z\<R \z\=R 
and thus, for some constants Cj we have 

M 2 < CiR 2 + C 2 + {C 3 R 2 + C 4 )M (154) 

whence 



M < C 5 R 2 + C 6 (155) 

If /in is entire it then follows that /in is a quadratic polynomial in 
£. But it is straightforward to check that does not, generically, 
admit quadratic solutions. Thus the radius of analyticity of /in is 
finite, say R , and 13 

sup {z-\ho{z) - 1)} < R \M Ro + 1) and M Ro < C 5 R 2 + C 6 

\z\<R 

(156) 

13 An upper bound for Ro can be found by comparing /io"(0) with its estimate 
from Cauchy's formula and l|154|l 
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2. /iq is uniformly continuous on Br . Indeed, if € £>_r we have 
by (jHSJ) and (11561 that /i 2 , is in fact Lipschitz: 

1^(0-^(01 < Connie (157) 

3. If £o £ 5^i? an d ^o(Cq) 7^ then £o is a regular point of eq. (|51)l 
and thus ho is analytic at £o- 

4. J/£ s G 0Br is a singular point of ho and — A^.7 1 + ^3 + d^ s ^ 
then £ s is a square root branch point of ho, i.e. ho(£) = — 4s) 1 / 2 ) 
where h\ is analytic at zero. From parts 2 and 3 above, ho(£ s ) = 0. 
It is convenient to look at the equation for £(/to) derived from ((Ml) : 

= ^ . m = t 

dh iM~ 1 + d 1 + d 2 £)ho + (-As^- 1 + d 3 + d 4 

(158) 

whose unique solution is analytic near zero. The claim now follows 
by noting that f'(0) = and £"(0) = 2(-A 2 £: 1 + d ?J + d£ s ). 

5. We now restrict the analysis to a smaller but generic set of coeffi- 
cients. We denote by K s the following subset of parameters (see (|79j) 
and (|HH) ) 



K s = {(d,j) = (dj, 7j)j=i,..., n G C 2n : ho not entire and 

P(£) = -2A 2 + A 2 di£ + (f 3 4 2 has distinct roots } (159) 

We show in parts 6 through 8 that if £ s G 9Bji is a singular point 
of ho and -P(£ s ) = 0, then a generic small variation of (0,7) in i^T s 
makes P(£ s ) 0) and by part 4, £ s becomes a square root branch 
point of ho- 

We thus assume that P(£ s ) = 0. The substitution £ — £ s = i in 
(JBU) gives 

2/jo/i'o = B x h + tB 2 (160) 

where 



Bi = 2\ 2 {t + 6)" 1 - d x + (t + 6)c/ 2 ; # 2 = *tf% g rf 3 

and the roots of P(£) are £1 7^ £ s . 

We now study /io in the following smaller region. Let to be small 
on the segment [0, — £ s ]. Choose Z\ = {t : \t — to\ < \to\} C Br q — £ s 
to be a disk tangent at t = to i?/j — £ s which does not contain 
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the points £ s — £ and — £ s . Then /io is analytic in Z\ and continuous 

in Z\\{0) (while hn is continuous in Z\), and lim hn(t) = 0. We 

A3t-*0 

assume t = is a singularity of ho- 

6. There exists a sequence {t n } n in A with t n — > suc/i that 
lim^oo h (t n )/t n = L with 2L 2 - LBi(O) - B 2 (0) = 0. 

We first show that ho — * 0, then prove ho/t is bounded below and 
above, and finally that ho/t has a limit. 

(a) . We estimate 

.Mt = max |/in(s)l 
l--*l<l*l 

for t £ (0, to) from (|16U|) written as 

/»§(*) = f B 3 (s)h (s)ds + t 2 B 3 (t) 
Jo 

where B 3 ^ are analytic in A and continuous on A. Then 

M 2 < \t\M t max |B 3 | + \t 2 \ max |B 4 | 
A A 

and thus Ait < | i | , for some ifi > and all t £ [0, to]. 

(b) Cauchy's formula on the circle \s — t\ = \t\ implies |/i (t)| < K\. 

(c) Equation IjlbUj) written in the form t/ho = (2h' — B^B^ 1 
implies now \t/ho\ < < oo. In conclusion, 



t 

h 



£ [K 2 ,Ki] for all t E [0,t ] 



To conclude the proof of the statement at the beginning of this part, 
the function y = ho/t, which is analytic in A and continuous on 
A \ {0} satisfies the equation 

2ty' = -2y + B x + B 2 /y 

which can be written as 

i y' 



2t P (y)+tf(t,y) 



(161) 



with Po(y) = —2y + Bi(0) + B 2 (0)/y. Assume, to get a contradiction, 
that for some e > we had |P (y(*))l > e for t £ (0,ti] C (0,t }. 
Since 



+ th(t,y) (162) 



Po(y) + tf(t,y) P (y) 



Formation of singularities in nonlinear differential systems 



63 



with fi(t,y) bounded on (0, t{\ if t\ is small, we get by integrating 
(Q2TJ) on [t,ti] C (0,ii] 

i ln(iAi) = FoMt)) - Fo(y(t x )) + jf ai/00/l(*, 2/00)^ 

where i*o is a primitive of l/Po> anci thus Fo(y(t)) is bounded for £ E 
(0,ii). Hence the r.h.s. of (j!62|) is uniformly bounded for t E (0, t±), 
which is a contradiction, given the l.h.s. of ([162)) . 
7. With S(t) defined by h (t) = Lt(l + <5(t)) we have S(t) -> as 
i — > 0, i E (0, to]. We show this by a contraction argument. 5 satisfies 
the equation 



tf'= (-2+^-\s + B(t,S) (163) 



where 



#(M) 



gi(0) 
2L 



1 + 5 



+ t 



Bi(t) 



2Lf 



B 1 (0) + B 2 (t) 



B 2 (oy 



2LH{\ + 5) 



or in integral form, 



6(t) = (JS)(t) = (l-Y S(t n )+t b ' t s-^BfaSisfids (164) 

\ 71 / J t n 

where {t n } n is the sequence found in 6, such that S(t n ) — ► 0. We 
have two cases, according to whether is positive or negative 

(3i(£>i) = is nongeneric). 

(a) In the case M(bi) > we define J on the space A of analytic 
functions in B n and continuous in B n , where B n is the ball having 
[0,t n ] as a diameter. Let A r = {5 E A : ||<5||oo < r} and r n = 2|5(t n )|. 

(b) For 3?(&i) < we instead define J on the space A' of analytic 
functions in B' n and continuous in B' n , where B' n is the ball having 
[t n ,t n -%] as a diameter. Let „4.' r = {5 E .4' : ||£||oo — r } an d r n = 
2\5{t n )\+2(t n _ l -t n ). 

Lemma 23 ^aj If $t{b{) > 0, /or n /arye enough, J : A Tn i— ► A Tn 
is contractive. Therefore, as n is large, we have \5{t)\ < 2\5(t n )\ on 
[0,t n ] so that 8(t) ast^O in (0,t ]. 

(b) 7/5ft(6i) < 0, forn large enough, J : A', i— > A', is contractive. 
Therefore, as n is large, we have \8(t)\ < 2\5(t n )\ + £ n _i — t n on 
[t n ,t n —i] so that, again, 5(t) — ► as t — ► in (0, to]- 
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Proof. A straightforward calculation. 

8. Now we bootstrap the information that 6(t) —* to sharpen the 
characterization of 5(t) for small t. 

Lemma 24 (a) If dt(bi) < then 5(t) is analytic in t at zero, and 
<5(0) = 0, thus 5{t) = 0{t). 

(b) If Sft(6i) > then 5(t) = 0(t bl ) + 0(t) for small t. 

Proof, (a) Since 6(t) — > we have 



S(t)=J (tzb 3 (S(tz),tz) +b 2 5 2 {zt)^z- bl - l dz (165) 

which for small t is manifestly contractive in a small sup ball in a 
space of analytic functions in a neighborhood of t = 0. 

(b) Fixing some n, from (|164j) and Lemma l28l we see that 

W)\< l -^1* + (t n - *)*P[|oo 

□ 

(In fact it is not difficult to show that in case (b), 5 can be written 
as an analytic function in the two variables t and t bl .) 
With u = we get the equation in variations 



uti + h u' = (2A 2 r X + a 2 )u + 1 (166) 

whence u is analytic as long as ho ^ 0. This equation being linear we 
can use classical Frobenius theory and, in the only interesting case 
R(6i) > we have u(£ ) = (Ai j2 - 2X 2 ^ 1 - aa)" 1 ^ 0. Thus an 
arbitrarily small variation of A 2 makes /io(£o) 7^ 0. 



6. 9. The recursive system for F m 

In applications it is usually more convenient to determine the func- 
tions F m recursively, from their differential equation. Formally the 
calculation is the following. 

The series F = J2m>o x ~~ m ^m{€) is a formal solution of @; substi- 
tution in the equation and identification of coefficients of x~ m yields 
the recursive system (|56|). (|57|). To determine the F m 's associated to 
y we first note that these functions are analytic at £ = (cf. Theo- 
rem^). Denoting by F m j, j = 1, .., n the components of F m , a simple 
calculation shows that (J56|) has a unique analytic solution satisfying 
F ,i(0 = £ + O(e) and F 0>j (0 = O^ 2 ) for j = 2, For m = 1, 
there is a one parameter family of solutions of (|57|) having a Taylor 
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series at £ = 0, and they have the form -Fi,i(£) = ci£ + 0(£ 2 ) and 
Fij(£) = 0(£ 2 ) for j = 2, ...,n. The parameter c\ is determined from 
the condition that (|57[) has an analytic solution for m = 2. For this 
value of ci there is a one parameter family of solutions F2 analytic at 
£ = and this new parameter is determined by analyzing the equa- 
tion of F3. The procedure can be continued to any order in m, in the 
same way; in particular, the constant c m is only determined at step 
m + 1 from the condition of analyticity of F m+ i . 



6.10. Sketch of a classical proof of Theorem^ 

It is also interesting to mention a direct, classical proof (i.e. not in- 
volving results of exponential asymptotics) of Theorem ^ (Since we 
do not rely on this more involved approach, we only give a brief out- 
line of this proof.) 

Having determined the initial conditions for F m as above, equa- 
tions (|56[). (|57|) can be transformed to integral equations possessing 
unique analytic solutions F m for small £. 

To show (|2T)|) let y(x) be a solution of (j2J such that y(x) — » in 
Strans- Denote 

/ N 
R N (x) = x N+1 y{x) - x- m F m (C(x)) 

V m=0 

(the remainder of y(x) with respect to the truncated expansion). 
Then R^r satisfies the differential equation 

R N = B N (x,R N ) (167) 
where 

B N (x, Rjy) = d y g(0, F )R N + F(x)R N + u N {x) + F^ N {x, R N ) 

(168) 

with F^ N {x,R N ) = 0(R 2 N ). 

Let RnJi j = 1, denote the components of Rjy and Enj be 
the components of E^r- 

Equation (|167|) can be written in the integral form 



art 

dx 



N 



+ 



A + x~ l (A - (N + 1)1) 



R NJ (x) = e -* x x N+1 - a i / e^ s s- N - 1+ ^E N (s,R N (s))ds (169) 

J 

(for j = l,...,n). For an appropriate (rather delicate) choice of the 
initial points x® and of the contours of integration in ()169j) the in- 
tegral operators defined by the r.h.s. of (|169|) are contractive for N 
sufficiently large, hence (j!69|) has a unique analytic solution R^v- 
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